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CHAPTER  I 


EXECUTIVE  SUMMARY 


1.1  INTRODUCTION 

Cementitious  materials  such  as  mortar  and  concrete  are  brittle  and 
have  an  inherent  weakness  in  resisting  tension.  They  are  known  to  crack 
under  low  levels  of  tensile  strains.  The  addition  of  discontinuous  fibers  to 
such  matrices  leads  to  a  drastic  improvement  in  their  toughness.  It  is 
generally  agreed  that  the  fiber  contribute  primarily  to  the  post-cracking 
response  of  the  composite  by  bridging  the  cracks  and  p'oviding  resistance 
to  crack  opening.  The  transmission  of  forces  between  the  fibers  and  the 
matrix  is  achieved  through  interfacial  bond,  defined  as  the  shearing  stress 
at  the  interface  between  the  fiber  and  the  surrounding  matrix.  Bond  has 
been  recognized  as  a  major  factor  in  composite  action.  There  is  hardly  any 
property  of  the  composite  that  is  not  dependent  on  bond.  Thus  there  exists 
a  genuine  need  to  provide  a  fundamental  undestanding  of  bond  and  bond 
mechanims  in  fiber  reinforced  cement  composites. 

The  nature  of  bond  in  today's  fiber  reinforced  cementitious 
composites  is  very  complex  because  of  the  presence  and  the  combined 
action  of  several  bond  components.  These  include:  1)  physical  and 
chemical  adhesion  between  fiber  and  matrix,  2)  the  mechanical  component 
of  bond  such  as  in  deformed,  crimped  and  hooked  fibers,  3)  fiber  to  fiber 
interlock,  or  entanglement,  which  exists  in  SIFCON  even  before  addition 
of  the  matrix,  and  4)  friction  which  is  greatly  influenced  by  confinement. 
The  addition  of  latex  to  a  cement  matrix  may  add  significantly  to  the 
magnitude  of  the  "adhesion"  component  of  bond  especially  in  the 
precracking  state.  The  shear  stress-slip  relation  of  a  smooth  fiber  may  be 
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substantially  different  from  that  of  a  hooked  or  deformed  fiber.  Both  may 
lead  to  the  same  maximum  shear  yet  have  drastically  different  slip 
characteristics  and  pull-out  energies.  All  the  above  factors  must  be 
accounted  for  to  better  undestand  the  behavior  of  fiber  reinforced  cement 
composites  and  develop  models  to  predict  their  mechanical  and  fracture 
properties. 

1.2.  OBJECTIVES 

The  main  objective  of  this  research  is  to  investigate  the  various 
fundamental  mechanisms  of  bond,  generate  a  related  comprehensive 
experimental  data  base,  and  develop  rational  analytical  models  to  describe 
bond  in  fiber  reinforced  cement  based  composites.  A  particular  focus  is 
placed  at  characterizing  bond  by  a  bond  shear  stress  versus  slip 
relatiohship.  Such  a  relationship  is  considered  to  be  a  constitutive 
property  of  the  interface,  and,  for  given  conditions,  is  considered  to  be 
location  independent.  It  is  also  the  objective  of  this  research  to  better 
assess  the  shear  stress  distribution  along  the  fiber  as  well  as  the 
distribution  of  normal  stresses  in  the  fiber  and  the  matrix. 

1.3  SUMMARY 

Following  is  a  brief  summary  of  what  was  achieved  in  this  study; 

1.  A  review  of  the  state-of-the-art  on  bond  in  conventional  reinforced 
and  prestressed  concrete  was  carried  out.  Particular  attention  was 
given  to  the  study  and  determination  of  a  bond  shear  stress  versus 
slip  relationship  (Chapter  II)  which  is  believed  to  be  a  characteristic 
property  of  the  interface  to  the  same  extent  as  a  stress-strain  curve  in 
compression  or  tension.  It  is  observed  that  the  derivative  of  the  slip 
between  the  reinforcemei  .  and  the  surrounding  matrix  is  the 
difference  between  the  strain  in  the  fiber  and  that  in  the  matrix  at  a 
given  point  along  the  reinforcement.  Experimental  approaches  to  the 
determination  of  a  shear  stress-slip  relationship  are  described,  and 
existing  models  for  bond  shear  stress  distribution  are  evaluated, 
discussed  and  compared  (Chapter  III).  The  most  important  factors 
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affecting  bond  are  reviewed.  It  is  observed  that,  while  the  concept 
of  a  bond  shear  stress  versus  slip  relationship  has  been  applied  to 
conventional  reinforcement,  such  a  relationship  has  so  far  not  been 
developed  for  fiber  reinforced  concrete. 

2.  A  state-of-the-art  review  on  bond  in  fiber  reinforced  cement 
composites  is  presented  in  Chapter  IV.  The  review  addresses  both 
experimental  work  and  analytical  studies.  Factors  affecting  bond 
and  bond  mechanisms  in  cement  based  composites  are  discussed. 
Observations  which  are  either  agreed  upon  or  a  subject  of 
disagreement  between  researchers  are  pointed  out. 

3.  A  basic  analytical  study  of  a  typical  pull-out  test  is  developed  in 
Chapter  V.  The  test  consists  in  a  fiber  embedded  in  a  matrix  and 
being  pulled  out  at  its  free  end.  The  study  is  subdivided  into  two 
parts,  one  dealing  with  what  is  defined  as  the  primal  problem  and  the 
other  dealing  with  the  dual  or  reciprocal  problem.  The  primal 
problem  focuses  on  predicting  the  pull-out  load  versus  end  slip 
response  of  a  pull-out  test  based  on  the  knowledge  of  an  idealized 
bond  shear  stress  versus  slip  relationship.  For  small  values  of  slip, 
the  idealized  relatiohship  is  assumed  to  be  linear  perfectly  frictional, 
and  characterized  by  three  parameters,  namely  a  shear  modulus,  a 
maximum  shear  strength  which  includes  adhesion  and  friction,  and  a 
frictional  shear  stress  which  may  be  different  from  the  shear 
strength.  The  dual  problem,  which  is  the  reciprocal  of  the  primal 
problem,  predicts  a  bond  shear  stress  versus  slip  relationship  based 
on  an  experimentally  obtained  pull-out  curve.  For  large  slips,  in 
both  problems,  decaying  friction  was  considered.  It  was  modeled 
using  the  concept  of  misfit  which  accounts  for  the  influence  of 
normal  pressure  due  to  shrinkage  of  the  matrix,  and  including  the 
effect  of  fiber  Poisson’s  ratio.  The  solution  of  the  dual  problem 
forms  the  basis  for  the  interpretation  of  the  results  of  the 
experimental  program  undertaken  in  this  investigation. 

4.  A  comprehensive  experimental  program  on  pull-out  and  pull-through 
tests  is  described  in  Chapter  VI.  .  The  pull-out  load  versus  end  slip 
response  was  carefully  measured.  Parameters  included  three 
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different  types  of  fibers  (smooth,  deformed,  hooked),  three  different 
mortar  matrices  with  low,  medium,  and  high  strengths,  one  slurry 
paste  matrix,  and  additives  such  as  latex,  fly  ash,  and  microsilica. 
Also  the  medium  from  which  the  fiber  was  pulled  out  included  a 
control  mix  without  fibers,  mixes  with  1%,  2%  and  3%  fibers  by 
volume,  and  a  SIFCON  matrix  containing  about  11%  fibers  by 
volume.  For  smooth  fibers,  five  different  diameters  and  three 
different  embeddment  lengths  were  investigated.  Chapter  VI  also 
contains  comparisons  between  experimental  pull-out  load  versus  slip 
curves  and  curves  predicted  by  the  analytical  model  developed  in 
Chapter  V.  Typical  bond  shear  stress  versus  slip  curves  are  derived 
using  the  analytical  approach  described  in  Chapter  V  and  are  input  in 
the  solution  equations  of  the  primal  problem  to  predict  typical  pull¬ 
out  load  versus  slip  curves  and  compare  them  with  the  experimentally 
observed  curves. 

5.  Two  different  and  independent  analytical  models  were  developed  in 
Chapters  VII  and  VIII.  These  models  are  also  independent  of  the 
pull-out  problem  described  in  Chapter  V  and  VI.  Both  models 
assume  that  the  composite  has  a  square  packing  distribution  of  fibers 
and  that  a  representative  unit  made  out  of  a  prism  of  matrix 
containing  one  fiber  can  be  analysed  in  lieu  of  the  composite  as  a 
whole.  The  first  model  is  based  on  an  assumed  bond  stress  versus 
slip  relationship  that  consists  of  a  linear  ascending  portion  followed 
by  a  constant  purely  frictional  portion.  It  leads  to  prediction 
equations  for  the  stress  distribution  in  the  fiber,  the  matrix  and  the 
interface.  Two  numerical  examples  illustrating  the  application  of  the 
procedure  are  described  in  Chapter  VII.  The  second  model 
developed  in  Chapter  VIII  is  based  on  a  mechanism  of  force  transfer 
between  the  fiber  and  the  matrix.  The  mechanism  is  similar  to  the 
truss  analogy  method,  except  that  the  matrix  is  assumed  uncracked 
and  the  external  tensile  force  is  applied  to  the  fiber.  Thus  the  model 
applies  particularly  well  to  the  case  of  continuous  fibers  or  long 
fibers  bridging  more  than  one  crack.  The  model  leads  to  prediciton 
equations  for  the  interfaciai  shear  stresses  as  well  as  the  normal 
stresses  in  the  fiber  and  in  the  matrix.  It  also  led  to  the  analytical 
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prediction  of  a  bond  modulus  as  well  as  the  prediction  of  the 
precracking  strength  of  the  composite.  Chapter  VIII  also  contains  a 
parametric  evaluation  of  various  composite  properties  using  the 
model. 

14  CONCLUDING  REMARKS 

Specific  and  detailed  conclusions  related  to  the  experimental 
investigation  as  well  as  to  the  analytical  models  developed  in  this  study  can 
be  found  at  the  end  of  Chapters  V  to  VIII.  The  following  conclusions 
address  more  generally  the  overall  objective  of  this  research. 

1.  For  conventional  reinforcing  bars  and  prestressing  strands,  the  bond 

shear  stress  versus  slip  relationship  can  be  determined 

experimentally  by  instrumenting  the  bar  along  its  length  using  for 
instance,  strain  gages.  Such  an  approach  cannot  be  used  in  the  case 
of  small  fibers.  In  this  study,  bond  shear  stress  versus  slip 
relationship  were  derived  indirectly  from  a  pull-out  test  in  which  the 
load  and  the  slip  were  measured  simultaneously. 

2.  A  methodology  was  developed  to  derive  a  bond  shear  stress  versus 
slip  relationship  for  various  types  of  smooth  fibers.  A  typical 
relationship  would  be  elastic  linear  up  to  a  maximum  shear  stress, 
■fnn**  then  drops  suddently  to  a  maximum  frictional  stress,  Xf,  then 
decays  exponentially  to  zero  with  increasing  slip.  It  is  observed  that 
for  hooked  fibers,  the  relationship  can  be  obtained  from 
superimposing  two  relationships,  one  characteristic  of  a  smooth  fiber 
without  the  hook  and  the  other  characteristic  of  the  hook  alone.  In  a 
hooked  fiber,  the  hook  tends  to  straighten  under  pull-out  leading  to  a 
characteristic  load  versus  slip  response.  However,  the  slip  at  which 
the  load  is  maximum  is  substantially  larger  than  that  for  which  the 
load  is  maximum  for  a  smooth  fiber.  Deformed  fibers  have 
mechanical  deformations  uniformely  distributed  along  their  length. 
The  fibers  used  in  this  study  were  dented  along  segments  about  0.2 
in  long  each.  Thus  they  had  a  section  that  varied  every  0.2  in.  They 
too  seem  to  respond  through  the  surperpositon  of  two  mechanisms, 
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one  due  to  the  adhesion  similar  to  that  of  a  smoo'h  fiber,  and  the 
other  due  to  the  mechanical  bearing  against  the  surrounding  matrix 
induced  by  the  variation  in  cross  section. 

For  small  values  of  slip,  such  as  encountered  up  to  the  peak  null-out 
load  in  a  typical  pull-out  test,  the  bond  shear  stress  versus  slip 
relationship  can  be  assumed  to  have,  as  a  first  approximation,  a 
constant  frictional  property.  That  is,  the  decay  in  friction  with  the 
slip  is  ignored.  Such  an  assumption  greatly  simplifies  the  analytical 
model  and  leads  to  a  negligeable  error. 

The  medium  (slurry,  plain  matrix,  fiber  reinforced  matrix)  from 
which  a  fiber  pulls-out  influences  the  pull-out  load  versus  slip  curve. 
It  was  observed  that  the  presence  of  fibers  tends  to  improve 
primarily  the  post  peak  response  in  fiber  reinforced  mortar,  and  the 
peak  load  as  well  as  the  post  peak  behavior  in  SIFCON.  The 
addition  of  latex  to  the  matrix  leads  to  a  significant  improvement  in 
the  peak  load  (up  to  four  times),  but  has  no  effect  on  the  post-peak 
response.  The  addition  of  microsilica  does  not  seem  to  improve  the 
bond  characteristics  of  smooth  fibers  beyond  what  is  predicted  from 
the  induced  improvement  in  matrix  strength.  The  addition  of  fly  ash 
to  the  matrix  led  to  very  small  improvement  in  the  overall  response. 

Everything  else  being  equal,  the  slip  at  maximum  load  for  hooked  or 
deformed  fibers  is  one  to  two  orders  of  magnitude  larger  than  that  of 
a  smooth  fiber.  This  is  due  to  the  mechanical  contribution  of  the 
hook  in  hooked  fibers  and  the  non-uniform  section  in  deformed 
fibers.  As  a  result,  the  pull-out  work  up  to  the  peak  load  for  hooked 
and  deformed  fibers  can  be  ten  to  twenty  times  that  of  a  smooth 
fiber. 

Everything  else  being  equal,  the  pull-out  work  of  a  hooked  or  a 
deformed  fiber  up  to  complete  pull-out  can  be  three  to  four  times  that 
of  a  smooth  fiber. 
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7.  A  deformed  fiber  shows  in  the  post  peak  load  response  a  cyclic 
behavior  that  is  characteristic  of  the  way  the  non  uniform  section 
moves  in  its  surrounding  tunnel  of  matrix. 

8.  Observed  values  of  maximum  pull-out  loads  on  single  fibers  lead  to 
the  following  ranges  for  the  maximum  bond  stress  calculated  as  the 
maximum  pull-out  load  divided  by  the  embedded  surface  area:  from 
150  to  400  psi  for  smooth  fibers,  from  500  to  1000  psi  for  hooked 
fibers,  and  from  400  to  950  psi  for  deformed  fibers.  It  should  be 
observed  that,  although  differences  between  smooth  fibers  and 
deformed  and  hooked  fibers  are  significant,  they  tend  to  be  smaller 
when  the  whole  composite  is  considered.  This  is  because  the  pull¬ 
out  load  per  fiber  decreases  when  the  number  of  fibers  pulling  out 
from  the  same  area  increases.  The  addition  of  latex  to  the  matrix  led 
to  bond  stresses  as  high  as  about  1400  psi.  However,  the  frictional 
stress  was  similar  to  that  obtained  without  latex. 

9.  Observed  frictional  shear  stress  values  for  smooth  fibers  ranged  from 
about  170  to  290  psi  for  the  range  of  variables  tested  in  this  study. 

10.  For  smooth  fibers,  bond  shear  modulus  (obtained  from  the  slope  of 
the  initial  portion  of  the  shear  stress  versus  slip  curve)  varied  from 
about  1x10^  to  9x10^  lb/in3. 

1.5  RECOMMENDATIONS  FOR  FUTURE  WORK 

Following  is  a  list  of  recommendations  for  future  research  related  to 
bond  in  fiber  reinforced  cement  composites.  The  list  is  by  no  means 
exhaustive  but  reflects  several  areas  that  were  touched  upon  in  this 
investigation. 

1.  There  is  need  to  better  understand  frictional  decay  with  the  extent  of 
slip  and  to  model  the  mechanisms  involved.  A  model  based  on 
damage  mechanics  may  be  most  approriate.  This  should  help  in 
better  predicting  the  descending  branch  of  the  load  versus  slip 
response  in  a  typical  pull-out  test.  The  analytical  work  may  be 
coupled  with  an  experimental  investigation  using,  for  instance,  a 


comprehensive  microscopic  analysis  of  the  history  of  debonding  with 
slip. 

There  is  a  need  to  better  understand  the  mechanical  component  of 
bond  for  a  deformed  or  a  hooked  fiber  and  to  model  such  behavior. 
It  is  observed  that  the  maximum  contribution  of  the  mechanical 
component  of  bond  occurs  much  after  the  maximum  adhesion-friction 
is  developed.  This  contribution  seems  to  significantly  improve  the 
pull-out  work  up  to  the  peak  load. 

The  phenomenon  of  fiber  entanglement  or  interlock  found  in  highly 
reinforced  fiber  cement  composites  such  as  SIFCON,  should  be  more 
thoroughly  addressed  both  experimentally  and  analytically. 

The  models  developed  in  Chapters  VII  and  VIII  in  which  a  tensile 
prism  containing  one  fiber  is  analysed,  assumed  that  the  load  is 
applied  to  the  fiber  and  transferred  to  the  matrix.  The  reciprocal 
problem  should  be  addressed,  that  is,  to  apply  the  load  to  the  matrix 
and  study  the  stress  transfer  to  the  fiber.  The  two  solutions  should 
help  in  modeling  the  behavior  of  the  composite  before  and  after 
cracking.  Such  a  study  can  also  be  coupled  with  a  finite  element 
solution  to  check  the  validity  of  the  model. 

The  effect  of  strain  rate  on  the  shape  of  the  bond  shear  stress  versus 
slip  curve  should  be  investigated.  Such  a  study  would  provide  a 
better  understanding  of  the  effects  of  strain  rate  on  the  mechanical 
properties  of  the  composite. 

There  is  need  to  investigate  the  modulus  of  elasticity  of  the 
composite  and  to  develop  a  rational  prediction  model  whereas  the 
bond  shear  stress  versus  slip  relationship,  which  characterizes  the 
interface,  is  accounted  for  to  the  same  extent  as  the  properties  of  the 
matrix  and  the  fiber. 

The  theory  developed  in  Chapter  V  regarding  the  pull-out  process 
should  be  used  in  a  fracture  mechanics  framework  to  predict  the 
post-cracking  response  of  the  composite.  The  shear  stress  versus 
slip  relationship  can  be  suitably  integrated  with  the  constant  crack 


opening  angle  assumption  to  predict  the  onset  of  fracture  and  the 
fracture  energy. 
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CHAPTER  II 


BOND  SHEAR  STRESS  -  SLIP  RELATIONSHIP 


2.1  GENERAL 

The  general  relationship  between  local  bond  stress  and  slip  for  dif¬ 
ferent  types  of  tension  reinforcement  is  illustrated  in  Fig.  2. 1.  All  types  of 
reinforcement  display  a  similar  behavior  up  to  the  point  where  the  adhesion 
breaks  (ascending  branch).  Plain  bars  will  slide  out;  the  failure  is  hence  of 
a  sliding  nature.  As  for  ribbed  bars,  further  loading  mobilizes  the 
mechanical  interlocking  between  the  lugs  and  the  concrete.  The  effects  of 
lug  spacing  and  depth  on  the  stress-slip  relationship  will  be  discussed 
farther  ahead.  This,  in  effect,  creates  circumferential  tensile  cracks  in  the 
surrounding  concrete  (see  Fig.  2.2),  and  thus  reduces  the  stiffness  of  the 
bond  stress-slip  relationship  (i.e.  the  slope  of  the  ascending  portion  of  the 
curve  becomes  smaller).  The  existence  of  these  internal  cracks  has  been 
proven  experimentlally  by  Goto  [18J.  For  unconfined  concrete  or  in  the 
case  of  a  small  volume  of  surrounding  concrete,  the  internal  cracks 
propagate  to  the  surface  resulting  in  bond  failure.  On  the  other  hand,  if  the 
concrete  is  well  confined,  the  load  can  be  increased  beyond  the  creation  of 
these  cracks,  up  to  a  point  where  the  maximum  bond  stress  (i.e.  the  bond 
strength)  is  reached.  At  this  stage,  shear  cracks  develop  as  shown  in  Fig. 
2.3.  After  the  maximum  bond  stress  is  reached,  further  slip  causes  a 
reduction  in  the  bond  resistance  until  bond  stresses  all  become  of  a  fric¬ 
tional  nature.  The  friction  that  takes  place  is  actually  between  the  rough 
concrete  at  the  cylindrical  surface  where  the  shear  cracks  developed  [14]. 

It  is  worth  noting  that  there  is  a  general  tendency  to  believe  that  the 
shear  stress-slip  relationship  is  of  a  local  nature,  meaning  that  this  rela- 
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Longitudinal  stction  of  axially  loaded  specimc 


Cross  section 


tiqNening  force  on  bar 
(due  to  wedge  action  and 
deformation  of  teeth  of 
comb -tike  concrete) 


Fig.  2.2 


Internal  bond  cracks  and  forces  acting  on  concrete  (From 
Ref.  18). 
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Fig,  2.3  '  Shear  cracks  in  concrete  keys  between  lugs  (From  Ref.  14). 
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Fig.  2.4  -  Local  bond  stress-slip  relationships  for  reinforcing  steel  bars 
after  different  researchers  (From  Ref.  14). 
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tionship  is  not  unique,  but  depends  instead  on  the  location 
[1,16,23,25,45,55],  However,  Edwards  and  Yannopoulos  [12]  stated  that 
"the  maximum  bond  stress  value  would  not  vary  with  distance  from  a 
loaded  end,  or  a  crack  face,  in  a  concrete  member."  Mirza  and  Houde  [37] 
reported  similar  findings. 

2.2  MODELS  AND  STUDIES  OF  STRESS-SLIP  RELATIONSHIPS 

Extensive  experimental  studies  were  undertaken  in  recent  years  to 
determine  the  characteristic  behavior  of  local  bond  stress-slip  relationship 
of  deformed  reinforcing  steel  bars  [11,14,37,45,46],  and  several  mathe¬ 
matical  models  were  proposed.  However,  depending  on  the  case,  the 
characteristic  behavior  of  the  bond  stress-slip  relationship  shows  signifi¬ 
cant  scatter.  This  scatter  can  be  illustrated  in  Fig.  2.4. 

Nilson  [45]  carried  out  tests  on  6x6x18  in.  prism  specimens  to  study 
the  local  bond  stress-slip  relationship  on  concrete-steel  interfaces.  #8  bars 
were  used  in  these  tests.  All  specimens  were  subjected  to  concentric 
tension  forces  acting  on  both  ends  of  the  embedded  bar.  As  was  mentioned 
earlier,  it  was  concluded  from  this  study  that  the  bond  stress-slip 
relationship  varies  along  the  embedment  length,  and  is  thus  not  unique. 
The  following  equations  were  suggested  by  Nilson  to  relate  the  bond  stress 
to  the  slip  and  the  distance  from  the  loaded  end: 

=  3100  (1.43  c  +  1.5)  Sx  Vrl  (2.1) 

^  (1.43  c  -(■  1.5)  Vf^  (2.2) 

where  is  the  local  bond  stress  in  psi,  c  and  represent  the  distance 
from  the  loaded  end  (in.)  and  the  slip  respectively,  and  f'c  is  the  concrete 
compressive  strength  in  psi.  The  slip  is  the  relative  displacement  between 
the  concrete  and  the  steel  bar. 

Comparisons  between  Eq.  2.1  and  the  experimental  resu'ts  for  dif¬ 
ferent  values  of  c  are  shown  in  Fig.  2.5. 
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Mirza  and  Houde  [37]  tested  62  pull-out  tension  specimens  rein¬ 
forced  with  steel  bars  ranging  from  #4  to  #8.  Different  concrete  strengths 
and  different  size  prisms  were  studied.  The  authors  found  out  that  the  bond 
stress  is  independent  of  the  distance  from  the  loaded  end.  This  is  in 
contradiction  with  what  was  stated  by  Nilson.  Mirza  and  Houde  suggested 
the  following  expression  for  the  bond  stress-slip  relationship: 

=  1.95x106  8^-2.35x109  5^2  +  1  39x1012  5^3.0.33x1015 

(2.3) 

where  is  the  local  bond  stress  in  psi,  normalized  to  ^|T~^  (f^-  =  5000  psi) 
and  is  the  local  slip  in  inches.  Concrete  strengths  other  than  5000  psi 

can  be  accounted  for  by  multiplying  the  right  hand  side  of  Eq.  2.3  by 
\  5000- 

Recently,  an  extensive  study  was  conducted  by  Eligehausen,  Popov 
and  Bertero  [14]  to  determine  the  local  bond  stress-slip  relationship  under 
generalized  excitations.  Pull-out  tests  were  conducted  on  concrete  speci¬ 
mens  with  small  embedment  lengths  (5  times  the  bar  diameter).  Different 
bar  diameters  dj,  (#6  to  #10)  and  different  concrete  strengths  (4000  to  8000 
psi)  were  used.  In  addition,  the  effect  of  bar  spacing  and  the  amount  of 
confining  reinforcement,  to  simulate  the  conditions  prevailing  in  an  actual 
joint,  were  also  studied.  The  authors  observed  that  confinement  helps 
restrain  the  propagation  of  the  secondary  internal  tensile  cracks  around  the 
bar  lugs,  and  thus  increase  the  bond  strength  significantly.  Specimens  with 
no  confining  reinforcement  developed  a  peak  bond  stress  of  about  6  MPa  as 
opposed  to  the  13.5  MPa  developed  in  confined  specimens.  It  was  also 
observed  that  the  bar  diameter  has  little  influence  on  the  maximum  bond 
strength.  Moreover,  it  was  observed  that  for  a  given  slip,  the  bond 
strength  increases  linearly  with  The  following  expression  was 

suggested  for  the  ascending  portion  of  the  local  bond  stress-slip  curve: 


t 

^max  1 J 


(2.4) 
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where  : 

=  local  bond  stress 

'^max  “  maximum  bond  strength 

Sjj  =  local  slip 

Si  =  local  slip  corresponding  to 

The  average  observed  values  [14]  of  ^1  were  respectively 

13.5  MPa  and  1.0  mm  for  a  30  MPa  concrete  strength.  Typical  plots  of  the 
bond  stress-slip  relationship  for  the  monotonic  loading  are  shown  in  Fig, 
2.6. 

Stocker  and  Sozen  [58]  carried  out  pull-out  tests  on  1/2"  prestressing 
strands  with  short  embedment  length  (0.5-2  in.)  and  different  concrete 
strengths.  The  average  value  of  maximum  bond  strength  in  their  study  was 
about  0.4  ksi  for  a  2500  psi  concrete  strength,  0.56  ksi  for  6000  psi 
concrete  and  0.58  ksi  for  a  7500  psi  concrete. 

Edward  and  Picard  [11]  tested  a  total  of  12  specimens  to  study  the 
bond  characteristics  of  1/2"  prestressing  strands.  The  embedment  length  of 
the  strands  was  1.5"  and  the  average  concrete  cube  strength  of  all  the 
specimens  was  9640  psi.  It  was  observed  that  the  average  value  of  the 
maximum  bond  strength  decreases  with  increasing  concrete  cover.  This 
observation  was  attributed  to  the  settlement  of  the  concrete  under  the  strand 
which  was  held  in  a  horizontal  position  during  casting.  They  concluded 
that  the  bond-slip  relationship  for  strands  is  almost  linear  up  to  the  peak 
bond  stress  and  remains  constant  and  equal  to  its  maximum  value  there¬ 
after,  provided  no  longitudinal  cracks  appear  in  the  specimens.  The 
maximum  bond  stress  in  this  test  varied  between  0.33  and  0.67  ksi,  and  the 
stiffness  of  the  bond  stress-slip  curve  varied  between  665  and  879  k/in^. 
Typical  plots  of  bond  stress-slip  relationships  in  their  test  is  shown  in  Fig. 
2.7. 


Bond  stress-slip  relationships  for  reinforcing  steel 
(From  Ref.  14). 
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Fig.  2.7  -  Local  bond  stress-slip  relationship  for  prestressing  strands 
(From  Ref.  1 1). 
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2.3  EXPERIMENTAL  APPROACH:  STRESS-SLIP  THROUGH  STRAIN 
MEASUREMENT 

Another  approach  to  the  problem  is  the  experimental  way,  which  is 
easy  conceptually,  but  rather  difficult  to  perform,  especially  in  the  case  of 
fiber  reinforcement.  It  was  used  by  Nilson  [45]  to  obtain  stress-slip 
curves. 

The  experimental  determination  of  the  force-displacement  relation¬ 
ship  in  bond  has  two  components:  the  bond  stress  must  be  found,  and  the 
slip  between  steel  and  concrete  must  be  measured.  Since  both  these  quan¬ 
tities  vary  with  distance  along  the  bar  at  any  load,  it  is  desirable  to  estab¬ 
lish  each  as  a  continuous  function  of  distance. 

The  average  unit  bond  stress,  expressed  in  units  of  force  per  unit 
interfacial  area,  is  an  appropriate  parameter.  This  stress  is  given  by: 


dT 

V  dx 

Y  dx 


(2.5) 

(2.6) 


Eq.  (2.6)  can  be  easily  restated  in  terms  of  the  tensile  strain  e,  in  the 


steel: 


u 


where: 


X 

dT 

¥ 


As  Eg  dEg 
\{r  dx 


(2.7) 


average  bond  stress  at  location  considered 
change  in  bar  tension  in  length  dx 
bar  perimeter 


bar  cross-sectional  area 


~ 

dfs  ”  change  in  steel  tensile  stress  in  length  dx 

Ej  =  modulus  of  elasticity  of  steel 

Bond  stress  at  any  location  along  the  interface,  at  any  load,  is  thus 
proportional  to  the  slope  of  the  steel  strain  distribution  curve  at  that  point 
and  that  load.  Eq.  2.7  assumes  that  steel  is  still  in  the  elastic  range.  In 
reinforced  concrete,  this  fact  is  usually  true  because  reinforced  concrete 
members  are  usually  designed  in  such  a  way  that  under  working  loads,  the 
stress  level  in  steel  does  not  exceed  50%  of  yield.  This  assumption  is  also 
sound  in  the  case  of  steel  fiber  reinforced  concrete  since  it  has  been 
established  that  steel  fibers  usually  pull-out  before  they  yield. 

To  get  the  shear  stress-slip  relationship  using  this  method,  all  is 
needed  is  the  strain  distribution  for  both  steel  and  concrete  along  the  steel, 
be  it  a  reinforcing  bar  or  a  fiber.  Indeed,  bond  slip  can  be  found  indirectly 
as  the  difference  between  the  steel  displacement  and  that  of  concrete,  at  any 
location  on  the  interface,  with  respect  to  some  transverse  reference  plane. 
In  case  of  reinforced  concrete,  this  plane  should  be  the  median  of  the 
portion  of  the  bar  between  two  primary  cracks.  In  the  case  of  fiber 
reinforced  concrete,  if  the  fiber  is  not  crossed  by  any  crack,  the  reference 
plane  would  simply  be  the  median  of  the  fiber,  otherwise  this  plane  would 
be  the  plane  of  the  crack.  The  displacement  function  for  the  concrete  and 
the  steel  may  be  obtained  by  integration  of  the  strain  functions. 

Nilson  [45]  obtained  the  steel  strain  distribution  as  follows.  The 
reinforcing  bar  was  sawed  longitudinally  on  a  diametral  plane,  and  slots 
milled  along  the  centerline  of  each  cut  surface.  Electric  resistance  gages 
were  then  attached  at  close  intervals  along  the  bottom  of  each  slot,  after 
which  the  two  halves  of  the  bar  were  rejoined  using  epoxy  cement.  As  for 
the  concrete,  he  used  electric  resistance  strain  gages  cemented  between  two 
thin  polyester  resin  blocks.  The  outer  surface  of  the  resin  is  coated  with 
coarse  grit  material  to  provide  good  bonding  with  the  concrete.  Such  gages 
were  embedded  in  test  specimens  at  intervals  along  the  interface,  and  1/2 
in.  from  it. 


In  Fig.  2.8,  the  variation  of  steel  strain  is  plotted  as  a  function  of 
distance,  and  so  is  the  concrete  strain.  Integration  of  each  strain  function 
along  the  bar  length  provides  the  displacement  of  the  steel  and  the  concrete, 
measured  from  the  reference  plane. 

With  the  bond  stress  calculated  from  the  slope  of  the  steel  strain 
curve  using  Eq.  2.7,  and  the  bond  slip  found  by  numerical  integration,  one 
point  on  the  bond  stress-slip  curve  is  established.  By  increasing  the  load 
incrementally,  subsequent  points  may  be  obtained  for  larger  bond  stress 
and  slip,  hence  defining  the  complete  curve  at  a  given  location. 


CHAPTER  III 


EXISTING  MODELS  FOR  BOND  SHEAR  STRESS  DISTRIBUTION 


3.1  GENERALITIES 

In  composite  materials,  the  study  of  bond  is  an  unavoidable  step  to¬ 
wards  the  understanding  of  the  way  forces  are  transmitted  from  the  matrix 
to  the  reinforcement.  The  idea  of  reinforcing  cementitious  matrices  stems 
from  this  concept.  Since  these  matrices  fail  at  much  lower  tensile  strains 
than  steel  does,  once  the  matrix  (concrete  or  mortar)  cracks,  the  tension 
that  cannot  be  taken  by  the  matrix  is  transmitted  to  the  reinforcement 
(reinforcing  bars,  prestressed  strands  or  discontinuous  fibers).  If  there 
was  no  bond  between  the  matrix  and  the  reinforcement,  i.e.  a  condition  of 
free  fnctionless  •'lippagc,  then  the  composite  would  fail  as  soon  as  the 
matrix  cracks,  and  the  reinforcement  would  not  have  taken  any  share  of  the 
tension.  Bond  is  needed  for  tension,  as  well  as  for  compression,  although 
the  force  transmission  and  the  mode  of  failure  are  different  for  the  two 
cases. 


In  the  case  of  fiber-reinforced  concrete,  bond  is  important  in  many 
respects.  .There  is  hardly  any  formula  giving  the  tensile  strength  of  a  fiber- 
reinforced  composite  that  does  not  contain  a  bond  term.  The  bond  is  also 
important  for  the  understanding  of  strain  rate  effects  on  the  tensile 
properties  of  fiber  reinforced  concrete.  This  bond  is  usually  the  average 
bond  stress  along  the  interface  which  cannot  be  determined  accurately 
unless  the  bond  stress  distribution  along  the  interface  is  known.  Various 
approaches  have  been  used  by  researchers  to  determine  bond  stress 
distributions.  In  this  chapter,  three  different  analytical  methods  used  in  the 
determination  of  the  bond  stress  distribution  along  reinforcing  bars  in  rein- 
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forced  concrete  are  reviewed.  §3.5  briefly  discusses  a  new  approach  to  the 
bond  action  between  concrete  and  a  typical  deformed  steel  bar.  All 
methods  presented  are  briefly  compared  to  each  other  in  §3.6  with  regards 
to  accuracy,  practicality,  and  simplicity. 

Usually,  it  is  the  average  shear  stress  along  the  interface  rather  than 
the  stress  distribution  that  is  needed.  However,  when  the  shear  stress 
distribution  along  the  interface  is  known,  a  simple  integration 
(mathematical  or  numerical)  gives  the  value  of  the  average  bond.  The 
following  three  sections  outline  three  different  methods  for  finding  the 
stress  distribution  along  an  interface. 

3.2  ANALYTICAL  METHOD  BASED  ON  ASSUMED  BOND  SHEAR 
STRESS-SLIP  RELATIONSHIP 

The  idea  behind  this  method  is  simple  conceptually,  and  the  mathe¬ 
matics  involved  are  not  particularly  cumbersome.  This  method  requires  an 
analytical  relationship  between  the  bond  shear  stress  and  the  slip.  It  also 
assumes  that  the  bond  shear  stress-slip  relationship  is  unique.  While  this 
assumption  may  not  be  exact,  as  was  pointed  out  earlier  in  Chap.  2,  the 
fact  remains,  however,  that  this  method  is  still  valid,  especially  when  the 
shear  stress  distribution  is  used  for  averaging  the  bond,  as  long  as  an 
average  stress-slip  relationship  is  used.  The  detailed  mathematical  aspects 
of  this  method  as  applied  to  fiber-reinforced  cementitious  composites  are 
developed  in  Chap.  5.  The  essence  of  this  method,  as  well  as  the  under¬ 
lying  rationale  are  briefly  described  next. 

It  is  first  noted  that  the  rate  of  change  of  the  slip  along  the  rein¬ 
forcement-matrix  interface  is  the  difference  between  the  strain  in  the  rein¬ 
forcement  and  the  strain  in  the  matrix.  This  same  relationship  can  be 
expressed  in  function  of  the  stresses  in  the  reinforcement  and  the  matrix, 
assuming  Hooke's  law  applies  for  both  of  them.  The  stress  in  the  matrix 
can  be  related  to  the  stress  in  the  reinforcement  in  a  linear  manner,  so  that 
the  rate  of  change  of  the  slip  with  distance  can  be  expressed  as  a  function 
of  the  stress  in  the  reinforcement.  The  second  derivative  of  the  slip  with 
respect  to  distance  is,  therefore,  proportional  to  the  rate  of  change  of  the 
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stress  in  the  reinforcement  with  distance,  which  is  directly  related  to  the 
shearing  stress  through  the  bond-slip  (x-S)  relationship  assumed  known. 
Using  the  (x-S)  relationship,  the  problem  is  reduced  to  a  second  degree 
differential  equation,  the  solution  of  which  provides  the  slip  distribution 
with  distance.  This  slip  distribution,  along  with  the  bond  shearing  stress- 
slip  relationship  gives  the  bond  shear  stress  distribution. 

This  method  does  not  involve  any  mathematical  approximations  so 
long  as  the  differential  equation  obtained  can  be  analytically  solved,  which 
is  not  necessarily  true.  Furthermore,  this  approach  yields  an  analytical  ex¬ 
pression  for  the  bond  stress,  rather  than  discrete  values  at  select  points 
along  the  interface.  The  shear  stresses  at  different  points  along  the 
interface  can  be  obtained  by  plugging  the  value  of  the  distance  from  the 
point  of  interest  to  the  reference  point  into  the  expression  of  the  shear 
stress. 


The  above  described  method  has  at  least  one  disadvantage,  in  that  it 
cannot  accommodate  bond  shear  stress-slip  relationships  that  cannot  be 
expressed  analytically.  Therefore,  (x-S)  relationships  expressed  graphi¬ 
cally  are  to  be  modeled  mathematically  first,  before  this  method  can  be 
applied.  Such  a  modeling  will  inevitably  introduce  an  error,  the  effect  of 
which  might  be  significant  depending  on  the  level  of  simplification 
adopted.  It  is  worth  mentioning  that  this  method  has  been  widely  used. 

3.3  NUMERICAL  METHOD  BASED  ON  EXPERIMENTALLY  DERIVED 
BOND  STRESS-SLIP  RELATIONSHIP 

Tassios  and  Yannopoulos  [62]  derived  a  computerized  method  ca¬ 
pable  of  determining  in  a  direct  and  analytical  way  the  stress  distributions 
in  the  reinforcement,  in  the  matrix,  as  well  as  the  bond  stress  distribution 
along  the  interface.  It  can  also  determine  all  the  corresponding  deformation 
characteristics.  Examples  of  applications  regarding  pull-out  tests  under 
monotonic  loading,  axially  loaded  elements  under  both  monotonic  and 
cyclic  loading,  and  flexural  members  under  reversible  loading  were  given. 

Fig.  3.1  shows  the  interrelationship  between  the  stress  in  the 
reinforcement,  the  stress  in  the  matrix,  the  shearing  stress  and  the  slip. 
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Local  Bond  Stress 
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Local  Slip 
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For  the  axially  loaded  element  shown  in  Fig.  3.2, 
ementary  equations  may  be  derived: 


AOc 


Ao, 


AS 


4  T 
V 
A. 


Ax 


Ac 


(es  -  ec)Ax 


the  following  el- 

(3.1) 

(3.2) 

(3.3) 


where  the  subscripts  s  and  c  refer  to  the  steel  reinforcement  and  the  con¬ 
crete  respectively,  c  is  the  normal  stress,  e  is  the  strain,  A  is  the  area,  S  is 
the  slip,  X  is  the  bond  stress  along  the  interface,  and  \\f  is  the  reinforcement 
perimeter. 

It  is  assumed  here  that  there  is  a  unique  relationship  between  local 
bond  t  and  local  slip  S  at  every  point  of  a  bar  inserted  in  a  mass  of 
concrete.  An  axially  tensioned  reinforced  concrete  element  after  cracking  is 
considered. 

A  typical  pullout  clement  (Fig.  3.3)  is  studied  analytically  by  di¬ 
viding  its  length  into  a  number  of  parts  of  sufficiently  small  length  Ax,  and 
applying  the  basic  relationships  of  Eqs.  3.1  -  3.3  at  each  consecutive 
section.  In  the  case  of  axially  loaded  elements,  only  half  of  their  length 
needs  to  be  considered  because  of  symmetry. 


The  following  algorithm  is  applied  to  compute  the  values  of  Oj,  Oq, 
S,  and  t  at  section  i-^l,  knowing  the  corresponding  values  at  section  i  and 
assuming  that  x  is  uniformly  distributed  over  the  elementary  length  Ax: 


^s,i+l  “ 


4  Xj 

Os.i - ^Ax 

V 

®c,i  ■  *  ®s.i  ) 


(3.4) 

(3.5) 
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Stresses  and  strains  on  a  finite  length  of  a  reinforced 
concrete  element  (Reproduced  from  Ref.  62). 
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With  Es^i+i  and  j  corresponding  to  Oj  ,+  (  and  Ogj  respectively  in 
the  Oj  -  Eg  curve,  and  Ej.  and  i  corresponding  to  (Jc,i+i  and  a^.  j  re¬ 
spectively  in  the  o^.  -  e^  curve: 

Ax  Ax 

^i+l  ~  Sj  -  (  Ej  +  £5  j  ) -^  +  (  E(.j+1  +  ^C.i  ^  ”2”  (3.6) 

where  corresponds  to  Sj+i  in  the  (x  -  S)  curve. 

The  boundary  conditions  for  the  pullout  element  are  shown  in  Fig. 
3.3.  Using  the  above  algorithms  and  boundary  conditions,  the  values  of 
^s.i’  ^c,i»  Sj,  and  Xj  can  be  computed  at  each  section  for  any  given  back 
slip  Sjj,  starting  from  the  back  face  (Sj  =  S^)  and  ending  at  the  front  one. 

As  can  be  seen  from  the  above  description,  the  numerical  method  is 
flexible,  easy,  and  short  if  a  computer  is  used.  Furthermore,  it  can  accom¬ 
modate  all  kinds  of  stress-slip  relationships,  no  matter  how  complicated 
they  are,  without  modeling  anything.  Finally,  the  error  introduced  can  ^>6 
reduced  by  reducing  the  value  of  Ax.  The  above  method  is  very  powerful 
and  useful,  and  can  accommodate  more  complex  situations,  such  as  cyclic 
loading. 

The  main  shortcoming  in  this  method  is  the  back  slip  that  cannot 
be  determined.  The  authors  did  not  suggest  any  way  to  determine  the  value 
of  Sb-  Another  possible  disadvantage  of  this  method  is  that  the  shear 
stress  along  the  interface  (or  axial  stress  in  the  reinforcement  or  matrix) 
cannot  be  found  at  any  point,  but  only  at  specific  points  or  "stations". 
This  means  that  the  bond  stress  curve  obtained  is  only  defined  by  discrete 
points,  joined  by  straight  lines,  or  perhaps  by  a  best  fit  curve.  Another 
source  of  error  is  that  the  integration  needed  to  obtain  the  average  bond 
stress,  has  to  be  done  numerically.  The  error  can  be  reduced  by  including 
more  stations  per  interface. 


3.4  ANALYTICAL  METHOD  BASED  ON  ASSUMED  BOND  STRESS 
DISTRIBUTION 
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surement  of  local  bond  stress  and  local  slip  along  a  stressed  reinforcing  bar 
of  a  tension  member  ia  difficult  and  very  sensitive  to  experimental  error. 
Moreover,  the  relationship  between  bond  stress  and  local  slip  has  been 
observed  to  vary  from  section  to  section  [45].  The  authors  of  this  work 
were  basically  trying  to  find  the  shear  stress  distribution  at  the  interface  of 
the  reinforcing  bar  and  the  surrounding  concrete,  as  well  as  the  normal 
stresses  in  the  bar  and  the  concrete  for  the  bar  shown  in  Fig.  3.4.  They 
introduced  a  differential  equation  relating  the  second  derivative  of  the  local 
slip  to  the  local  bond  stress,  along  with  five  boundary  conditions,  namely 
that: 


1. 

2. 

3. 

4. 

5. 


.  s 

Atx  =  2* 

Sx  =  0. 

At  X  =  0  , 

dSx_  , 
dx  ~  ■‘'s 

At  X  =  Lj  , 

II 

o 

At  X  =  0  , 

*  s 

At  X  =  2  ' 

o 

II 

where  is  the  local  slip,  s  is  the  length  of  the  specimen,  is  the  strain  in 
the  steel,  and  Lj  is  the  transfer  length,  defined  as  the  embedment  length 
required  to  satisfy  the  condition  that  at  x  =  L(,  the  strains  in  the  steel  and 
the  concrete  are  equal  to  each  other. 

The  authors  then  noted  that  if  a  stress-slip  relationship  is  to  be  as¬ 
sumed,  and  if  it  is  assumed  that  the  local  bond  stress  is  linearly  dependent 
on  the  local  slip,  then  the  resulting  solution  would  be  of  the  form: 

Sx  =  A  e''*  -*■  B  e-l^*  (3.7) 

where  A  and  B  are  constants  to  be  determined  from  the  boundary  condi¬ 
tions,  and  k  is  a  constant  dependent  on  the  material  properties  and  cross- 
sectional  dimensions. 


eoNO 

STRESS 
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The  problem  lies  in  the  fact  that  Eq.  3.7  has  only  two  unknown 
coefficients,  whereas  five  boundary  conditions  are  to  be  satisfied  simulta¬ 
neously.  Two  observations  were  made  by  Somayaji  and  Shah  namely  that  : 
1)  the  relationship  between  the  local  bond  stress  and  the  local  slip  is  not 
unique,  2)  it  is  experimentally  difficult  to  obtain  a  precise  relationship 
between  the  local  bond  stress  and  local  slip  that  can  be  mathematically 
modeled.  These  problems  led  them  to  the  following  alternative  approach. 


Instead  of  ass  U  nil  lig  a  relationship  between  the  bond  stress  and  the 
slip,  an  exponential  bond  stress  distribution  is  assumed  leading  to  the 
following: 

d^S, 


dx2 


A  e*  B  e*’^  +  C 


(3.8: 


Integrating  the  above  equation  twice  leads  to  the  following  expres¬ 
sion  for  the  local  slip  S^: 


Ae*  +  Be-*  +  Cy+Dx-t-E 


(3.9) 


where  the  five  coefficients  A,  B,  C,  D,  and  E  can  be  found  from  the  five 
boundary  conditions  mentioned  earlier.  After  determining  the  coefficients, 
one  can  evaluate  the  strains  in  steel  and  concrete  at  any  section  by  succes¬ 
sively  differentiating  Eq.  3.9,  and  using  the  relationships  between  the 
strains  in  the  reinforcement  and  the  matrix,  and  the  local  slip.  In  solving  a 
typical  pull-out  problem,  the  five  coefficients  were  found  to  depend  on  the 
transfer  length  Lp  as  well  as  the  crack  spacing  or  the  segment  length  s.  It 
was  also  noted  that  in  the  course  of  the  development  of  this  theory,  a  linear 
relationship  between  the  transfer  length  and  transfer  load  was  assumed, 
which  assumption  is  not  unreasonable. 

This  approach  is  unique,  and  Somayaji  and  Shah  are  the  only  in¬ 
vestigators  to  have  adopted  it.  The  authors  believe  that  the  method  can  be 
substantially  improved  if  Eq.  3.8  is  replaced  by: 


d2s. 

dx2 


A'  el^*  +  B'  e-‘‘* 


(3. 10) 
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Then,  by  integrating  twice,  the  following  expression  for  the  local 
slip  is  obtained: 

=  A  ekx  +  B  e->^^  +  Cx  +  D  (3.11) 

where  A  =  A'  /  and  B  =  B’  /  k^. 

k  in  these  expressions  is  a  function  of  the  cross-section,  and  should 
in  fact  reflect  the  stress-slip  relationship. 

This  method  overall  is  simple  and  easy  to  implement.  However,  it 
lacks  power  and  accuracy.  As  can  be  noted,  the  expressions  obtained  de¬ 
pend  only  on  the  boundary  conditions.  In  other  words,  the  shear  stress 
distribution  obtained  using  this  method  is  independent  of  the  stress-slip 
relationship,  meaning  that  the  stresses  in  the  steel  and  the  concrete,  as  well 
as  the  shear  distribution  along  the  interface  are  independent  of  the  bond 
properties  of  the  interface.  It  is  true  that  because  of  the  analytical  nature  of 
this  method,  the  shear  stress  can  be  found  "exactly"  at  every  point,  but 
again,  this  accuracy  is  in  contradiction  with  the  rough  assumption  made  as 
a  starting  point  that  the  slip  distribution  satisfies  Eq.  3.9.  Finally  it  should 
be  mentioned  that  Somayaji  and  Shah  were  able  to  derive  a  stress-slip  curve 
from  the  assumption  made  in  Eq.  3.9. 

3.5  YANKELEVSKY’S  THEORITICAL  APPROACH 

So  far,  all  the  methods  outlined  in  this  chapter  have  in  common  that 
the  bond  shear  stress  distribution  is  dictated  by  the  relationship  between 
the  stress  and  the  slip.  Of  course,  the  stress  distribution  and  the  stress- 
slip  curve  are  interrelated  in  that  one  can  always  be  derived  from  the  other. 
However,  assuming  one  in  order  to  obtain  the  other  is  not  the  only  valid 
approach.  Indeed,  Yankelevsky  [65]  adopted  a  completely  different  path  to 
study  the  bond  between  a  deformed  steel  reinforcing  bar  and  the 
surrounding  concrete. 

Yankelevsky' s  work  was  centered  on  the  problem  of  axially  loaded 
specimens,  where  a  single  reinforcing  deformed  bar  is  encased  concentri¬ 
cally  in  a  long  concrete  cylinder  with  the  bar  ends  exposed.  Goto  [18] 


performed  an  interesting  series  of  tests  in  which  such  specimens  were 
loaded  in  tension  and  found  that  shortly  after  primary  lateral  cracks 
formed,  internal  cracking  began  (Fig.  2.2).  These  numerous  cracks  form 
cones  with  their  apexes  near  bar  lugs.  The  angle  of  internal  cracking  is 
generally  within  the  range  of  45*  to  80*  to  the  bar  axis  and  many  crack 
angles  reach  60*.  These  cracks  were  found  to  start  at  a  steel  stress  level  of 
about  100  MPa  (14  ksi). 

The  model  in  question,  which  is  closely  related  to  Goto's  experi¬ 
mental  observations,  represents  the  interaction  between  uncracked  concrete 
and  a  deformed  bar  through  a  mechanical  system;  conceivably  it  could  be 
modified  to  apply  to  cracked  concrete  in  between  two  cracks. 

Fig.  3.5  shows  a  concrete  cylinder  having  a  length  L  and  a  diameter 
D,  in  which  a  deformed  steel  bar  is  centrally  located.  A  tensile  force  Fq  is 
applied  to  the  steel  bar  ends.  A  differential  element  dx  of  the  steel  bar  is 
located  to  the  left  of  the  axis  of  symmetry.  The  applied  tensile  force  pulls 
to  the  left  all  the  segments  located  to  the  left  of  the  axis  of  symmetry.  It  is 
assumed  that  the  surrounding  concrete  resists  the  pull  by  inclined 
compressive  forces  dC  (Fig.  3.6),  thus  reducing  the  axial  force  F  in  the  bar 
along  the  x  coordinate. 

The  structural  model  through  which  the  resisting  force  dC  is  trans¬ 
ferred  to  the  concrete  is  shown  in  Fig.  3.6.  Concrete  segments  in  com¬ 
pression  AB  and  A'B'  are  inclined  to  the  steel  bar  segment  at  angle  a  and 
react  on  the  rod  by  a  force  dC  per  unit  length  in  the  hoop  direction.  The 
circumference  BB'  is  the  locus  of  segment  centroids  of  a  typical  disc  EE' 
(Fig.  3.6).  The  centroid  B  (or  B')  in  Fig.  3.6  is  the  location  of  the  resul¬ 
tant  tensile  force  to  equally  distributed  tensile  stress  in  concrete  in  the  axial 
direction.  The  compressive  force  dC  per  unit  length  of  the  circumference 
BB'  in  the  inclined  segment  AB  (Fig.  3.6)  at  B  is  in  equilibrium  with  a 
radial  force  dN  per  unit  length  of  the  circumference  BB'  in  the  disc  EE'  and 
with  the  longitudinal  tensile  force  per  unit  length  dQ.  By  such  a 
mechanism,  the  tensile  force  F  in  the  steel  bar  is  reduced  at  a  rate  of  dF/dx 
and  is  increased  at  an  equal  rate  of  dT/dx  in  the  concrete.  T  is  the  total 
tensile  force  and  equal  to  27i  r  Q,  where  the  centroid  location  r  (Fig.  3.6)  is 
given  by: 


Schematic  representation  of  Yankelevsky' s 
(Reproduced  from  Ref.  65). 


specime 


Structural  model  of  Yankelevsky  (Reproduced  from  Ref.  65) 
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(3. 12) 


For  the  common  cases  where  d«D,  r  can  be  taken  as  D  /  3. 


From  this  point,  static  equilibrium  and  compatibility  are  used  to 
relate  the  forces  together.  A  differential  equation  in  F  (the  axial  force  in  the 
steel)  is  obtained  and  solved.  The  solution  obtained  was: 

F  =  C,  C2  e-P«  .  Fo  <313) 

where  n  is  the  modulus  ratio  (n  =  /  E^)  and  P  is  a  constant  dependent  on 

the  modulus  of  elasticity  of  steel,  on  that  of  concrete,  as  well  as  the 
dimensions  of  the  specimen. 

It  is  found  that  the  third  term  in  the  right  hand  side  of  Eq.  3.6  (say 
Fj)  is  the  portion  of  the  total  force  Fq  which  is  carried  by  the  steel  bar  ac¬ 
cording  to  its  relative  stiffness,  the  other  portion  (say  Fc  =  '  ^s)  being 

carried  by  the  concrete.  Strain  in  the  steel  bar  at  X=0  (Fig.  3.5)  has  its 
highest  value,  and  strain  in  the  concrete  there  is  zero.  As  distance  from  bar 
end  increases,  the  strain  in  the  concrete  increases  and  that  in  the  steel 
decreases.  At  a  distance  I  from  the  bar  end,  the  strains  in  the  concrete  and 
steel  become  equal.  In  that  section,  the  force  in  the  steel  is  Fj,  whereas  the 
force  in  the  concrete  is  T  =  F^.,  and  bond  stresses  approach  zero. 

The  two  constants  Cj  and  C2  arc  found  from  the  following  boundary 
conditions: 

(1)  At  X  =  0,  F  =  Fo  (3.14) 


(2)  At  X  =  I ,  the  shear  stress  is  zero,  hence: 


IE 

dx 


0 


(3.15) 


Using  these  two  conditions,  Eq.  3.13  becomes: 
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— ? -  ePx  +  p  - 1 -  e“P*  +  F 

1  +  e-2pi  1  +  e-231  ^ 


(3.16) 


For  high  values  values  of  I  (practically  for  (31  >  4),  the  first  term  has 
a  negligible  contribution  and  Eq.  3. 16  simply  becomes 


F.  e-P*  +  F, 


(3. 17) 


From  Eq.  3.17,  and  using  the  fact  thai  Fq  =  F  +  T  =  F^  +  F^.,  t.he 
following  is  found: 
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(3.18) 

(3. 19) 


where  fj  is  the  stress  in  the  steel,  and  t  is  the  shearing  stress  at  the  inter¬ 
face. 


Yankelevsky  verified  these  expressions  againt  results  of  tests  run  at 
Berkeley  (Fig.  3.7),  and  the  results  compared  very  well  with  his  theory. 
Indeed,  this  approac'.  is  powerful,  and  the  assumptions  made  are  sound. 
However,  the  approach  is  limited  to  the  uncracked  concrete  case. 


3.6  GENERAL  OVERVIEW 


As  far  as  accuracy  is  concerned,  the  new  theory  presented  in  §3.5 
seems  to  be  the  most  accurate,  without  being  inconsistent.  The  analytical 
method  presented  in  §3.2  is  apparently  exact;  it  is,  however,  based  on  an 
analytical  stress-slip  relationship  that  is  probably  not  very  accurate,  nor  is 
it  always  available.  This  fact  introduces  an  inconsistency  between  the 
mathematical  rigidity  of  the  derivation  of  the  bond  distribution,  and  the 
inaccuracy  of  the  stress-slip  relationship,  on  which  the  results  are  based. 
The  numerical  method  presented  in  §3.3  is  the  most  powerful,  while 
practical  and  realistic,  for  it  can  accommodate  a  wide  variety  of  cases, 
including  a  graphical  stress-slip  relationship  (thus  no  mathematical  model 


Fig.  3.7 
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Predicted  strain  along  bar  compared  to  Berkeley  test  results 
(From  Ref.  65). 


Oo  *10  a 


ig.  3.8 


Splitting  forces  and  bond  splitting  failure  with  deformed  bars 
(From  Ref.  30). 
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is  needed)  and  a  cyclic  case  of  loading.  The  assumed  bond  distribution 
method  of  §3.4  is  not  realistic,  and  lacks  consistency  and  soundness. 

As  far  as  practicality  and  ease  of  use  is  concerned,  the  new  theo¬ 
retical  method  in  §3.5  is  the  easiest  to  use  since  it  involves  simple  math¬ 
ematical  expressions. 

The  authors  have  extended  this  new  approach  to  the  case  of  fiber- 
reinforced  cementitious  composites  (see  Chap.  8)  as  one  of  two  different 
methods  of  predicting  the  interfacial  bond  distribution  in  FRC. 

3.7  FACTORS  AFFECTING  BOND  IN  REINFORCED  CONCRETE 
3.7.1  General 

As  was  emphasized  earlier,  it  is  very  important  to  understand  bond  to 
be  able  to  explain  the  behavior  of  composites,  especially  under  tensile  type 
loading.  The  composites  in  question  can  be  anything  typically  represented 
by  a  weak  brittle  matrix  and  a  stronger  ductile  reinforcement.  While  the 
main  focus  of  this  study  pertains  primarily  to  bond  in  fiber-reinforced 
cementitious  composites,  bond  in  reinforced  concrete  and  prestressed 
concrete  cannot  be  ignored,  and  ought  to  be  addressed,  if  briefly.  While  it 
is  true  that  reinforced  and  prestressed  concrete  and  fiber  reinforced  con¬ 
crete  fall  under  the  category  of  cementitious  composites,  they  display  to  a 
large  extent  a  different  behavior,  notably  in  the  case  of  bond  mechanisms. 

Factors  affecting  bond  in  reinforced  concrete  are  numerous.  Some  of 
them  are  obvious,  while  others  are  harder  to  comprehend.  The  problem  of 
bond  in  reinforced  concrete  is  far  from  being  fully  understood,  even  less  so 
for  the  case  of  fiber-reinforced  composites.  Moreover,  there  is  not  always 
a  consensus  between  researchers  as  to  what  factors  affect  bond  and  in  what 
way.  In  this  work,  only  three  factors  affecting  bond  will  be  discussed, 
even  though  more  factors  are  recognized.  Those  three  factors  were  chosen 
because  of  their  importance,  and  because  they  have  been  addressed  in  the 
past.  They  are:  the  texture  of  the  steel  reinforcement  (plain  versus  de¬ 
formed),  the  concrete  cover,  and  the  confinement.  One  section  will  be 
devoted  to  each,  and  finally,  §3.7.5  will  summarize  the  conclusions. 
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3.7.2  Texture  of  The  Steel  Reinforcement 

This  is  the  most  important  and  the  most  significant  factor  affecting 
bond  properties  of  reinforced  concrete,  but  also  of  other  cementitious 
composites.  One  can  sense  the  difference  between  smooth  and  deformed 
reinforcing  bars  in  bond.  It  is  clear  that  deformed  bars  will  display  a  better 
bond  performance  than  smooth  (plain)  bars.  It  is  also  important  to 
understand  that  the  behavior  of  deformed  bars  is  completely  different  from 
that  of  plain  bars. 

The  classical  concept  of  bond  is  that  it  consists  of  adhesion  between 
the  reinforcement  and  the  surrounding  concrete  and  bond  failure  is  caused 
by  the  failure  of  this  adhesion.  Traditional  design  methods  and  code  rules 
concerning  bond  and  anchorage  of  reinforcement  are  based  upon  this 
concept  and  result  in  provisions  for  allowable  bond  stress  values  or 
corresponding  development  lengths  necessary  for  developing  the  full 
design  strength  in  the  bars. 

This  concept  is  originally  based  on  bond  tests,  usually  pull-out  tests, 
with  smooth  reinforcing  bars,  and  in  such  cases  the  concept  corresponds 
relatively  well  to  the  real  function.  For  deformed  bars,  tests  of  this  kind 
and  this  concept  of  bond  are  however  quite  wrong,  especially  for  bond 
failure  in  an  ordinary  beam  case  where  the  bar  is  placed  near  a  surface  with 
normal  concrete  cover.  In  this  case  the  failure  is  preceded  and  finally 
caused  by  longitudinal  cracking  due  to  the  compressive  forces  radiating  out 
with  a  certain  inclination  a  from  the  surface  deformation  lugs  (see  Fig. 
3.8).  These  forces  in  turn  produce  tensile  stresses  in  the  surrounding 
concrete  and  cause  successive  splitting  cracks  following  some  critical 
splitting  surface  line  between  the  reinforcing  bars  and  the  surface  of  the 
concrete  member.  The  final  failure  takes  place  when  these  cracks  have 
propagated  all  along  the  anchorage  zone  along  the  bar.  Depending  on  the 
arrangement  of  the  reinforcement  and  its  location  relative  to  the  concrete 
surface,  the  covering  layer,  etc...,  these  splitting  lines  can  form  different 
patterns  (Fig.  3.9).  This  concept  of  splitting  line  was  introduced  by 
Losberg  and  Olson  [30]  who  found  that  the  bond  force  is  closely  related  to 
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Splitting  line  patterns  in  some  cases  of  reinforcement 
arrangen-f^nt  (From  Ref.  30). 


(a) 


(b) 


ig.  3. 10  -  Single  rib  and  single  concrete  tests  (From  Ref.  31). 
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the  length  u  of  the  splitting  line  (Fig.  3.9)  and  the  concrete  tensile  strength 
(splitting  strength)  fd- 

In  general,  bond  can  be  thought  of  as  the  shearing  stress  or  force 
between  a  bar  and  the  surrounding  concrete.  Bond  is  made  up  of  three 
components: 

(a)  Chemical  adhesion 

(b)  Friction 

(c)  Mechanical  interlock  between  concrete  and  steel 

Bond  between  concrete  and  plain  bars  consists  primarily  of  the  first 
two  elements,  although  there  is  some  mechanical  interfacial  friction  due  to 
the  roughness  of  the  bar  surface.  Deformed  bars,  however,  depend  fore- 
mostly  on  mechanical  interlocking  for  superior  bond  properties.  Never¬ 
theless,  friction  and  chemical  adhesion  are  not  necessarily  negligible  in  the 
case  of  deformed  bars,  but  are  somewhat  secondary. 

To  obtain  a  clear  understanding  of  the  bond  of  a  deformed  bar  in 
concrete,  the  bond  forces  as  well  as  the  slip  will  be  examined.  Unless  the 
strain  in  the  concrete  and  the  steel  is  the  same  and  constant  over  a  length,  a 
deformed  bar  attempts  to  move  or  slip  in  relation  to  the  surrounding 
concrete.  Initially,  chemical  adhesion  combined  with  mechanical  interac¬ 
tion  prevents  significant  slippage.  After  adhesion  is  destroyed,  and  slip 
occurs,  the  ribs  of  a  bar  restrain  this  movement  by  bearing  against  the 
concrete  between  the  ribs;  this  concrete  will  be  called  the  concrete  key 
(Fig.  3. 10).  Friction,  which  would  occur  after  slip  is  induced  in  the  case 
of  plain  bars,  does  not  occur  here  because  of  the  presence  of  the  ribs. 

Slip  of  a  deformed  bar  can  occir  in  two  ways: 

(a)  the  ribs  can  push  the  concrete  away  from  the  bar  (wedging  action) 

(b)  the  ribs  can  crush  the  concrete 

It  was  determined  experimentally  from  pull-out  tests  on  bars  with 
only  a  single  concrete  key  (Fig.  3.10)  that  the  movement  of  the  ribs  is 
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about  the  same  for  all  ribs  with  face  angles  greater  than  40'  (the  face  angle 
of  ribs  is  measured  with  respect  to  the  bar  axis).  These  tests  included 
specimens  with  ribs  having  face  angles  of  90*  (and  even  105')  in  which 
case  the  rib  cannot  push  the  concrete  outward,  hence  preventing  any 
wedging  action).  The  friction  between  the  rib  face  and  the  concrete  is 
evidently  sufficient  to  prevent  relative  movement  at  the  interface  when  face 
angles  are  larger  than  about  40  to  45‘.  It  follows  that  slip  is  due  almost 
entirely  to  the  crushing  of  the  porous  concrete  paste  (mortar)  in  front  of  the 
ribs  if  the  rib  face  angles  are  larger  than  about  40*.  It  is  understandable 
that  crushing  of  the  concrete  does  occur,  since  the  average  bearing  pressure 
exerted  by  these  ribs  is  very  large. 

Bars  with  ribs  having  face  angles  less  than  about  30'  exhibit 
markedly  different  load-slip  relationships  (Fig.  3.10).  Here  the  friction 
between  ihe  itb  face  and  the  concrete  is  not  sufficient  to  prevent  relative 
movement.  Thus,  slip  is  mainly  due  to  the  relative  movement  between  the 
concrete  and  steel  along  the  face  of  the  rib,  and,  to  a  lesser  extent,  to  some 
crushing  of  the  mortar. 

Slip  caused  by  relative  movement  (in  addition  to  that  caused  by 
crushing)  also  occurs  when  the  frictional  properties  of  the  rib  face  are  re¬ 
duced  by  grease.  In  other  words,  bars  with  ribs  having  face  angles  greater 
than  about  45*,  but  where  the  friction  is  impaired,  exhibit  load-slip  curves 
which  are  similar  to  that  for  a  rib  with  flat  rib  face  angle.  The  extent  to 
which  slip  properties  are  affected  by  grease  depends  nevertheless  on  the 
face  angle;  ribs  with  flatter  face  angles  are  understandably  more  affected  by 
poor  frictional  properties. 

For  the  more  common  case  of  good  frictional  properties  and  a  rib 
face  angle  greater  than  40*,  slip  occurs  by  progressive  crushing  of  the 
concrete  paste  structure  in  front  of  the  rib.  This  does  not  appear  to  pro¬ 
duce  significant  wedging  action  until  considerable  crushing  has  occured,  at 
which  time  a  wedge  of  crushed  concrete  (compacted  powder)  becomes 
lodged  in  front  of  the  rib  and  moves  along  with  it.  This,  in  effect,  pro¬ 
duces  a  rib  with  a  face  angle  of  30  to  40*.  Sudi  wedges  uf  cuuipacicd 
concrete  powder  have  been  observed  in  front  of  ribs  that  had  carried  high 
bond  forces;  Rehm  [53]  states  that  crushing  extends  in  front  of  the  rib  for  a 
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length  of  5  to  7  times  the  height  of  the  rib,  although  the  compacted  powder 
that  moves  with  the  rib  extends  at  most  2  times  the  height  of  the  rib. 

It  is  interesting  to  note  (Fig.  3.11)  that  the  bond  modulus  (slope  of 
bond-slip  curve)  upon  reloading  is  considerably  higher  than  that  displayed 
initially.  This  was  observed  by  Rehm  [53].  This  behavior  can  best  be 
explained  by  the  fact  that  during  the  second  loading,  the  rib  bears  against  a 
compacted  nonporous  crushed  concrete,  as  opposed  to  the  porous,  intact 
concrete  during  the  initial  loading. 

It  should  be  pointed  out  here  that  it  is  not  possible  to  produce  ex¬ 
tensive  crushing  and  wedge  action  in  front  of  every  rib  without  having 
transverse  and  longitudinal  cracking,  caused  by  the  forces  produced  by  the 
ribs  bearing  on  the  concrete. 

As  can  be  deduced  by  what  preceded,  the  bond  behavior  of  deformed 
bars  is  in  essence  very  different  from  that  of  smooth  bars,  and  surface 
deformations  improve  the  bond  performance  of  bars.  Also,  it  is  the  angle 
of  the  rib,  rather  than  its  size,  that  dictates  the  general  behavior  of 
deformed  bars.  Many  researchers  were  thus  misled  all  along  when  they 
assumed  that  a  deformed  bar  with  closely  spaced  lugs  has  a  higher  bond 
strength  than  a  similar  one  with  a  larger  rib  spacing. 

3.7.3  Concrete  Cover  Effect 

Another  important  factor  affecting  bond  in  reinforced  concrete  is  the 
thickness  of  concrete  binding  the  steel  bar.  Conceptually,  it  is  easy  to 
visualize  how  the  cover  makes  a  difference  in  the  bond  strength  and  prop¬ 
erties  of  the  interface.  The  concrete  surrounding  a  steel  bar  in  tension  is 
subjected  to  compressive  stresses.  When  the  concrete  cover  is  very  small 
for  instance,  then  the  concrete  might  fail  under  the  compressive  forces  that 
are  acting  on  it.  However,  if  the  concrete  cover  is  very  large,  say  in¬ 
finitely  large,  then  this  does  not  mean  that  the  bond  at  the  interface  of  the 
concrete  and  the  steel  is  infinitely  strong. 

The  problem  has  hence  to  be  viewed  in  the  following  manner.  If  the 
cover  is  too  small,  then  the  concrete  governs  the  bond,  and  the  failure  of 
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the  interface  is  that  of  the  concrete.  Beyond  a  certain  value  for  the  cover, 
the  bond  is  dictated  by  the  strength  of  the  adhesion  at  the  interface.  In 
other  words,  even  though  the  concrete  surrounding  the  steel  will  remain 
intact,  the  bond  between  the  concrete  and  the  steel  will  eventually  fail.  No 
discussion  of  this  threshold  value  can  be  found  in  the  available  literature  on 
bond,  hence  the  weakness  of  most  bond  theories. 

Orangun,  Jirsa,  and  Breen  [46]  derived  the  following  equation  from 
a  non-linear  regression  analysis  of  test  data; 

=  1.2  -i-  3.0;^  +  50^  (3.20) 

where  is  the  bond  stress  capacity  (psi),  c  is  the  smallest  concreic 

cover  or  half  the  bar  spacing,  d^  is  the  bar  diameter,  L  is  the  embedment 
length,  and  f  ^  is  the  concrete  cylinder  strength  (psi). 

The  major  flaw  of  this  equation  is  that  it  predicts  an  infinite  bond 
strength  for  a  steel  reinforcing  bar  embedded  in  a  very  huge  matrix  of  con¬ 
crete  (since  the  cover  then  tends  to  infinity).  It  also  predicts  an  infinitely 
large  bond  strength  for  small  embedment  lengths. 

Kemp  and  Wilhelm  [28]  presented  the  following  equation: 

^bv  =  (^2.64  2.37  ~  j  (3.21) 

where  tbv  the  average  bond  stress  at  first  visible  crack  (psi),  c  is  the 
concrete  cover  (in),  db  is  the  bar  diameter  (in),  and  f  ^  is  the  concrete 
compressive  strength  (psi). 

This  concept,  directly  relating  the  bond  strength  (or  average  bond)  to 
the  concrete  cover  lacks  in  that  it  only  works  for  small  covers.  For  larger 
covers,  it  is  the  adhesion  between  the  steel  and  the  concrete  that  governs 
the  strength  of  the  interface 

In  conclusion,  the  concrete  cover  influences  bond  in  a  sense  that  it 
needs  to  be  large  enough  not  to  create  problems  at  the  interface  (namely 
causing  a  splitting  mode  of  failure).  However,  the  equations  presented  in 
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this  section  cannot  be  used  for  all  cases,  but  need  to  be  modified  to  ac¬ 
commodate  the  fact  that  beyond  a  certain  value,  concrete  cover  has  hardly 
any  effect  on  bond. 

3.7.4  Effect  of  Confinement  on  Bond 

The  confinement  effect  can  be  not  so  different  from  that  of  the 
concrete  cover.  When  an  external  pressure  is  applied  to  the  concrete  sur¬ 
rounding  the  steel  subjected  to  tension,  the  compressive  capacity  of  the 
concrete  increases  by  up  to  20%  of  its  original  value,  which  makes  the 
interface  stronger  in  resisting  splitting.  In  this  respect,  the  effect  of  con¬ 
finement  is  exactly  the  same  as  that  of  the  concrete  cover. 

However,  the  confinement  reinforces  the  bond  capacity  of  the  inter¬ 
face  no  matter  what  the  concrete  cover  is,  even  when  the  mode  of  failure  is 
not  of  the  splitting  type.  Indeed,  after  the  adhesion  between  the  concrete 
and  the  steel  breaks,  frictional  stresses  develop  along  the  interface.  These 
stresses  are  directly  related  to  the  normal  stresses  acting  on  the  steel.  This 
means  that  for  the  same  slip,  more  frictional  stresses  would  be  needed  to 
develop  along  the  interface,  which  implies  that  the  tensile  capacity  of  the 
system  is  increased. 

Finally,  it  can  be  stated  that  confinement  of  the  concrete  around  the 
steel  subjected  to  tension  is  helpful  in  increasing  the  bond  capacity  of  the 
interface  and  the  tensile  strength  of  the  system.  Such  a  confinement  is 
usually  provided  by  tranversc  steel  reinforcement  which,  by  preventing  the 
concrete  from  expanding,  induces  a  confining  pressure. 


CHAPTER  IV 


OVERVIEW  OF  BOND  IN  FIBER-REINFORCED  CONCRETE 


4.1  GENERAL 

Bond  in  fiber-reinforced  concrete  is  a  discipline  in  itself.  However, 
it  is  very  useful  to  look  at  bond  in  reinforced  concrete  before  attempting  to 
solve  the  bond  problem  with  fibers  since  fibers  feature  a  completely 
different  behavior  than  conventional  reinforcement. 

The  main  difference  between  fiber-reinforced  concrete  and 
conventional  reinforced  concrete  lies  in  the  fact  that,  unlike  reinforcing 
bars,  fibers  are  not  continuous.  Moreover,  since  many  fibers  are  near  each 
other,  their  interaction  must  be  considered.  In  the  more  extreme  case  of 
Slurry  Infiltrated  Fiber  Concrete  (SIFCON),  interaction  may  include  fiber 
entanglement.  The  fiber  entanglement  modifies  the  behavior  of  the  whole 
composite.  A  long  and  single  steel  fiber  embedded  in  concrete  and 
subjected  to  a  direct  tensile  load  (typical  pull-out  test),  will  behave  like  a 
steel  reinforcing  bar  in  the  same  condition  and  undergoing  the  same  kind  of 
loading.  However,  even  in  the  case  of  normal  fiber-reinforced  concrete 
(as  opposed  to  SIFCON),  a  single  fiber,  in  the  presence  of  the  many  other 
discontinuous  fibers  in  a  concrete  matrix,  behaves  differently.  The  very 
fact  that  to  date,  no  complete  analytical  model  of  the  behavior  of 
discontinuous  fibers  in  a  cementitious  matrix  has  been  developed  yet, 
shows  how  intricate  the  problem  of  bond  in  fiber-reinforced  concrete  is. 

The  analytical  work  in  fiber-reinforced  concrete  is  not  the  only 
problem  a  researcher  might  encounter.  Indeed,  the  experimental  work  can 
be  at  least  equally  as  streneous.  The  best  example  is  the  problem  of 
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measuring  strains  and/or  displacements  along  the  interface  of  a  steel  fiber 
with  concrete.  It  seems  almost  impossible  to  make  any  actual  measure¬ 
ment,  mainly  because  of  the  size  of  the  fiber,  not  to  mention  the  interaction 
among  fibers.  One  possible  way  of  solving  the  problem  will  be  to  use  the 
following  procedure: 

(1)  To  analytically  model  the  behavior  of  a  single  fiber  in  a  concrete 
matrix. 

(2)  To  use  results  of  experiments  to  evaluate  some  of  the  parameters 
involved  in  the  theory. 

(3)  To  generalize  the  theory  to  the  case  of  multiple  fibers,  and  possibly 
SIFCON  in  a  semi-empirical  way. 

(4)  To  use  the  model  thus  obtained  in  conjunction  with  other  theories 
to  predict  the  overall  behavior  of  the  specimen  under  tensile 
loading. 

This  study  will  only  address  the  first  part  of  the  problem,  i.e.  the 
analytical  model  for  a  single  fiber  will  be  developed,  and  then  experimental 
data  would  be  used  to  substantiate  it. 

The  remaining  part  of  this  chapter  will  deal  with  the  differences 
between  fibers  and  rebars  in  bond,  as  well  as  the  properties  of  bond  in 
fiber-reinforced  concrete. 

4.2  TYPES  OF  BOND  IN  FIBER-REINFORCED  CONCRETE 

Depending  on  the  type  of  stress  transferred  accross  an  interface,  two 
types  of  bond  are  identified:  the  tensile  bond  and  the  shear  bond. 

The  tensile  bond  resists  displacement  caused  by  forces  acting  per¬ 
pendicular  to  the  interface.  This  kind  of  bond  is  of  no  major  interest  in 
this  study,  and  will  thus  be  overlooked  in  the  remaining  part  of  this  report. 

The  shear  bond  controls  the  transfer  of  stresses  parallel  to  the  lon¬ 
gitudinal  axis  of  the  reinforcement  (rebar  or  fiber).  In  an  uncracked  com- 


posite,  the  shear  bond  transfers  stresses  from  the  matrix  back  into  the  re¬ 
inforcement.  When  the  matrix  cracks  and  the  load  is  carried  by  the  fibers 
bridging  the  crack,  the  shear  bond  enables  the  load  to  be  transferred  back 
into  the  uncracked  parts  of  the  matrix.  Shear  bond  also  resists  the  pull-out 
of  the  reinforcement  from  the  matrix  and  is  therefore  one  of  the  main 
factors  influencing  the  mechanism  and  mode  of  failure  of  a  composite. 
Two  types  of  shear  bond  are  capable  of  transferring  stresses  acting  along 
the  interface:  the  elastic  shear  bond  which  includes  adhesion  if  any,  and  the 
frictional  shear  bond.  When  the  elastic  shear  stress  at  the  interface 
between  the  reinforcement  and  the  matrix  exceeds  the  bond  strength  of  the 
interface,  bond  becomes  frictional  in  nature  and  the  adhesion  is  broken. 

Elastic  Shear  Bond:  If  an  elastic  shear  bond  exists  at  an  interface  and 
the  shear  stress  does  not  exceed  the  strength  of  the  bond,  then  the  shear 
stresses  are  directly  related  to  the  relative  slip  between  the  reinforcement 
and  the  matrix.  Bartos  [4]  stated  that  when  elastic  bond  is  prevailing,  "the 
longitudinal  displacements  of  the  fibre  and  matrix  at  the  interface  remain 
compatible."  A  compatibility  of  displacements  would  imply  zero  slip, 
which  would  then  mean  that  no  shear  stresses  are  developed  at  the 
interface,  hence  implying  an  infinitely  large  bond  modulus. 

Frictional  Shear  Bond:  The  frictional  shear  bond  resists  displace¬ 
ments  along  the  interface  parallel  to  the  length  of  the  fiber,  and  may  or  may 
not  depend  on  the  local  slip.  If  the  frictional  bond  is  perfectly  plastic 
(theoretical  case),  then  the  stress  will  be  completely  independent  of  the 
slip.  The  frictional  bond  is  characterized  by  an  initial  value  of  the  bond, 
referred  to  as  the  bond  strength  which  has  to  be  exceeded  for  the 

elastic  bond  to  break  and  for  the  frictional  stresses  to  take  place. 

4. 3  BOND  IN  FIBER  REINFORCED  CONCRETE 
4.3.1  Uncracked  Composite 

When  a  fiber-reinforced  specimen  is  subjected  to  a  tensile  load,  the 
tensile  stresses  are  transmitted  from  the  matrix  to  the  fibers  and  the  fibers 
undergo  tension.  Depending  on  the  strength  of  the  bond  between  the  fiber 
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and  the  matrix,  the  bond  may  or  may  not  fail  under  a  given  load.  It  is 
possible  that  fhe  following  happens  before  ciacking  of  the  matrix: 

1-  The  load  is  so  low  (or  the  chemical  adhesion  is  so  strong)  that  the 
adhesive  bond  does  not  break,  and  there  is  no  relative  displacement 
between  the  fiber  and  the  matrix. 

2-  The  load  breaks  the  chemical  adhesion  between  the  matrix  and  the 
fiber;  however  the  bond  strength  of  the  fiber  is  not  attained,  so  that  a 
condition  of  elastic  bond  prevails  along  the  interface. 


4.3.'*.  Cracked  Composite 

Bond  remains  a  very  important  factor  in  determining  the  properties 
of  a  composite  even  in  its  post-cracking  stage,  but  its  function  changes. 
The  change  is  most  significant  in  the  case  of  discontinuous  fibers.  In  the 
pre-cracking  stage,  when  the  elastic  shear  bond  was  intact,  the  maximum 
shear  flow  occured  at  the  ends  of  the  fibers.  The  stresses  were  transferred 
from  the  matrix  to  the  fibers.  In  the  post-cracking  stage,  the  fibers 
bridging  a  crack  are  pulled  out  and  the  interfacial  shear  stress  changes;  the 
maximum  stress  now  occurs  at  the  point  where  the  fibers  enter  the  matrix 
i.e  at  the  crack.  If  the  frictional  shear  bond  is  perfectly  plastic  before  the 
start  of  the  cracking,  it  is  likely  to  stay  this  way  in  the  post  cracking  stage. 

4.4  FIBER  DEBONDING 

Fundamental  information  on  the  nature  and  strength  of  the  interfacial 
bond  between  the  fiber  and  the  matrix  in  a  fiber-reinforced  composite 
material  is  usually  obtained  using  the  single-fiber  concentric  pull-out  test 
specimen  and  configuration  shown  schematically  in  Fig.  4.1.  Stress  dis¬ 
tributions  in  this  specimen  are  similar  to  those  within  and  on  the  surface  of 
a  fiber  normal  to  and  bridging  a  matrix  crack  in  a  composite.  The  sim¬ 
ilarity  holds,  assuming  that  the  failure  is  by  matrix  tensile  fracture  fol¬ 
lowed  by  fiber-matrix  debonding  and  fiber  pull-out,  i.e.  the  common  fail¬ 
ure  mode  for  discontinuous  fiber  composites  with  a  brittle  or  clastic  matrix. 
However,  the  values  and  the  exact  stress  distribution  are  obviously 
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Fig.  4. 1 


Typical  single-fiber  concentric  pull-out  test  specimen. 
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different  for  a  single  fiber  does  not  behave  as  it  does  in  the  presence  of  a 
large  number  of  neighboring  fibers. 

Resistance  to  the  debonding  and  pull-out  processes  is  primarily  a 
function  of  the  fiber-matrix  interfacial  bond  shear  strength  and  the  inter¬ 
facial  bond  area.  Theoretical  analyses  of  this  resistance  were  initially 
developed  for  ductile  or  plastically  deforming  matrices  and  a  uniform 
interfacial  shear  stress  along  the  length  of  the  embedded  fiber  was  as¬ 
sumed.  Under  this  assumption,  sudden  failure  of  the  interfacial  bond 
occurs  when  its  shear  strength  is  reached  and  the  debonding  load,  that  is 
the  peak  load  observed  in  the  pull-out  test,  is  directly  proportional  to  the 
embedded  fiber  length. 

However,  theoretical  analyses  of  the  stress  conditions  in  a  pull-out 
test  specimen  with  a  brittle  or  elastic  matrix  show  that  the  distribution  of 
shear  stress  along  the  interface  is  not  uniform,  and  that  the  relationship 
between  the  debonding  or  the  pull-out  load  and  the  embedded  length  of 
fiber  in  the  specimen  is  more  complex  than  indicated  above.  Furthermore, 
the  variation  of  the  peak  pull-out  load  with  the  embedded  length  is  nonlin¬ 
ear,  and  consequently,  it  cannot  be  inferred  that  the  true  value  of  the 
interfacial  bond  strength  can  be  obtained  simply  from  this  load  and  the 
nominal  interfacial  area.  In  fact,  as  will  be  shown  in  Chap.  6,  the  peak 
pull-out  load  is  related  to  the  frictional  bond  shear  stress  as  well  as  the 
bond  strength. 

What  preceded  indicates  that  even  though  the  pull-out  load  is  not  lin¬ 
early  related  to  the  embedment  length,  the  embedded  length  has  to  be 
included  in  any  bond  model. 

4.5  FACTORS  AFFECTING  BOND  IN  FIBER-REINFORCED  CONCRETE 

The  problem  of  bond  in  fiber-reinforced  concrete  is  hardly  under¬ 
stood  at  all.  All  ideas  presented  so  far  are  incomplete.  Experimental 
invetigations  have  been  carried  out  [10,  19,  35, .40,  50,  59],  but  as  will  be 
seen  in  what  follows,  corresponding  results  were  not  always  consistent. 
Following  is  a  brief  presentation  of  some  of  the  findings. 
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4.5.1  Experimental  Investigation 

Naaman  and  Shah  [40]  reported  that  the  bond  efficiency  in  a  pull-out 
test  of  steel  fibers  inclined  with  respect  to  the  line  of  stress  is  at  least  as 
good  as  that  of  fibers  parallel  to  the  direction  of  stress.  They  also  found 
that  fiber  efficiency  decreases  with  the  number  of  fibers  being  pulled  out 
simultaneously  from  the  same  area. 

Gokoz  and  Naaman  [16]  noticed  that  for  steel  fibers  where  a  pull-out 
mode  of  failure  was  observed,  the  postpeak  response  is  almost  insensitive 
to  the  loading  velocity.  They  also  reported  that  the  peak  pull-out  load, 
which  is  indicative  of  the  elastic  bond  at  the  interface,  is  strain-rate 
sensitive. 

It  was  reported  by  Gray  and  Johnston  [19]  that  the  direction  of 
casting  has  a  substantial  influence  on  the  bond  strength.  They  stated  that 
vertically-cast  specimens  have  a  higher  interfacial  bond  strength  than 
horizontally  cast  specimens.  They  also  found  out  that  an  increase  in  the 
sand-cement  ratio  in  the  mortar  matrix  has  contrasting  effects  on  the 
strength  of  the  interfacial  bond  in  the  vertically  and  horizontally-cast 
specimens.  For  the  vertically-cast  specimens,  there  is  a  decrease  in  the 
apparent  interfacial  bond  strength.  They  explained  that  "this  reduction  in 
bond  strength  is  almost  certainly  a  result  of  a  change  in  the  relative  pro¬ 
portions  of  cement  paste  and  aggregate  in  the  matrix  with  a  change  in  the 
sand-cement  ratio."  For  the  horizontally-cast  specimens  there  is  an  in¬ 
crease  in  the  apparent  bond  strength  with  an  increase  in  the  sand-cement 
ratio.  This  was  explained  by  what  the  authors  called  the  water  gain,  i.e. 
"the  accumulation  of  water  under  the  embedded  fiber  due  to  the  tendency  of 
the  bleed  water  in  the  mortar  matrix  to  rise  during  and  after  consolidation." 
Water  gain  is  known  to  have  a  weakening  effect  on  the  cement  paste- 
aggregate  Dond  in  plain  concretes. 

It  was  also  noted  by  Gray  and  Johnston  [19]  that  the  average  inter¬ 
facial  bond  strength  increases  slightly,  but  not  significantly,  with  an  in¬ 
crease  in  the  fiber  withdrawal  rate,  which  is  not  in  contradiction  with  the 


findings  of  Gokoz  and  Naaman  [161  who  stated  that  bond  in  fibers  is  rate- 
independent. 

Maage  [35]  found  out  that  the  bond  properties  between  steel  fibers 
and  cement-based  matrices  are  of  mechanical  nature,  where  the  anchoring 
effect  is  essentially  more  important  than  the  adhesive  effect.  He  also  said 
that  the  mean  pull-out  load  per  fiber  is  unaffected  by  the  number  of  fibers. 
This  finding  is  questionable,  since  this  would  imply  that  N  fibers  have  N 
times  more  strength  than  one  single  fiber,  which  is  not  the  case  as  reported 
by  Naaman  and  Shah  [40].  Maage  ascerts  that,  "based  on  the  weakest  link 
theory,  it  should  be  reasonable  that  the  pull-out  load  per  fiber  would 
decrease  when  the  number  of  fibers  in  the  area  increased."  This  of  course 
is  an  indication  of  the  variability  often  encountered  in  fiber  bond  tests. 

Pinchin  [50],  who  carried  out  tests  on  wires,  stated  that  a  com¬ 
paction  of  the  concrete  surrounding  a  fiber  woula  actually  increase  the 
frictional  bond.  This  property  is  discussed  in  Chap.  3  for  the  case  of 
reinforced  concrete.  He  also  determined  that  the  pull-out  load  increases 
linearly  with  the  confinement.  Another  interesting  finding  is  that  deformed 
fibers  are  less  sensitive  to  pressure  than  plain  (straight)  fibers.  He  also 
showed  analytically  that  the  pull-out  load  is  proportional  to  the  fiber-matrix 
misfit,  which  he  defined  as  the  difference  between  the  radius  of  the  wire 
and  that  of  the  hole  in  the  matrix  in  the  absence  of  the  wire. 

Burakiewicz  [10]  reported  that  the  shape  of  the  load-displacement 
curve  depends  on  the  fiber  type,  and  that  hooked  fibers  show  smaller 
scatter  in  bond  strength  than  all  the  other  types  he  tested.  He  also  ob¬ 
served  that  the  pull-out  of  hooked  and  indented  fibers  requires  more  energy 
than  that  of  plain  fibers,  which  implies  that  deformed  fibers  have  a  higher 
bond  strength  than  straight  fibers.  He  also  stated  that  there  is  no 
significant  influence  on  the  bond  strength  of  vibration  and  orientation  of 
fibers  during  setting  and  hardening  of  the  matrix,  which  is  in  contradiction 
with  the  findings  of  Gray  and  Johnston  [19].  He  finally  noted  that  bond 
strength  seeir-  to  depend  on  the  rate  of  loading  of  the  matrix,  which 
contradicts  the  findings  of  Gokoz  and  Naaman  [16]. 
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Gopalaratnam  and  Abu-Mathkour  [66]  investigated  experimentally  the 
effect  of  the  following  three  parameters  on  the  fiber  pull-out 
characteristics:  fiber  embedment  length,  fiber  diameter  and  matrix  quality. 
They  observed  that  the  average  bond  strength  is  inversely  related  to  the 
embedment  length  and  that  the  average  bond  strength  of  an  interface 
increases  with  an  increase  in  fiber  diameter.  They  also  observed  that  the 
peak  pull-out  load  for  identical  fibers  being  pulled  out  from  matrices  of 
increasing  compressive  strength  does  not  necessarily  increase.  Their 
reasoning  is  that  the  frictional  bond  strength  may  be  unrelated  to  the  matrix 
''ompressive  strength. 

4.5.2  Analytical  investigation 

Few  studies  have  addressed  the  analytical  modeling  aspects  of  bond 
in  cement  composites.  Most  often  in  such  cases,  some  aspects  of  the 
problem  were  overlooked. 

Romualdi  and  Batson  [54]  assumed  that  "in  the  absence  of  a  crack, 
the  strain  in  the  concrete  due  to  remote  tension  is  equal  to  the  strain  in  the 
wires,"  hence  implying  zero  slip  between  the  wires  and  the  surrounding 
concrete. 

Aveston,  Mercer,  and  Sillwood  [3]  assumed  frictional  shear  bond 
with  no  debonding,  hence  ignoring  the  elastic  bond  that  develops  for  shear 
stresses  not  exceeding  the  bond  capacity  of  the  interface. 

Kar  and  Pal  [27],  on  the  other  hand,  assumed  a  linear  shear  stress 
distribution  along  the  fiber,  with  a  maximum  at  the  point  where  a  crack 
intersects  the  fiber. 

Swamy  and  Mangat  [61]  made  a  similar  assumption  concerning  the 
linearity  of  the  shear  stress  distribution  at  the  interface.  However,  it  can 
be  shown  that  a  linear  bond-slip  relationship  does  not  lead  to  a  linear  bond 
shear  stress  distribution  at  the  interface. 
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Larson  and  Bayasi  [69]  developed  a  special  test  to  measure  bond 
strength  of  carbon  fiber  embedded  in  a  cementitious  matrix.  They  observed 
that  bond  strength  increases  with  the  water-binder  ratio  of  the  matrix  which 
is  attributed  to  an  improvement  in  workability.  They  found  out  that  the 
bond  strength  decreases  with  addition  of  microsilica  and  increases 
drastically  with  the  addition  of  latex. 

Parameswaran  and  Rajagopalan  [48]  recognized  the  fact  that  there  is 
no  strain  compatibility  between  the  matrix  and  the  reinforcement  once 
debonding  occurs  but  didnot  extend  this  assumption  to  the  pre-debonding 
phase.  Indeed,  at  any  point,  the  existance  of  a  mechanical  shear  stress 
must  always  be  accompanied  with  a  corresponding  slip. 

More  recently,  a  study  by  Gopalaratnam  and  Shah  [17]  addressed  the 
problem  of  the  tensile  response  of  steel  fiber-reinforced  mortar.  Bond  was 
incorporated  in  the  study  of  the  pull-out  behavior  of  fibers.  However,  in 
an  attempt  to  apply  shear  lag  theory  to  the  problem  of  fiber  pull-out,  they 
ignored  the  normal  tensile  stresses  in  the  matrix,  hence  creating  a  serious 
flaw  in  the  analytical  part  of  the  study. 

Nammur  ct  al.  [43]  performed  a  fundamental  study  of  the  pull-out 
problem  and  related  the  pull-out  curve  to  the  interfacial  bond-slip  rela¬ 
tionship.  They  derived  a  method  of  analytically  predicting  the  bond-slip 
curve  from  an  experimental  pull-out  curve.  The  method  is  presented  with 
some  improvement  in  detail  in  Chap.  5. 

Mandel,  Wei  and  Said  [67]  developed  a  procedure  to  determine  the 
coefficient  of  adhesion,  and  the  energy  release  rate  per  unit  area  of  crack 
surface  at  the  crack  front  required  for  unstable  growth  of  a  crack  along  the 
fiber-matrix  interface.  This  procedure  uses  the  results  of  the  pull-out  test. 

Li  and  Wang  [70]  developed  a  theoritical  model  to  predict  the  load- 
crack  seperation  relationship  for  synthetic  fiber  reinforced  concrete.  The 
distinctive  feature  of  the  model  lies  in  the  recognition  that  the  fiber  matrix 
bond  strength  at  any  point  on  the  fiber  is  a  function  of  the  slippage 
accumulation  at  that  fiber  segment. 


Shah  and  Jenq  [71]  investigated  the  interfacial  bond  properties 
between  fibers  and  matrix  through  analysis  of  pull-out  tests.  The  bonding 
between  fibers  and  matrix  was  assumed  to  be  perfect  before  the  pull-out 
load  is  applied.  Griffith  energy  criterion  was  used  to  govern  the  crack 
propagation  at  the  interfacial  region  where  debonding  begins.  They  also 
developed  a  fracture  model  that  predicts  the  progressive  failure  of  the 
interfacial  bond.  However,  the  proposed  model  doesnot  account  for 
friction  in  the  debonded  zone. 

Gao,  Mai,  and  Cotterell  [72]  studied  the  interface  between  the  matrix 
and  the  fiber  with  a  new  fracture  mechanics-based  debonding  criterion  by 
including  the  friction  in  the  debonded  region.  In  their  study  they  defined 
several  basic  concepts  such  as  the  threshold  value  of  the  interface  pressure, 
friction  zone  as  well  as  strong  and  weak  fibers. 

Gopalaratnam  and  Chang  [73]  formulated  the  pull-out  problem  in  one 
dimension  with  a  view  to  focus  on  the  nonlinear  interfacial  response.  In 
general,  their  work  was  narrowed  to  the  effect  of  the  softening  of  the 
interface  on  the  pull-out  strength,  composite  ductility,  energy  absorption 
and  the  stability  of  the  debonding  process. 

What  preceded  gives  an  indication  about  what  could  affect  the  bond 
performance  of  a  fiber  in  a  concrete  matrix.  It  also  shows  how  far  away 
the  current  state  of  the  art  of  bond  in  fiber-reinforced  concrete  is  from 
being  scientifically  reliable  and  accurate.  As  was  mentioned  earlier,  this  is 
partly  due  to  the  difficulties  involved  in  performing  experiments  and 
making  measurements  when  small  discontinuous  fibers  are  used. 


CHAPTER  V 


ANALYTICAL  STUDY  OF  PULL-OUT  PROBLEM 


5.  1  OBJECTIVE 

The  main  objective  of  the  experimental  part  of  this  study  is  to 
generate  bond  shear  stress-slip  curves  from  experimentally  obtained  pull¬ 
out  curves  for  various  fibers  and  matrix  types.  Furthermore,  an  attempt  is 
to  be  made  at  experimentally  characterizing  the  thus-obtained  bond  shear 
stress-slip  curves. 

An  analytical  study  of  the  mechanics  of  the  pull-out  test  is  first  pre¬ 
sented,  so  as  to  make  it  possible  to  relate  a  typical  pull-out  curve  (P-A)  to 
a  bond  shear  stress-slip  (x-S)  relationship.  A  bond-slip  relationship  is 

first  assumed,  and  a  pull-out  curve  is  hence  derived.  The  problem  of 
obtaining  (x-S)  from  a  given  (P-A)  curve  would  be  the  reciprocal  to  the 

above  mentioned  problem. 

5.2  STATEMENT  OF  THE  PROBLEM 

The  analysis  that  follows  pertains  to  a  pull-out  test  whereby  a  tensile 
force  P  is  applied  to  the  tip  of  a  fiber  embedded  over  a  length  f  in  a 
cementitious  body  (Fig.  5.1). 

Initially,  it  is  assumed  that  the  relationship  between  the  bond  shear 
stress  at  the  interface  between  the  fiber  and  the  matrix,  and  the  relative  slip 
between  the  same  components  is  as  shown  in  Fig.  5.2.  A  more  accurate 
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Fig.  5.  1 


Pull-out  test  confi 


Fig.  5.2  -  Assumed  bond  she 
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relationship  will  be  introduced  later  in  this  chapter.  Thus,  it  is  assumed 
that  the  bond-slip  curve  is  linear  elastic  up  to  the  point  where  the  bond 
strength  of  the  interface  is  reached,  beyond  which  stage  purely 

frictional  conditions  prevail,  with  a  constant  frictional  shear  stress  equal  to 
Tf.  It  is  also  assumed  that  tf  cannot  exceed  Tmax-  value  of  Tf 

ranges  between  zero  and  While  it  is  understood  that  a  constant  bond 

stress  equal  to  Xf  cannot  be  sustained  forever,  it  is  furthermore  assumed 
that  slips  causing  the  bond  stress  to  diminish  significantly  beyond  Xf  will 
not  be  attained. 


5.3  MATHEMATICAL  DERIVATION 
5.3.1  Basic  Equations 

The  free  body  diagram  of  an  infinitesimal  segment  of  fiber  (Fig.  5.3) 
leads  to  the  following  equation: 

dF  -  X  y  dx  =  0  (5.1) 

where  F  is  the  local  force  in  the  fiber  at  a  distance  x  from  the  embedded 
end  of  the  fiber,  y  is  the  perimeter  of  the  fiber,  and  x  is  the  local  shear 
stress  at  the  interface  between  the  fiber  and  the  matrix. 


Eq.  5. 1  can  be  rewritten  as: 


dx 


X  y 


(5.2) 


The  product  x  y  is  the  shear  force  per  unit  length,  or  the  shear  flow  t 
at  the  interface.  Hence: 

t  =  X  y  (5.3) 


Furthermore,  the  local  tensile  force  F  in  the  fiber  can  be  related  to 
the  local  strain  in  the  fiber  Ef  thru  the  following  equation: 


F 


Af  Ef  Ef 


(5.4) 


-64- 


where  Af  and  Ef  arc  respectively  the  area  and  the  elastic  modulus  of  the 
fiber,  and  Ef  is  the  local  strain  in  the  fiber. 

Moreover,  in  the  elastic  region,  the  local  shear  stress  X  can  be  related 
to  the  local  slip  S  through  the  bond  modulus  k  as  follows: 

X  =  K  S  (5.5) 

=  K  (5f  -  5m)  (5.6) 

where  5f  and  5^  arc  the  local  displacement  in  the  fiber  and  the  matrix  re¬ 
spectively. 

To  satisfy  static  equilibrium,  the  total  force  P  applied  at  the  free  tip 
of  the  fiber  has  to  be  resisted  at  any  section  by  the  local  force  in  the  fiber 
F,  as  well  as  the  local  force  in  the  matrix  T,  or: 

P  =  F  -t- T  (5.7) 

=  Af  Ef  Ef  +  Am  Em  Em  (5-8) 

where  Am  and  Em  are  respectively  the  area  and  the  elastic  modulus  of  the 
matrix,  and  Em  is  the  local  strain  in  the  matrix. 

Differentiating  Eqs.  5.2  and  5.6  leads  to: 


d^F 

dx2 

dx 

dx 


y  dx 
dx 


dx  ) 


(5.9) 

(5.10) 


Combining  Eqs.  5  9  and  5. 10,  and  recognizing  the  fact  that  the  first 
derivative  of  the  local  displacement  function  is  the  local  strain,  it  can  be 
proven  that: 


d2F 

dx2 


V  <  (Ef  -  Em) 


(5. 11) 


From  Eq.  7.8,  Em  can  be  expressed  as  a  function  of  Ef  as  in: 
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1 


-m 


‘tn  *^m 


p—  (P  -  Af  Ef  Ef) 

*-»m 


(5.12) 


Replacing  in  Eq.  5.11  by  the  right  hand  side  of  Eq.  5.  12  leads  to: 
d^F  V  K 


dx2 


where; 


[  Am  Em  Ef  -  (P  -  Af  Ef  Ef)] 


■KP  +  K[Am  Em  +  AfEf]  Ef 


-K  P  +  K  [  Am  Em  +  Af  Ef] 


-K  P  +  K  Q  F(x) 


F(x) 
Af  Ef 


(5.13) 

(5.  U) 

(5. 15) 

(5. 16) 


and 


y  1C 
Am  ^m 


Am  +  Af  Ef 
Af  Ef 

A  F 

1  +  >  gm 

Af  Ef 


(5.17) 


(5.18) 


Eq.  7.16  is  a  second  degree  differe'’'’al  equation  in  F  of  the  form; 
d2F 


dx2 


-  F  =  -K  P 


(5.19) 


in  which 


\  =  VKQ 

The  solution  to  this  differential  equation  is  of  the  form; 
F(x)  =  A’  e^*  +  B’  +  Q 


(5.20) 


(5.21) 


The  unknown  coefficients  A'  and  B'  are  determined  from  the  follow¬ 
ing  boundary  conditions; 
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1)  F(0)  =  0 

2)  F(0  =  P 


(5.22) 

(5.23) 


Eqs.  5.22  and  5.23  are  equivalent  to  the  following  two  equations  in 
A'  and  B' : 


A'  +  B’  +  ^  =  0 


A’  e^^  +  B'  e-^^  +  ^  =  P 


(5.24) 

(5.25) 


The  solution  to  the  system  of  equations  described  in  Eqs.  5.24  and 


5.25  is: 


A' 


and  B' 


I  -  e-2>^^ 
P 


(1-  ^)  e-^^  +  ^  e-2^^ 


_1_ 

Q 


1  -  e-2^^  *  Q 


A'  B' 

Introducing  the  coefficients  A  and  B,  equal  to  -p-  and  -p- 


respectively,  we  get: 
A 


1 


and  B 


1  -  e-2>^f 
1 

1  -  e-2>^^ 


(1-  5)  e 


Q 


>.r  +  J- 

Q  '  J 

(5.26) 

(5.27) 

Using  these  expressions  for  A  and  B,  the  force  F(x)  in  the  fiber,  and 
the  intcrfacial  shear  flow  t(x)  can  be  respectively  expressed  as: 


F(x) 

t(x) 


P  ^A  e^*  +  B  e-^*  +  ^ 
dx 

P  X  (a  e^*  -  B  c-^*) 


(5.28) 


(5.29) 
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As  can  be  seen  from  Eq.  5.29,  the  interfacial  shear  flow  t  is  a  direct 
function  of  the  load  P.  Furthermore,  the  shear  flow  is  a  direct  function  of 
the  shear  stress  as  in  Eq.  5.3.  Also,  the  maximum  shear  stress  (or  flow) 
for  a  given  load  will  always  occur  at  the  point  where  the  fiber  penetrates 
the  matrix,  or  at  x  =f . 


5.3.2  Critical  Force 


Given  the  bond-slip  relationship,  there  will  be  a  critical  force  Peru 
that  will  induce  a  shear  stress  at  x  =  f  equal  to  tmax’  o'"  ^  shear  flow  equal 
to  tmax  where: 


tfnax  “  ’^max  V 


(5.30) 


To  find  Peril’  the  shear  flow  at  x  =  f,  or  t(0  has  first  to  be 


t(0  = 

p  X  (a  ( 

-  B  e 

■^0 

P 

1  -  e'2X 

'•5) 

1  (1  +  e*2X^)-t-  j 

[5: 

1  2e-^^j 

Setting  t(0  equal  to 

^max 

yields  the  value  of  F 

> 

cm 

j;  thus; 

tmax 

r 

1  -  e-2^‘ 

\ 

Peril  ~ 

X 

(1  +  e'2Xf)+  j 

[5] 

2e'^^ 

evaluated; 

(5.31) 

(5.32) 


(5.33) 


When  a  force  P  ^  Pcrit>  elastic  bond  conditions  prevail  at  the  inter¬ 
face,  and  no  debonding  occurs;  in  other  words,  the  fiber  remains  fully 
bonded  to  the  surrounding  matrix. 


5.3.3  Elastic  Displacement 

The  displacement  of  the  free  end  of  the  fiber  (at  x  =  f)  can  be  eval¬ 
uated  using  the  following  equation: 
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(5.34) 


To  carry  out  the  integral  in  Eq.  5.34,  F(x)  is  replaced  by  the  right 
hand  side  of  Eq.  5.28.  A  is  found  to  be  as  follows: 


A 


P  (Q  -  2)  n  -  e-^^N 
A„  E„  (l  + 


(5.35) 


As  can  be  seen  from  Eq.  5.35,  there  is  a  linear  relationship  between 
the  applied  force  P  and  the  displacement  of  the  fiber  free  end.  Moreover, 
in  a  pull-out  curve,  the  slope  of  the  ascending  linear  portion  of  the  curve  is 
given  by: 


^  Efn  1  +  e~^^ 

(Q  -  2)  I  - 


6) 


5.3.4  Debonding  Zone 

When  the  applied  force  P  exceeds  Pcrit.  ^  region  identified  as  "zone 
of  debonding"  w’’l  develop  and  grow  as  the  applied  force  P  increases.  In 
other  words,  two  intcrfacial  zones  will  adjacently  coexist,  one  that  is 
bonded  and  one  that  is  debonded,  i.e.  one  where  debonding  has  occurred 
because  the  shear  stresses  have  exceeded  Xmax-  Fig.  5.4  shows  typical 
bond  shear  stress  distribution  for  cases  with  different  values  of  P.  The 
forces  resisted  by  these  two  individual  zones  will  be  identified  as  bonded 
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force  P|j  and  debonded  force  respectively.  To  satisfy  static  equilibrium, 
it  can  be  inferred  that  for  any  load  P  >  Pcr^, 

P  =  Pb  +  Pd  (5.37) 

Based  on  the  assumed  bond  shear  stress  -  slip  relationship  shown  in 
Fig.  5.2,  the  interfacial  shear  stress  prevailing  in  the  zone  of  debonding  is 
constant  and  equal  to  tf.  This  means  that  the  normal  force  distribution  in 
the  fiber  is  linear,  decreasing  at  the  rate  of  tf  per  unit  length,  where: 

tf  =  XfVjr  (5.38) 

Ir  this  study,  the  length  of  the  zone  of  debonding  will  be  noted  by  u. 
The  length  of  the  bonded  zone  is  hence  (f  -  u).  Over  this  bonded  length, 
the  same  shear  stress  distribution  prevails  as  in  the  case  where  P  <  Pcrif 
except  that  the  force  is  P'  =  (P  -  tf  u),  and  the  length  is  (f  -  u).  Therefore, 
the  shear  flow  distribution  along  the  interface  of  the  bonded  length  can  be 
expressed  as: 

tb(x)  =  (P  -  tf  u)  A.  (a*  e^*  -  B*  e-^*)  (5.39) 

where 

A*  =  1  ■  e.2M(..)  (('-  5)  *  h 

®  ^  1  .  e.2».((.u)  (■('■  q)  ='*•<'■“*  -  q]  (S  '*!) 

The  relationship  between  P  and  u  can  be  obtained  by  setting  the  value 
of  the  shear  flow  at  x  =  (f  -  u)  equal  to  tn,ax.  or: 

(P  -  tf  u)  X  (a*  e^(^  -  u)  -  B*  e-^(^  -  “))  =  t^ax  (5.42) 

or,  replacing  A*  and  B*  by  their  respective  values  given  in  Eqs.  5.40  and 
5.41,  P  can  be  found  explicitly  as: 
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tf  u  + 


"max 


1  .  q-2X{[  -  u) 


2 
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^  cMi  -  u)  + 


1  -  (l  +  -  “)) 

(5.43) 


Since  the  force  resisted  by  the  debonded  zone  can  be  easily  found  as: 

Pd  =  tf  u  (5.44) 

then.  Pi,,  which  is  equal  to  the  difference  between  the  total  force  P  (given 
in  Eq.  5.43)  and  Pd  (as  per  Eq.  5.44),  can  hence  be  found  as: 


"max 


1  .  q-2X([  -  u) 


A  -  u)  +  1^1  -  (1  +  -  u)) 


(5.45) 


To  determine  the  distribution  of  the  force  in  the  fiber  in  the  bonded 
zone,  Fi,,  the  shear  flow  distribution  in  that  zone  is  integrated: 


Pb(x)  = 


A 

i 


tb(y)  dy 


(5.46) 


A 

(P  -  tf  u)  X  J  (a""  e^y  -  B*  e'^y)  dy 


=  (P  -  tf  u)  (a*  e^*  +  B*  e-^’^  -  A*  +  B*)  (5.47) 

Naturally,  Eq.  5.47  holds  for  values  of  x  between  0  and  (f  -  u).  x  = 

0  at  the  fiber's  embedded  end,  and  x  =  (f  -  u)  at  the  point  of  demarcation 
between  the  bonded  and  debonded  zones.  Fig.  5.5  shows  the  interfacial 
bond  shear  stress  and  the  fiber  normal  force  distribution  under  conditions  • 

of  partial  debonding. 

The  force  distribution  in  the  debonded  zone  can  be  expressed  as: 


Fd(z)  = 


P  -  tf  z 


(5.48) 


x  =  0  x=^-u  \  =  [ 


5.5  -  Bond  stress  and  normal  force  distributions  under  partially 
debonded  conditions. 


Fig.  5.6  - 


Pull-out  test  under:  (a)  pre-peak  conditions,  and  (b)  dynamic 
pull-out  conditions. 
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Equation  5.48  holds  for  values  of  z  between  0  at  the  free  end,  and  u 
at  the  point  of  demarcation  between  the  bonded  and  debonded  zones.  • 


To  find  the  displacement  at  the  free  end,  Eq.  5.34  is  still  valid,  ex¬ 
cept  that  the  integral  will  be  carried  out  separately  over  each  of  the  two 
zones,  thus: 
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-P  [ 
A  F 


A  F 


J  Fd(^) 


dx  -t- 


VO 


l-\i 


dx 

/ 


(5.49) 


in  which  Fji(z)  and  Fb(x)  are  as  given  in  Eqs.  5.48  and  5.47  respectively. 
The  integral  in  Eq.  5.49  yields: 
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P(Q-l)u  - 


tfu2 

2 


(Q-2) 


^;^l-e-^(^-“)  >  Q-2 
(^■'f'’)(l+e-M^u)J  ^ 


tfUf 


^mFm 


(5.50) 


5.3.5  Dynamic  Mechanism  of  Pull-Out 

The  value  of  tne  length  of  the  debonded  zone  u  goes  from  zero  at  the 
onset  of  debonding  to  C  at  complete  debonding  (purely  frictional  bond  shear 
stresses  at  the  interface).  For  each  value  of  u,  one  value  for  the  pull-out  ^ 

load  and  the  corresponding  end  slip  exists. 

Once  u  has  reached  C,  a  dynamic  mechanism  of  pull-out  develops. 

The  relationship  between  the  pull-out  load  and  the  slip  can  now  be  derived.  • 

If  the  rigid  body  displacement  of  the  fiber  is  referred  to  as  v  (Fig.  5.6), 
then  the  pull-out  force  is: 

P  =  tf(r-v)  (5.51)  ^ 


The  force  F  within  the  fiber  at  a  distance  x  from  the  point  at  which 
the  fiber  penetrates  the  matrix  is: 


F(x)  =  tf(f-v-x) 


(5.52) 
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The  dynamic  pull-out  slip  A  is  now  equal  to  the  total  rigid  body 
•  movement  of  the  fiber  v,  added  to  the  elastic  elongation  within  the  fiber, 

or: 


A 


The  dynamic  pull-out  slip  A,  as  given  by  Eq.  5.53,  is  almost  equal  to 
the  fiber  rigid  body  movement  v.  This  suggests  an  almost  perfectly  linear 
relationship  between  the  applied  pull-out  force  P,  and  the  dynamic  pull-out 
slip  A,  of  the  form: 

P  =  tf(r-A)  (5.54) 

As  is  demonstrated  in  Chap.  6,  puP-out  tests  conducted  in  the  lab 
suggested  a  descending  branch  that  is  not  linear,  but  decays  more  rapidly 
(Fig.  5.7),  which  suggests  a  decay  in  the  the  frictional  bond.  In  other 
words,  the  assumption  that  the  bond  shear  stress  remains  constant  after  the 
bond  strength  is  reached  is  not  correct.  However,  it  is  acceptable  at 
small  values  of  slip. 

An  alternative  convenient  way  to  look  at  the  descending  branch 
would  be  to  assume  that  the  bond  shear  stress  versus  slip  curve  is  as 
shown  in  Fig.  5.8.  Therefore,  it  is  assumed  that  at  some  point  after  the 
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(5.53) 
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ig.  5.7  - 


Schematic  representation  of  a  typical  pull-out  curve. 


Mg.  5.8a  - 


Alternative  bond  shear  stress  versus  slip  relationship  with 
frictional  decay. 


Fig.  5.8b  - 


SLIP 

O  V 

Alternative  bond  shear  stress  versus  slip  relationship  with 
frictional  decay  when  smaller  scale. 
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bond  strength  X^ax  reached,  the  frictional  shear  stress  starts  decreasing 
with  the  slip.  In  what  follows,  it  will  be  assumed  that  the  slip  for  which 
the  friction  starts  decaying  is  the  slip  Aq  that  corresponds  to  the  case  of 
full  debonding,  which  is  the  onset  of  the  dynamic  mechanism.  Aq  is 
assumed  small  in  comparison  to  the  embedded  length.  Under  this 
assumption,  the  equations  given  in  this  chapter,  that  pertain  to  both  the 
fully  bonded  case  and  the  the  case  of  partial  bonding  are  still  applicable. 
At  this  stage,  a  state  of  complete  failure  of  the  elastic  bond  between  the 
fiber  and  the  matrix  exists.  Thus,  the  only  phenomenon  that  describes  the 
resistance  of  the  fiber  against  pull  out  is  friction. 


5. 3. 5.1  Normal  Contact  Pressure.  The  frictional  stress  is  equal  to  the 
normal  contact  pressure  on  the  steel  fiber  multiplied  by  a  coefficient  of 
friction;  ^ 

td  =  ^  Pn  (5.55) 

in  which  Pjvj  is  the  normal  contact  pressure  between  the  fiber  and  the 

matrix.  # 


From  examining  the  free  body  diagram  of  the  embedded  fiber  in  Fig. 
5.6,  one  can  write  from  equilibrium  that: 

dF  =  Td  V  dx  (5.56) 

where  dF  is  the  differential  of  the  local  force  in  the  fiber  at  a  distance  x. 

So  if  we  divide  Eq.  5.56  by  the  cross  sectional  area  of  the  fiber  we 

get: 

dff  =  (5-57) 

Substituting  the  right  hand  side  of  Eq.  5.55  in  Eq.  5.57  leads  to: 


dff 


2  H  Pn 
Tf 


dx 


(5.58) 


-78- 


But  from  the  shrink  fit  theory  [Timoshinko],  the  interfacial  contact 
pressure  in  a  shrink  fit  configuration  (Fig.  5.9)  with  no  load  on  the 
fiber  and  for  rf  <<rn^  is  given  by: 


'N 


-mr 


^m)  ^  -  Vf) 


-m 


(5.59) 


where  e^r  is  the  radial  shrinkage  strain  in  the  matrix  in  addition  to,  if 
present,  the  strain  due  to  any  externally  applied  confining  load. 

The  radial  shrinkage  strain  in  the  matrix  can  be  expressed  as: 


e 


mr 


5 


(5.60) 


5  in  Eq.  5.60  represents  the  matrix-fiber  misfit.  The  misfit  was 
defined  by  Pinchin  [50]  as  "the  difference  between  the  fiber  and  the  hole 
radius  in  the  absence  of  the  fiber"  (Fig.  5.9). 


5  =  rf-ro 


(5.60-a) 


5. 3. 5. 2  Poisson's  Effect.  However,  when  the  fiber  is  loaded  longitudinally 
with  a  stress  ff,  it  will  be  subjected  to  a  Poisson’s  contraction  Efr  such  that: 


This  strain  £f,  will  obviously  reduce  the  interfacial  contact  pressure 
caused  by  the  original  matrix  strain,  thus; 

£cff  ~  ^mr  ■  ^fr  (5.62) 

=  (5.63) 

substituting  the  right  hand  side  of  Eq  5.63  in  Eq.  5.59,  and  solving  the 
resulting  differential  equation,  the  stress  in  the  fiber  can  be  expressed  as 
follows: 


shrinkage 


Hole  after 
shrinkage 


(in  absence  of  the  fiber) 


Fig.  5.9  -  Fiber  -  matrix  misfit 
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(5.64) 


where  x  is  the  embedded  length  of  the  fiber. 

Multiplying  Eq.  5.64  by  the  cross  sectional  area  of  the  fiber  yields 
an  equation  for  the  pull-out  load  as  a  function  of  the  embedded  length  x; 
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(5.65) 

5.  3. 5. 3  Effect  of  Decay  in  Misfit  5.  It  is  observed  that  5,  the  fiber-matrix 
misfit,  deteriorates  (i.e.  decreases)  as  the  fiber  is  pulled  out.  This  is  due 
to  a  combined  action  of  abrasion  and  compaction  of  cement  and  sand 
particles  surrounding  the  fiber  since  the  steel  fiber  is  harder  than  the 
matrix. 

Several  attempts  were  made  in  this  study  to  model  the  decay  in 
misfit.  A  linear  and  a  parabolic  function  were  considered.  The  best  results 
were  obtained  when  the  decreasing  trend  of  the  misfit  was  assumed  ex¬ 
ponential,  as  suggested  by  comparison  with  experimental  pull-out  curves. 
The  following  expression  for  the  fiber-matrix  misfit  5  is  proposed  in  this 
study: 

5 

in  which: 
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(5.66) 


relative  slip  of  the  fiber  under  conditions  of  full 
debonding 


the  end  slip 
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6(5  =  the  initial  fiber-matrix  misfit  =  rf-ro 

4  =  a  dimensionless  constant  to  give  the  analytical 

descending  branch  the  same  asymptotic  value  as  the 
experimental  one. 

The  initial  fiber-matrix  misfit  is  calculated  by  equating  the  pull¬ 
out  load  that  corresponds  to  an  embedded  length  x  =  fin  Eq.  5.65.  That  is 
P  is  replaced  by  tf.  f.  Hence: 
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Ag  can  be  obtained  from  Eq.  5.50,  with  a  P  =  tf  f,  and  u  =  f,  or: 


(Q  -  2)  a 

2  Anf)  Eg, 


(5.68) 


Fig.  5. 10  shows  a  typical  plot  of  the  variation  of  the  misfit  as  related 
to  the  relative  slip  between  the  fiber  and  the  matrix.  This  plot  pertains  to  a 
smooth  fiber  of  diameter  0.01  inches  pulled  out  from  a  high  strength 
cementitious  matrix  containing  2%  by  volume  steel  fibers  of  the  same  kind 
as  the  one  being  pulled  out.  Additional  information  on  this  series  can  be 
looked  up  in  chapter  6,  tables  6.3  through  6.5. 

In  pull-out  tests  the  embedded  length  of  the  fiber  is  equal  to  the 
difference  between  the  original  embedded  length  and  the  additional  end 
slip,  i.e.  X  =  f  -  (  A  -  Ag  ).  On  the  other  hand,  in  pull-through  tests  the 
embedded  length  of  the  fiber  is  constant  and  equal  to  the  original  embedded 
length,  i.e.  x  =  f 


5. 3. 5.4  Frictional  Shear  Value.  The  above  described  methodology  makes 
it  possible  to  find  an  equivalent  value  for  the  frictional  shear  bond  for  any 
given  pull-out  load,  or  end  slip  A,  since: 
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P  =  V  'td(x)  X 


or 


Xd(x)  = 
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y  (r  -  A  +  Ao) 


The  value  of  P  obtained  from  Eq.  5.65  can  be  divided  by  the 
corresponding  embedded  length,  thus  giving  an  equivalent  value  for  the 
frictional  shear  bond  for  a  given  slip,  hence: 
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(5.69) 


which  is  valid  for  A  >  Ag. 

A  typical  plot  of  Xj  versus  x  where  x=/'-A+Ao  is  shown  in  Fig.  5. 1 1. 
The  initial  frictional  stress  Xf,  is  equal  to  X(i(Ao),  hence: 
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Hence,  the  decaying  frictional  stress  Xj  can  be  expressed  as: 
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5.4  PREDICTION  OF  PULL-OUT  CURVE  (PRIMAL  PROBLEM) 

For  a  given  bond  shear  stress  versus  slip  relationship  of  the  type 
discussed  above  (Fig.  5.8),  a  complete  pull-out  curve  can  be  predicted. 
The  predicted  pull-out  curve  (Fig.  5.7)  is  divided  into  three  distinct  regions 
and  can  be  obtained  using  the  following  procedure: 

1.  Pre-critical  region:  it  is  a  linear  portion;  the  critical  point  (Pcrif^crit) 
is  enough  to  describe  this  linear  zone,  which  extends  from  the  origin 
to  the  critical  point. 

2.  Partial  debonding  region:  for  each  value  of  the  length  of  the 
debonding  zone  u,  one  point  on  the  curve  can  be  found  by  computing 
the  corresponding  pull-out  force  and  end  slip  (Eq..  5.43  and  5.50). 
As  many  points  as  needed  within  this  range  can  be  found,  so  long  as 
the  value  of  u  is  bound  between  0  and  the  fiber  length  [. 

3.  The  pull-out  region:  for  each  value  of  the  end  slip  A  (Aq  <  A  <  O.or 
(0  ^  (A-Aq)  ^  0.  a  value  for  the  pull-out  load  can  be  determined 
from  Eq.5.65,  hence  leading  to  one  point  on  the  pull-out  curve. 
When  (A-Aq)  reaches  the  value  of  f,  the  force  would  be  equal  to 


zero. 


-85- 


5.5  PREDICTION  OF  BOND-SLIP  CURVE  (DUAL  PROBLEM) 

Given  an  experimental  pull-out  curve  for  a  given  fiber,  the  bond 
shear  stress  versus  slip  curve  can  be  theoretically  obtained,  assuming  that 
it  is  of  the  type  described  in  Fig.  5.8,  i.e.  linear  elastic  ascending  branch, 
followed  by  a  purely  frictional  region,  and  then  a  deteriorating  frictional 
zone(see  Eq.  5.70).  The  whole  curve  can  hence  be  described  by  five 
parameter. ,  namely  the  bond  modulus  k,  the  bond  strength  Xmax- 
constant  frictional  bond  stress  tf,  the  value  of  the  slip  at  which  the  bond 
starts  to  deteriorate  Aq,  and  the  decaying  frictional  parameters  ^  and  r\, 
describing  the  deteriorating  frictional  zone. 

The  bond  modulus  k  is  determined  from  the  slope  of  the  linear  as¬ 
cending  portion  of  the  pull-out  curve,  which  can  be  obtained  graphically. 
Furthermore,  the  value  of  Q  (Eq.  5.18)  can  be  evaluated  from  the  physical 
and  mechanical  properties  of  the  fiber  and  the  matrix.  Judgment  and 

common  sense  are  to  be  used  in  evaluating  the  area  of  the  matrix  A^- 

Indeed,  only  a  fraction  of  the  matrix  cross  sectional  area  is  effective  if  the 
proportion  of  the  specimen’s  cross  sectional  area  to  that  of  the  fiber  is 

relatively  large.  However,  it  can  be  shown  that  the  solution  to  the  dual 

problem  is  insensitive  to  the  value  of  A^-  Once  Q  and  the  slope  j  are 

known,  X  can  be  solved  for  in  Eq.  5.36  by  iteration  or  by  some  numerical 
procedure.  Having  found  X,  and  using  Eq.  5.20,  K  can  now  be  solved  for: 

K  =  Q  (5.71) 

Finally,  the  value  of  the  bond  modulus  k  can  be  derived  from  Eq. 

5. 17: 

■c  =  -  A„  E„  K  (5.72) 

¥ 

The  next  step  would  be  to  evaluate  the  bond  strength  of  the  interface 
'^max  the  frictional  bond  if.  The  peak  load  Pp  and  the 

corresponding  end  slip  Ap  are  to  be  studied  first.  The  peak  load  will 
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generally  occur  under  partial  debonding  conditions.  The  value  of  u 
corresponding  to  that  case  is  the  value  of  u  that  would  maximize  P.  Hence, 
for  u  =  Up,  we  have: 


where  P  is  as  given  in  Eq.  5.43: 


P 


tf  u  + 


X 


.  -  u) 

[^1  -  (l  +  -  u)) 


To  make  the  differentiation  easier,  a  change  of  variables  is  in  order. 
Introducing  the  variable  X  defined  as: 


X 

tMc  ■  u) 

(5.73) 

then 

dX 

du 

X  e->-(^  -  iJ) 

(5.74) 

and 

dP  dX 

du 

dX  du 

thus 

(^].X 

=  0 

=  Xp 

(5.75) 

where 

P(X)  = 


^ln(X)  +  tf  ft- 


(5.76) 


1  -  X2 


1  -  ^ 


3)  *  X') 


'•max 


X  =  Xp  XX, 


(5.77) 


-2  2 

Q 


^  -  qPp 


Q 


'  ■  5)  ’‘z  h  ’‘p  ^ 


(5.78) 
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also  P(Xp) 


^ln(Xp)  +  tf  ^ 

^  Am 


^  X  + 

Q  Ap  + 


1  - 


1  +  x“^ 

pj 


(5.79) 


and  A(Xp) 


p(Q  ■  +!]  -  ^  +ff 


(Q 


A  F 


P-  tf 


/|n(2y_  ^  (Q  -_2) 


M. 


1  +  X  I  ■  'f  ^ 


(5.80) 


To  get  the  values  of  tf  and  tn,ax>  ^  system  of  three  non-linear 
equations  in  three  unknown  has  to  be  solved,  the  three  unknowns  being  tf, 
tmax’  21*^^  ^p-  f'^st  equation  is  obtained  by  setting  the  right-hand  side 
of  Eq.  5.78  equal  to  zero.  The  second  equation  is  obtained  by  setting  the 
right-hand  side  of  Eq.  5.79  equal  to  Pp,  the  maximum  recorded  load  on  the 
(P-A)  curve,  and  the  third  equation  consists  of  equating  the  right-hand  side 
of  Eq.  5.80  to  Ap,  the  end  slip  corresponding  to  Pp. 


The  corresponding  system  of  three  equations  in  three  unknowns  is 
thus  as  follows: 


(5.81) 


^ln(Xp)  -h  tf  f-K  ^ 


Pp  (5.82) 


# 


1 


I*W«] 


Pp  (Q  -  I)  ^ 


X  )  2  I  V-  -^'(Q  -  2)  -  tf  r 


\ 


Pm 


Pp*  tf 


/In(Xo) 


+/■ 


(Q  -  2)  fi  -  X, 
I  +  X„ 


'^m  Pm 


=  A, 


(5.83) 


the  value  of  Up  corresponding  to  the  value  found  for  Xp  has  to  be  between 
0  and  C.  The  above  three  equations,  when  solved,  lead  to  the  values  of  tf, 
^max  equivalently  Up.  The  numerical  solution  is  described  in 

Appendix  V-A. 


The  experimental  pull-out  curves  developed  in  the  course  of  this 
investigation  suggested  a  steep  initial  decay  in  the  post-peak  pull-out 
behavior,  that  is  when  A  is  close  to  Aq.  A  paiametric  study  of  Eq.  5.70 
indicated  that  the  factor  reflecting  the  rate  of  initial  decay  is  r|.  To  best 
simulate  the  real  experimental  curves,  a  good  value  for  was  found  to  be 
near  0. 2. 


To  evaluate  the  factor  one  point  on  the  experimental  pull-out  curve 
is  needed,  say  (P^,  A^).  Then,  from  Eq.  5.65,  and  using  the  values  of  Px 
and  Ax,  5  can  be  solved  for: 

A\1 

JJ 

(5.84) 

Tl’f  .  substituting  in  Eq.  5.66,  ^  can  be  solved  for: 


Px  Vj 
Ef  Jt  Tf 


1  -  Exp 


f  -2  Vf  li  (f 


V 


Aq) 


E,  .  (1 


Vf) 


V 


e-(0»  *  .  f.^.) 

1  c-(f-Ax-Ao)0-2 

4 


C-(Ax-Ao) 


(5.85) 


where  5  is  as  given  in  Eq.  5.84,  and  5^  is  as  given  in  Eq.  5.67. 

With  the  five  basic  parameters  tc,  Tf  ,  Aq,  and  ^  known,  the 

whole  bond  shear  stress  versus  slip  relationship  can  be  constructed. 
Examples  and  correlation  with  experimental  observations  are  presented  in 
Figs.  5.12-5.18.  Those  examples  show  the  application  of  the  analytical 
model  in  solving  the  dual  problem  as  explained  earlier  in  section  5.5  of  this 
chapter.  The  numerical  data  used  in  the  solution  is  presented  in  Appendix 
V-B.  Examples  on  the  application  of  the  proposed  model  in  solving  the 
primal  problem  i.e.  predicting  the  pull-out  curve  are  presented  in  Chapter  6 
(Figs.  6.22-6.27). 
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Fiber  Volume  Fraction 
Additives  1  None 
Fiber  Diameter  =  0.019’ 
Embedded  Length  =  1’ 


SLIP  (in.) 

Bond  shear  stress  versus  slip.  HISN  series,  full  scale. 
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Bond  shear  stress  versus  slip.  HISN  series,  small  scale. 
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5.13a 


SLIP  (in,) 

Bond  shear  stress  versus  slip,  HOSN  series,  full  scale. 
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5. 13b  -  Bond  shear  stress  versus  slip,  HOSN  series,  small  scale 


.0004 

SLIP 


.0006 

(in.) 


,0010 


SHEAR  STRESS  (psi)  ^  BDHD  SHEAR  STRESS 


Matrix  Strength  i  Lou 
Pull-out  Fiber  >  Smooth 
Fiber  Volume  Fraction  =  0’/. 
Additives  i  None 
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Embedded  Length  =  1* 
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Bond  shear  stress  versus  slip,  LOSN  scries,  full  scale 
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Bond  shear  stress  versus  slip,  LOSN  series,  small  scale. 
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5. 15b.  Bond  shear  stress  versus  slip,  H2SL  series,  small  scale. 
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Bond  shear  stress  versus  slip,  AISN  series,  full  scale 
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5.16b  -  Bond  shear  stress  versus  slip,  AISN  series,  small  scale. 
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APPENDIX  V-A 

APPLICATION  OF  NEWTON’S  METHOD  FOR  SOLVING  A  SYSTEM  OF 
NON-LINEAR  EQUATIONS  TO  THE  DUAL  PROBLEM  OF  PULL-OUT 


A.  1  GENERAL 


A  system  of  n  non-linear  equations  in  n  unknows  can  be  expressed  as 
follows: 


F(C)  =  0 


(A.  1) 


where  F  is  a  vector  of  functions  Fj,  F2,  •••Fn,  evaluated  at  the  solution 
vector  of  the  unknows  zj,  Z2,  ..  Zn,  or: 


r.;i 


and 


(  Zl  Z2  .  .  .  Zn  ) 


hence  Eq.  (A.  1)  is  equivalent  to: 
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^l(zi,  Z2.  -Zn)^ 

F2(zi,  Z2,  ...Zn) 

0 

Koj 

Newton  method,  discussed  in  the  following  secti  )n,  converges  to  the 
exact  solution  vector  provided  the  initial  guess  vector,  say  is  close 
enough  to  and  provided  the  Jacobian  F'  of  F  is  continuous  and  the 
matrix  F’(^)  is  invertible. 


A.  2  ALGORITHM 

The  algoritm  to  be  followed  in  Newton’s  method  is  to  start  from  the 
initial  guess  vector  and  successively  improve  it  till  the  improvement 
becomes  insignificant.  Given  a  guess  z^*"),  an  improved  solution 
over  z^"™)  can  be  obtained  by  the  following  equation: 

z(fn+l)  =  [F’(z('"))I->  FCz^f"))  (A. 3) 


or 


F'(z(n>))  (Az('"+l))=  F(z(>n)) 


(A.4) 


where 


^z(n™+l)  =  z("*+l)  -  zl*") 


(A.5) 


F’  is  the  Jacobian  of  the  vector  of  functions  Fj,  F2,  •••Fn- 
Therefore,  F'  is  an  nxn  matrix,  defined  by  its  F'jj  element,  that  is  the 
element  in  the  i^**  row  and  the  column,  as: 


(A.6) 
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Therefore,  the  way  to  solve  the  system  of  n  non-linear  equations  in  n 
unknowns  consists  of  defining  the  vector  F  and  the  corresponding  Jacobian 
matrix  F’ ,  and  then,  starting  with  an  initial  guess  vector  successively 
improving  the  solution  by  solving  Eq.  A. 4  for  and  then  finding  an 

improved  solution  in  the  following  equation: 

z(m+l)  =  ^z(m+l)  +  2('^)  (A. 7) 

A.  3  Application  of  Newton's  Method  to  the  Dual  Problem  of  Pull-out. 

It  can  be  recalled  from  Chap.  7  that  the  dual  pull-out  problem 
consists  of  solving  three  non-linear  equations  in  tf,  t^ax*  "^^e 

equations  in  question  are  rewritten  here: 

tf  e^(^  - 
- + 

X 


UELax. 


V 


^  Q-lXa  -  U)  .  4  A  e-X(r  -  U)  .  ^  ^ 


1  -  e*2^(f  -  u)+  ^  -  u)  +  ^1  .  XjJ 


=  0 


J 


(A.8) 


tf  u  + 


1  -  Q-2MC  -  u) 


X  2 


±  cMl  -  u)  +  J"!  .  (i  +  c-2X{[  -  u)) 


=  P, 


(A.9) 


Pp(Q-l)u  -  -V  (Q-2)  *  (Pp-.tu)(;4.,(t.„))  Y 


Am^m 


=  A, 


(A.  10) 


The  vector  of  functions  F  is  equal  to: 
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F2 

\^lj 

The  vector  of  unknowns  ^  is: 

r  u  \ 


c 


^max 

V  y 


The  functions  Fj,  F2,  and  F3  are  defined  respectively  as: 


^  1  ^max* 

^  -2 


X(/  -  u) 


^ma2L 


^  c-2X(C  -  u)  .  4  ^1.  e-^(^  ■  ^  ^ 


V 


1  .  e-2Mr  -  U)^  2  .  u)  ^  1^1  . 


J 


(A.  11) 


(A. 12) 


^max> 


tf  u  + 


^uiaa. 


1  .  q-2X(C  -  u) 


(A. 13) 


^  cMc  -  u)  + 


(1  ■  i)  (1  ^  e-2Mt  ■  u)) 


-  P, 


(A.  14) 


Pp(Q-l)a  -  ^  (Q-2)  +  ^  -  t,uf' 


(A.  15) 


The  Jacobian  matrix  F'  can  be  written  as  follows: 


-100- 


The  9  elements  of  the  Jacobian  matrix  F’  are  given  next: 


F' 


1,1 


aF, 


-  u) 

X 


(A.  16) 


(A.  17) 


(A.  20) 
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F' 


2.2  = 


aF2 
9t 


max 


]_  _ 1  -  -  u) _ 

^  1  qMC  -  u)  +  a  -  (l  +  e-2^(^  -  u)) 


F' 


3.2 


9F' 

Tt 


max 


(A.21) 


(A. 22) 


Since  the  partial  derivatives  involving  the  variable  u  are  rather 
complex  and  involved,  the  functions  Fi  and  F2  will  be  rewritten  as  follows 
first: 


tf  ^  tmax  4>(^) 

X  X  Tl(u) 


(A.  23) 


in  which 

<l)(u)  =  •  u)  .  4  J'l  .  ^je-Xd-u).  1  (A. 24) 

and 

n(u)  =  to  -  ■  “)+  ^  c-UI  ■  u)  +  ^1  .  ijJ(A.25) 

F2(>f.<max.»)  =  ‘f  “  ^  ^  '  Fp 

X  GJ(u) 


in  which: 

H(u) 

1  .  c-2X(r  u) 

(A.27) 

and 

0(u)  = 

^  cM(  -  u)  +  j"! 

-  ^)(l  +  *  “0 

(A.28) 
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F' 


1.3 


in  which 

4)’  (u) 


and 


n'(u) 


F' 


2.3 


in  which: 

^'(u) 


and 


GJ’(u) 


I  ^XiC  -  u)  +  4)'(u)  Ti(u)  -  0(u)  Ti'(u) 

X  ri2(u) 


u)  .  4  X  J'l 


Q-\([  -  U) 


2  X  (l  -  e-2^(^  -  u)  +  Q-XiC  -  u) 


1  -  e'2^(^  -  u)+  A  Q-X{[  -  u)  +  ("i  .  ± 


3u 

j  Ifnax  (u){5(u)  -  |J.(u)tO  (u) 

^  X  G52(u) 


-  2  X  c-'2-M(  -  u) 


2  X 
Q 


qMC  -  u)  +  2  X 


U) 


I 


► 


(A.29) 


(A. 29) 


* 

J 

)) 

(A.  30) 


(A. 31) 


(A.32) 


(A.33) 
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Pp  (Q  -  1)  -  2  tf  u  (Q  -  2)  -  tf 

_  2  Pn  e-^(^  -  (Q  -  2) 

(1  +  ■  u))2  ■ 

2  tf  u  e-^(^  ~  (Q  -  2) 

(1  +  e-^(^  -  u))2 


-  2  1  -  e-^(^  •  u) 
X  1  +  e-^(^  -  u) 


(A.34) 
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APPENDIX  V-B 

NUMERICAL  DATA  USED  IN  THE  SOLUTION  OF  THE  DUAL 

PROBLEM 
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TEST  RESULTS: 

.  Slope  (P/A)  :  44x10^  Ib/in 

.  Max.  pullout  load:  15.4  lb 
.  Slip  at  max  load  :  8.8x10  '*  in 
.  Point  on  the  descending  branch 


HISN  series 

.  High  strength  concrete  (8.6  ksi). 
.  1%  fiber  by  volume. 

.  No  additives. 

(4.3  lb,  0. 15  in) 


MEASURED  PARAMETERS: 

.  Fiber  diameter  :  0.019  in. 

.  Fiber  Embedded  length  :  1  in. 


ASSUMED  PARAMETERS: 

.  Fiber  modulus  of  elasticity  (Ef)  :  3x10^  psi. 

.  Matrix  modulus  of  elasticity  (E^)  :  3x10^  psi. 

.  Effective  matrix  area  surrounding  the  fiber  :  1  in^. 

.  Coooefficient  of  friction  (stell-mortar)  :  0.5. 

.  Fiber  poisson's  ratio  :  0.2. 

.  Matrix  poisson's  ratio  :  0.2. 

COMPUTER  PROGRAM: 


Solving  the  three  basic  nonlinear  equations  yielded  the  following: 


4  =0.786 

K  =  56.9xl0Mb/in3 

Xm  =  258  psi 

Xf  =  258  psi 
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TEST  RESULTS: 

.Slope  (P/A)  :  25x103  Ib/in 

.  Max.  pullout  load;  15.6  lb 
.  Slip  at  max  load  :  9.25x10  ''  in 
.  Point  on  the  descending  branch 


HOSN  series 

.  High  strength  concrete  (8.6  ksi). 
.  0%  fiber  by  volume. 

.  No  additives. 

(1.3  lb.  0. 1  in) 


MEASURED  PARAMETERS: 

.  Fiber  diameter  :  0.019  in. 

.  Fiber  Embedded  length  :  1  in. 


ASSUMED  PARAMETERS: 

.  Fiber  modulus  of  elasticity  (Ef)  ;  3x10^  psi. 

.  Matrix  modulus  of  elasticity  (E^,)  :  3x10^  psi. 

.  Effective  matrix  area  surrounding  the  fiber  :  1  in^. 

.  Coooefficient  of  friction  (stell-mortar)  :  0.5. 

.  Fiber  poisson's  ratio  ;  0.2. 

.  Matrix  poisson's  ratio  :  0.2. 


COMPUTER  PROGRAM: 

Solving  the  three  basic  nonlinear  equations  yielded  the  following: 


^  =  1-3 

ic  =  0.841x10®  Ib/in3 

tm  =  334  psi 

tf  =  257  psi 


SLIP  (in.) 

Bond  iheir  itrcit  versus  slip,  HOSN  series,  sniill  scile. 
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Slope  (P/A)  :  52.1x103  Ib/in  .  Low  strength  concrete  (4.8  ksi). 


Max.  pullout  load:  12.. 6  lb 
Slip  at  max  load  :  7.3x10  '*  in 


.  0%  fiber  by  volume. 
.  No  additives. 


Point  on  the  descending  branch  (2.97  lb,  0. 15  in) 


.  Fiber  diameter 


:  0.019  in. 


.  Fiber  Embedded  length  :  1  in. 


.  Fiber  modulus  of  elasticity  (Ef) 

.  Matrix  modulus  of  elasticity  (E^) 


:  3x10^  psi. 
:  3x10^  psi. 


.  Effective  matrix  area  surrounding  the  fiber  :  1  in^. 
.  Coooefficient  of  friction  (stell-mortar)  :  0.5. 


.  Fiber  poisson's  ratio 
.  Matrix  poisson's  ratio 


:  0.2. 
:  0.2. 


Solving  the  three  basic  nonlinear  equations  yielded  the  following: 
E  =0.91  SI - 


=  5.25x10^  lb/in3 


=  211  psi 
=  211  psi 


Motrix  Strength  >  Low 
PuU-out  Tibfr  »  Smooth 
Tibfr  Volume  Tractioo  ■  OX 
Addttlvtf  •  Mono 
Tibor  Dtomotor  ■  0019* 
Cnboddod  Longth  ■  T 


—  Actual  Dobon^ng 

Asconding  Brooch 


.0002  .0004 


SLIP  (In.) 

Bond  ihcir  itresi  vertui  slip.  LOSN  leriei,  imill  icile. 
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Slope  (P/A) 


67.8x103  Ib/in  .  High  strength  concrete  (8.6  ksi). 


Max.  pullout  load:  22. 1  lb 


2%  fiber  by  volume. 


Slip  at  max  load  :  15.9x10'^  in  .  Latex  as  an  additive. 
Point  on  the  descending  branch  (4.8  lb,  0. 1  in) 


Fiber  diameter 


0.019  in. 


Fiber  Embedded  length  :  1  in. 


.  Fiber  modulus  of  elasticity  (Ef) 

.  Matrix  modulus  of  elasticity  (E^) 


:  3x10^  psi. 
:  3x10^  psi. 


.  Effective  matrix  area  surrounding  the  fiber  :  1  in^. 
.  Coooefficient  of  friction  (stell-mortar)  :  0.5. 

.  Fiber  poisson's  ratio  :  0.2. 

.  Matrix  poisson's  ratio  :  0.2. 


Solving  the  three  basic  nonlinear  equations  yielded  the  following; 
t  =  0.938  h - - - — 


=  8.79xl0<*  lb/in3 
=  1390  psi 


=  259  psi 


Hotrix  str»ngth  .  High 
Pvjll-out  Tibvr  '  Snooth 
Tibfr  Voluht  frbcttorv  ■ 
Additivfs  '  Lat»« 

Tibtr  Dianttrr  ■  0019’ 
Enb»ad»d  L»n.;th  •  !• 

Actual  Dtboottng 

Asctndtno  Branch 


SLIP  (in.) 

Bond  iheir  iireii  veriui  ilip,  H2SL  series,  smill  scile. 
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Slope  (P/A)  :  55x10^  Ib/in 

Max.  pullout  load:  10.8  lb 


.  Average  strength  concrete  (7.4  ksi). 
.  1%  fiber  by  volume. 


Slip  at  max  load  :  6.34x10'^  in  .  No  additives. 


Point  on  the  descending  branch 


(5  lb.  0. 15  in) 


Fiber  diameter 


:  0.019  in. 


Fiber  Embedded  length  :  1  in. 


.  Fiber  modulus  of  elasticity  (Ef) 

.  Matrix  modulus  of  elasticity  (E^) 


:  3x10^  psi. 
:  3x10®  psi. 


.  Effective  matrix  area  surrounding  the  fiber  :  1  in^. 


.  Coooefficient  of  friction  (stell-mortar) 


Fiber  poisson's  ratio 
Matrix  poisson's  ratio 


Solving  the  three  basic  nonlinear  equations  yielded  the  following: 
£  =  0.0337 


=  5x10®  Ib/in^ 
=  200  psi 
=  180  psi 


Matriw  Strength  •  Avfragt 
PuU-Owt  rib*r  •  Smooth 
FJbrr  Volu«f  TroCtJO^'  •  1^ 
Additivvs  •  Nonf 
Tibfr  Diandttr  ■  0.019* 
Cnbrddtd  Length  •  l* 

Actw«(  Ourv»g  D#bond»9 
AfCtnding  Branch 


SLIP  On) 

Bond  their  stresi  versus  slip,  AISN  series,  tmill  tcile. 
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TEST  RESULTS: 

.  Slope  (P/A)  :  45x10^  Ib/in 

.  Max.  pullout  load:  12.6  lb 
.  Slip  at  max  load  :  7. 18x10-^  in 
.  Point  on  the  descending  branch 


AQSN  series 

.  Average  strength  concrete  (7.4  ksi). 
.  0%  fiber  by  volume. 

.  No  additives. 

(4. 9  lb,  0. 1  in) 


MEASURED  PARAMETERS: 

.  Fiber  diameter  :  0.019  in. 

.  Fiber  Embedded  length  :  1  in. 


ASSUMED  PARAMETERS: 

.  Fiber  modulus  of  elasticity  (Ef)  :  3x10^  psi. 

.  Matrix  modulus  of  elasticity  (E^)  :  3x10^  psi. 

.  Effective  matrix  area  surrounding  the  fiber  :  1  in^. 

.  Coooefficient  of  friction  (stell-mortar)  :  0.5. 

.  Fiber  poisson' s  ratio  :  0.2. 

.  Matrix  poisson' s  ratio  :  0.2. 


COMPUTER  PROGRAM: 

Solving  the  three  basic  nonlinear  equations  yielded  the  following: 


4  =*1 

ic  =  3.87xl0Mb/in3 

Xm  =210  psi 

Xf  =  210  psi 


a 

ui 
a: 

h-  — 

QC 

<  § 
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X 
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r 

a 
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Putl-out  rib#r  »  Smooth 
T'bfr  VoluH#  Traction  •  0^ 
Additives  •  Non# 
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_ -  Actual  OjrmQ  Otteon^O 
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CHAPTER  VI 


EXPERIMENTAL  INVESTIGATION  AND  RESULTS 


6.1  EXPERIMENTAL  PROGRAM 

Table  6. 1  summerizes  the  experimental  program  carried  out  in  this 
study.  A  total  of  fifty-two  series,  four  fibers  each,  were  tested.  The  pro¬ 
gram  was  designed  in  such  a  way  so  as  to  investigate  the  effects  of  the 
fiber  type,  matrix  strength,  fiber  reinforcement  ratio,  fiber  diameter,  fiber 
embedment  length  ,and  additives  on  the  interfacial  bond  in  cementitious 
composites.  Test  series  were  given  names  that  would  indicate  the  mix  and 
fiber  type  as  well.  A  typical  code  name  would  consist  of  four  alphanumeric 
characters.  The  first  character  is  a  letter  that  indicates  the  matrix  strength: 
H  for  high  strength,  A  for  average  strength,  and  L  for  low  strength.  The 
second  character  is  a  numeral  that  indicates  the  fiber  volume  fraction  in  the 
matrix,  such  as  0,  1,  2,  or  3.  The  third  character  is  a  letter  that  indicates 
the  fiber  type;  S  for  smooth,  D  for  deformed,  and  H  for  hooked.  Finally, 
the  fourth  and  last  character  indicates  the  additive  used:  F  for  fly  ash,  L  for 
latex,  M  for  microsilica,  and  N  for  no  additives.  A  typical  code  name  for  a 
specimen  could  be:  LOSN  for  a  Low  strength  matrix  with  0%  fiber 
reinforcement,  in  which  the  pull-out  fiber  is  Smooth,  and  where  No 
additives  are  used.  When  a  different  diameter  or  embedded  length  is  used  a 
remark  is  added  near  the  specimen  code. 
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Table  6. 1  -  Experimental  Program 


Vf 

Fiber  Type 

Additives 

Code 

0 

Smooth 

None 

HOSN 

0 

Smooth 

None 

HOSN,  L=0.5" 

0 

Smooth 

Latex 

HOSL 

0 

Hooked 

None 

HOHN 

0 

Hooked 

None 

HOHN  L=0.5" 

0 

Deformed 

None 

HODN 

1 

Smooth 

None 

HISN 

1 

Deformed 

None 

HIDN 

2 

Smooth 

None 

H2SN 

2 

Smooth 

Latex 

H2SL 

2 

Smooth 

Microsilica 

H2SM 

0 

Smooth 

Microsilica 

HOSM 

2 

Smooth 

H2SF 

2 

Deformed 

None 

H2DN 

3 

Smooth 

None 

H3SN 

3 

Deformed 

None 

H3DN 

0 

Smooth 

None 

AOSN 

0 

Smooth 

Latex 

AOSL 

0 

Smooth 

AOSF 

0 

Smooth 

None 

0 

Smooth 

None 

BiHMBia 

0 

Smooth 

None 

0 

Smooth 

None 

0 

Smooth 

None 

wagitiiB 

0 

Hooked 

None 

AOHN 

mm 

Hooked 

None 

AlHN 

0 

Hooked 

None 

AOHN  L=0.5" 

0 

Deformed 

None 

AODN 

0 

Deformed 

None 

A0DN,L=0.5" 

Average 

1 

Smooth 

None 

AISN 

Average 

2 

Smooth 

None 

A2SN 

Average 

3 

Smooth 

None 

A3SN 

Average 

1 

Deformed 

None 

AIDN 

Average 

1 

Deformed 

None  A1DN,L=0.5" 

Average 

2 

Deformed 

None 

A2DN 

Low 

0 

Smooth 

None 

LOSN 

Low 

0 

Hooked 

None 

LOHN 

Low 

1 

Hooked 

None 

LIHN 

Low 

3 

Hooked 

None 

L3HN 

Low 

0 

Deformed 

None 

LODN 

Low 

0 

Deformed 

Latex 

LODL 

Low 

0 

Deformed 

Microjilica 

LODM 

Low 

1 

Smooth 

None 

LISN 

Low 

2 

Smooth 

None 

L2SN 

Low 

2 

Deformed 

None 

L2DN 

Low 

3 

Deformed 

None 

L3DN 

SIFCON 

12 

Smooth 

None  L 

=  r’.D=o.or' 

SIFCON 

12 

Smooth 

None  L 

=  1.D=0.019" 

SIFCON 

12 

Smooth 

None  L 

=  1",D=0.03'’ 

SIFCON 

12 

Smooth 

None  L 

=  1'\D=0.04'’ 

0 

Deformed 

None 

AODN  * 

0 

Hooked 

None 

HOHN  * 

*  [Pull-out  Fiber  was  greased  on  the  surface]. 

6.2  TEST  SET-UP 

Each  test  scries  consists  of  a  block  of  matrix  7"x4"xr'  crossed  by 
four  pull-out  or  pull-through  fibers  (Fig.  6.1).  Hence  a  series  represents 
four  tests.  The  test  set-up  used  in  this  study  is  illustrated  in  Fig.  6.2.  As 
can  be  seen  in  the  figure,  the  top  end  of  the  fiber  is  held  by  a  specially 
designed  grip  attached  to  the  load  cell  of  an  INSTRON  machine,  while  the 
bottom  end  of  the  same  fiber  is  attached  to  a  miniature  LVDT,  connected  to 
the  bench  of  the  machine.  Also,  the  displacement  of  the  top  end  of  the 
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test  SET-UP 
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fiber  is  monitored  by  means  of  another  miniature  LVDT  installed  inside  the 
grip.  The  end  of  the  core  of  this  LVDT  touches  the  top  surface  of  the 
specimen.  The  upward  movement  of  the  crosshead  of  the  testing  machine 
applies  a  pull-out  force  on  the  top  end  of  the  fiber.  The  value  of  the  force, 
as  well  as  the  movements  at  both  ends  of  the  fiber  are  recorded  by  a  Data 
Acquisition  System,  and  stored  on  a  disk.  Although  only  the  readings  of 
the  top  LVDT  are  needed  to  obtain  the  ascending  branch  of  the  full  pull-out 
curve,  and  hence  the  bond  shear  stress-slip  relationship,  the  reading  of  the 
bottom  LVDT  is  useful  in  determining  the  load  at  which  full  debonding 
occurs.  The  LVDTs  used  are  Schaevitz,  250  MHR,  with  a  stroke  of  ±0.25 
in,  and  a  sensitivity  of  IxlO''*  in. 

When  the  smooth  or  deformed  fibers  were  used  (0.3  and  0.5  mm  in 
diameter  respectively),  a  100  lb.  load  cell  was  used.  The  hooked  fibers  ne¬ 
cessitated  the  use  of  a  1,000  lb.  load  cell. 

The  tests  were  all  crosshead  movement-controlled.  The  speed  of  the 
crosshead  movement,  and  hence  of  the  test  was  varied  during  the  test.  A 
typical  test  started  off  rather  slow  (0.002  in. /min)  for  the  ascending 
branch,  and  ended  at  a  much  higher  rate  (  0.1  in/min  )  for  the  steady 
portion  of  the  descending  branch  . 

6.3  VARIABLES  INVESTIGATED 
6.3.1  Fiber  Parameters 

6.3.  1.  1  Pull-out  Fiber  Type.  Three  different  types  of  fibers  were  pulled 
out,  namely:  smooth,  deformed,  and  hooked.  While  the  actual  length  of  the 
fibers  may  differ,  the  embedment  length  was  kept  equal  to  1"  for  most  of 
the  specimens.  As  seen  in  Table  6.1,  two  other  embedment  lengths  (0.5" 
and  1.5")  were  selectively  used.  The  information  about  most  of  the  pull-out 
fibers  used  is  summarized  in  Table  6.2.  Fig.  6.3  shows  the  shapes  of 
these  fibers  . 


Table  6.2  -  Fiber  Specifications. 


Fiber  Type 

Code 

Diameter  (in. ) 

S  mooth 

S 

0.019 

Deformed 

D 

0.018 

Hooked 

H 

0.03 

6.3. 1.2  Fiber  Diameter.  For  smooth  fibers,  four  different  fiber  diameters 

were  used  ,  0.01  ",  0.019"  [control  diameter  ]  ,  0.03  ",  and  0.04  "  .  In 

most  cases  the  embedment  length  was  1"  ,  and  the  matrix  was  plain  of 
average  strength  . 

6.3. 1.3  Fiber  Embedment  Length.  Three  different  embedment  lengths  were 
investigated  ,  namely  ,  0.5  ",  !"[  control  length  ]  ,  and  1.5".  For  the  0.5" 
and  1.5"  cases,  the  matrix  was  plain  mortar  of  average  strength  and  the 
diameter  was  0.019  "  for  smooth  fibers,  0.018"  for  deformed  fibers, and 
0.03"  for  hooked  fibers. 

6.3.2  Matrix  Parameters 

6.3.2. 1  Matrix  Strength.  Three  different  mix  proportions  were  used  to 
create  low  strength,  average  strength,  and  high  strength  matrices.  The 
sand-cement  ratio  was  kept  constant  at  2.0,  while  the  water-cement  ratio 
used  was  0.40,  0.50,  and  0.60  for  the  high  strength,  average  strength,  and 
low  strength  mixes  respectively.  Type  I  cement  and  silica  sand  -  ASTM 
C109  were  used.  For  the  SIFCON  matrix,  a  very  fine  silica  sand  (type  400) 
was  used.  The  mix  proportions  are  summarized  in  Table  6.3. 


Table  6.3  -  Mix  Proportions 


Matrix  Strength 

Code 

w/c  ratio 

Cement 

Sand 

High 

H 

0.4 

1 

2 

8650 

mm 

0.5 

1 

2 

7400 

Low 

n 

0.6 

1 

2 

4850 

SIFCON 

- 

0.52 

1 

1 

6000 

•  from  average  tests  on  3"x6"  cylinders. 


6.  3.2.2  Additives.  Three  different  additives  were  used  in  some  mixes  to 
see  whether  or  not  their  use  would  affect  the  bond  properties  of  the  inter¬ 
face.  The  additives  used  were  Latex,  Microsilica,  and  Fly  Ash.  No  two 
additives  were  used  simultaneously  in  the  same  mix.  Information  on  the 
proportion  of  each  additive  is  summarized  in  Table  6.4. 


Table  6.4  -  Additives. 


Additive 

Symbol 

%  of  cement  weight 

%  of  solids  in  Additive 

Latex 

L 

20 

50 

F 

20 

100 

Microsilica 

M 

20 

50 

6.  3. 2. 3  Fiber  volume  fraction.  Four  different  volume  fractions  of  fibers 
were  used  in  the  specimens:  0%  fiber  volume  fraction  (control 
unreinforced  mortar),  1%,  2%,  and  3%. 

In  addition,  four  Slurry  Infiltrated  Fiber  Concrete  (SIFCON) 
specimens  were  tested  to  study  the  effect  of  fiber  entanglement.  The  fiber 
volume  fraction  in  SIFCON  (Vf  =  11%)  is  usually  significantly  higher  than 
that  in  normal  fiber-reinforced  concrete  (up  to  3%).  In  all  cases  of 
reinforced  matrices,  the  same  type  of  fibers  (hooked,  30/50)  were  used  in 
the  mix  . 


6.4  TESTING  PROCEDURE 

Typical  pull-out  tests  were  performed  on  52  series  of  four  pull-out 
fibers  each.  As  was  mentioned  in  §6.2,  different  combinations  of  fibers, 
matrices,  and  additives  were  used.  During  these  pull-out  tests,  the  pull-out 
force,  the  end  slip,  as  well  as  the  movement  of  the  bottom  end  (free  end)  of 
the  pull-out  fiber  [for  pull-through  tests  only]  were  monitored.  The  pull¬ 
out  force  and  the  displacement  at  the  top  end  of  the  fiber  were  used  to 
develop  pull-out  curves,  i.e.  pull-out  force  versus  end  slip  curves. 
Typically,  for  each  series,  four  different  pull-out  curves  were  obtained. 
For  smooth  fibers  pull-out  and  pull-through  tests  were  conducted  at  a 
crosshead  speed  of  0.002  in/min  up  to  the  peak  load.  For  deformed  and 
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hooked  fibers  the  crosshead  speed  ranged  from  0.002  to  0.02  in/min  up  to 
the  peak  load.  The  speed  was  then  increased  up  to  0. 1  in/min  in  the  steady 
portion  of  the  descending  branch. 

For  any  given  pull-out  test,  both  the  pull-out  load  P  (lbs.)  and  the 
end  slip  A  (in.)  at  the  section  where  the  fiber  enters  the  matrix  were 
recorded  by  a  data  acquisition  system  and  stored  in  a  data  file  on  IBM- 
Diskettes.  The  data  were  processed  as  described  in  Fig.  6.4  .  The  pull-out 
load  versus  end  slip  relationship  was  then  developed  and  plotted  using  such 
data  files.  Plots  with  different  vertical  and  horizontal  scales  were  used 
depending  on  the  specimen  type. 

6.5  CURVE  AVERAGING  PROCEDURE 

One  representative  curve  was  needed  for  each  of  the  series  tested. 
The  criterion  used  to  chose  the  proper  averaging  method  was  that  the  aver¬ 
age  curve  has  a  peak  pull-out  load  Pp  equal  to  the  average  of  the  peak  pull¬ 
out  loads  recorded  for  the  individual  pull-out  tests.  Furthermore,  the 
average  curve  was  to  have  a  corresponding  end  slip  at  peak  Ap  equal  to  the 
average  end  slips  at  peak  for  the  individual  tests.  Finally,  equal  weights 
were  to  be  assigned  to  all  curves  used  in  averaging. 

The  average  curves  were  obtained  using  a  computer  program  written 
by  D.  E.  Otter  [47],  The  method  used  to  obtain  these  curves  was  as  fol¬ 
lows:  the  peak  pull-out  loads  for  all  tests  in  the  series  were  first  determined 
and  averaged,  and  so  were  their  corresponding  end  slips.  The  average  peak 
pull-out  load,  as  well  as  its  corresponding  end  slip  were  hence  determined. 
The  pull-out  loads  on  the  ascending  (pre-peak)  branch  of  each  pull-out 
curve  were  evaluated  at  50  equal  end  slip  intervals.  The  interval  was 
selected  in  terms  in  each  curve's  own  peak  displacement  .  The  average  of 
these  loads  at  each  interval  was  evaluated,  and  taken  as  the  average  load  at 
the  corresponding  end  slip  interval  based  on  the  peak  point  of  the  average 
curve.  The  same  interval  was  used  for  the  descending  (post-peak)  branch 
of  the  average  curve.  Examples  of  series  of  four  pull-out  tests  and  their 
average  curve  (in  dotted  line)  are  shown  in  Fig.  6.5. 
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end  SLIP  (in.) 


END  SLIP  (in.) 


Fig.  6.5  -  Pull-out  tests  of  the  four  fibers  and  the  average  curve 
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6.6  PULL-THROUGH  VERSUS  PULL-OUT  TESTS 

In  the  pull-through  tests  the  fiber  extends  below  the  bottom  surface 
of  the  specimen  so  that  the  embedment  length  remains  constant  during  the 
pull-out  process.  For  the  pull-out  tests,  the  fiber  does  not  extend  below  the 
bottom  surface  of  the  specimen  and  thus  the  emdedment  length  decreases 
during  pull-out  . 

The  differences  in  loads  measured  in  the  pull-out  and  pull-through 
tests  for  1"  embedment  length  were  not  significant  for  end  slips  less  than 
about  0.2"  .  To  show  the  difference  between  these  two  types  of  tests  at 
larger  slips,  pull-out  and  pull-through  loads  were  plotted  versus  the  end 
slip  for  the  entire  embedment  length  ,  i.e.  up  to  1"  embedment  length  . 
Figs.  6.6  and  6.7  show  some  typical  results.  It  can  be  observed  that  there 
is  some  decay  in  the  response,  even  when  the  embedded  length  remains 
constant. 

6.7  PULL-OUT  LOAD  VERSUS  END  SLIP  RELATIONSHIPS 
6.7.1  Smooth  Fibers 

For  smooth  fibers  the  load  was  generally  linear  up  to  a  load  termed 
Peril  in  Chapter  V-  ^^^e  value  of  however  was  close  to  Pp,  the  peak 
load.  After  Pp,  the  load  dropped  quickly,  the  higher  the  matrix  strength, 
the  faster  the  drop  (  up  to  60%  drop  in  the  peak  load  at  full  debonding  in 
some  series).  As  shown  in  Fig.  6.8a,  the  descending  branch  of  the  curve 
indicates  that  debonding  in  smooth  fibers  is  faster  than  for  other  types  of 
fibers.  The  rate  of  debonding  is  negative  (positive  curvature),  meaning  that 
as  the  slip  increases,  the  bond  decreases  drastically. 

In  the  case  of  pull-through  tests,  after  full  debonding  has  occurred, 
the  load  approaches  a  constant  value  at  small  slips.  However  large  slips 
tend  to  induce  decay  at  the  interface  between  the  fiber  and  the  matrix  and 
the  concrete  tunnel  around  the  fiber  surface  is  damaged.  In  the  case  of 
pull-out  tests,  frictional  forces  decrease  not  only  because  of  decay  but  also 
because  the  embedment  length  decreases. 


PULL-THROUGH  LOAD  (lbs.)  21  PULL-OUT  LOAD  (lbs.) 
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g.  6.6  -  Pull  -  out  test  of  a  smooth  fiber  up  to  1”  end  slip  . 


Fig.  6.7  -  Pull-through  test  of  a  smooth  fiber  up  to  1"  end  slip  . 
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^  Pc  r  1 1 


SMOOTH  FIBERS 


FULL  DEBONDING 


Fig.  6.8a 


Ascending  Branch 


Motrlx  Strength  >  High 
Pud-out  Fiber  i  Smooth 
Fiber  Volume  Fraction  =  S/i 
Additives  •  Latex 
Fiber  Diometer  =  0.018* 

_ Actual  Test  Data 

(Including  fiber  elongation) 

-  Without  Fiber  Elongation 


END  SLIP  (in.) 


Fig.  6.8b  -  Ascending  branch  of  pull-out  curve  of  snwwth  fiber  with  and  without 


elastic  strains. 
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In  the  ascending  branch  it  is  observed  that  the  elastic  strain 
contribution  to  the  measured  slip  was  significant  especially  for  small  fiber 
diameters  (0.01")  .  This  contribution  decreased  when  larger  diameters 
were  used  although  the  pull-out  load  increased  .  Fig.  6.8b  illustrates  the 
pull-through  curve  when  elastic  strains  are  accounted  for  ,  i.e,  the  elastic 
elongation  at  the  free  end  is  substracted  from  the  observed  slip. 

Another  source  of  error  is  due  to  the  fact  that  the  smooth  fibers, 
being  cut  from  spring-type  steel  wire,  were  not  perfectly  straight  ,  i.e.  (had 
some  curvature  ).  Such  imperfections  in  straightness  or  at  the  fiber  tip 
which  becomes  flat  due  to  cutting  did  influence  the  pull-out  load  versus 
slip  measurements  especially  at  large  slips  . 

Table  6.5  shows  some  values  of  (P/A),  the  initial  slope  of  the  load- 
slip  curve  for  different  types  of  specimens  in  which  the  pull-out  fibers 
were  smooth.  These  values  range  from  25,000  Ib/in  to  68,000  Ib/in, 
depending  on  the  type  of  specimen.  In  these  results,  the  value  of  ^  is 
corrected  to  account  for  the  elastic  strain  contribution. 


Table  6.5  -  Pull-out  Data 


Series 

-Vl03  Ibs/in) 

Pp  (lbs.) 

Ap  (10-4  in.) 

H2SF 

44.0 

13.4 

7.80 

LOSN 

52. 1 

12.6 

7.30 

AISN 

55.0 

10.8 

6.34 

H2SL 

67.8 

22. 1 

15.9 

H2SM 

30.0 

16.8 

9.48 

HOSN 

25.0 

15.6 

9.25 

L2SN 

31.0 

10.3 

6.01 

AOSN 

45.0 

12.6 

7.18 

HOSL 

40.8 

12.1 

7. 12 

H3SN 

35.0 

17.2 

4.33 

A2SM 

28.0 

12.5 

2.72 

H2SN 

40.0 

18.9 

18.2 
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The  readings  of  the  bottom  LVDT  showed  that  the  load  at  which 
the  fiber  starts  to  pull  out  is,  in  general,  slightly  less  than  the  maximum 
observed  load.  Table  6.6  summerizes  the  values  of  the  peak  load  and  the 
load  at  which  the  second  LVDT  starts  to  move  with  their  corresponding 
displacements  as  obtained  from  tests  not  including  the  effect  of  fiber 
elongation.  In  Table  6.6,  ^2  is  the  end  slip  of  the  first  LVDT  for  the  load 
P2  described  in  column  4.  Fig.  6.9  shows  the  ascending  branch  up  to  the 
peak  load  on  a  larger  scale  for  smooth  fibers  as  obtained  from  tests  without 
accounting  for  fiber  elongation. 


Table  6.6  -  Loads  and  displacements  of  smooth  fibers. 


Series 

Pmax(lbs.) 

Amax  (in.) 

P2  (lbs.) 

Az  (in.) 

AOSL 

21.7 

0.00382 

19.3 

0.00419 

AOSN 

13.4 

0.00294 

13.0 

0.00310 

AOSN  D=0.03 

20.6 

0.00260 

A2SN 

16.7 

0.00291 

A3SN 

18.8 

0.00275 

15.9 

0.00322 

AOSN  D=0.01 

10.6 

0.00717 

10.2 

0.00754 

AOSN  L=0.5" 

10.5 

0,00226 

9.8 

0.00226 

HOSM 

14.6 

0.0030 

13.7 

0.00322 

Sifcon  D=0.0l 

12.3 

0.00921 

11.5 

0.00921 

sifconD=0.019 

22.0 

0.00464 

16.0 

0.00515 

Sifcon  D=0.03 

30.2 

0.00159 

28.5 

0.00196 

To  account  for  fiber  elongation  during  pull-out,  the  elastic  elongation 
Ae  =  (P*L)/(Af*Es)  should  be  substracted  from  the  values  of  Amax  and 
A2.  For  example  ,  for  series  AOSL,  P=21.7  lbs.,  L=  0.5",  Af=  0.000283 
in2,  Es  =  29000  ksi  =>  =  0.00132  in,  thus  the  corrected  value  of  Amax 

is  0.00382  -  0.00132  =  0.0025  in.  It  should  be  noted  here  that  the  values 
of  A  given  in  Tables  6.7  and  6.8  do  not  include  the  effect  of  fiber 
elongation. 


PULL-THROUGH  LOAD  (lbs.)  PULL-DUT  LOAD  <lbs) 
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.000  .ooa  .004  .006  .008  .010 


END  SLIP  (In.) 


.000  .008  .004  .006  .008  .010 


END  SLIP  (in.) 


Fig.  6.9  -  Experinicntal  pull-out  curves  of  smooth  fibers  up  to  the  peak  load  . 
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6.7.2  Deformed  Fibers 

Typical  puil-out  load  versus  end  slip  curves  of  deformed  fibers  are 
shown  in  Fig.  6.10.  It  can  be  observed  that  the  ascending  branch  is  made 
out  of  two  parts.  The  first  part  is  linear  and  the  second  part  is  highly 
nonlinear.  It  is  observed  that  for  the  second  part,  the  resistance  to  pull-out 
is  due  not  only  to  the  interfacial  bond  between  the  fiber  and  the  matrix 
(similar  in  nature  to  that  found  in  the  case  of  smooth  fibers),  but  due  to  an 
additional  resisting  force  provided  by  the  mechanical  deformations  on  the 
fiber  surface.  Fig.  6.11  shows  the  ascending  part  of  a  pull-out  curve  on  a 
larger  scale  for  a  typical  deformed  fiber,  while  Fig.  6.12,  shows  the  pull¬ 
out  curve  for  a  deformed  fiber  up  to  an  end  slip  of  0.25"  .  Following  the 
peak  load,  the  pull-out  load  starts  to  drop.  However,  in  this  case  the  drop 
is  different  from  that  seen  in  the  case  of  smooth  fibers.  This  difference  is 
due  to  the  fact  that  for  deformed  fibers,  as  the  slip  increases,  the  smooth 
portion  of  fiber  surface  tends  to  debond  fast  (  as  in  the  case  of  smooth 
fibers)  while  the  mechanical  indentations  on  the  fiber  surface  try  to 
displace  the  concrete  encasing  these  deformations.  When  that  concrete  is 
displaced  or  crushed,  the  load  drops  drastically  to  values  close  to  zero. 
This  can  be  explained  by  the  fact  that  when  the  concrete  in  contact  with  the 
fiber  surface  has  been  highly  disturbed,  it  could  no  longer  provide  the 
frictional  resistance,  as  the  one  observed  in  the  case  of  smooth  fibers. 
However  as  soon  as  the  deformations  enter  another  section  of  the  tunnel  (of 
different  cross  section  ),  the  load  picks  up  again  and  a  cyclic  response  is 
observed.  This  is  illustrated  in  Figs.  6. 13  and  6. 14. 

Pull-out  tests  performed  on  deformed  fibers  showed  that  these  fibers 
would  reach  ultimate  strains  and  fail  at  loads  around  45  lbs.  Many  fibers 
failed  when  such  load  levels  were  reached  or  exceeded.  Deformed  fibers 
whose  diameter  is  0.018",  exhibit  very  high  stresses  before  debonding  and, 
as  a  result,  failure  may  occur  before  full  debonding  or  the  start  of  pull-out. 
An  optimum  case  will  be  when  the  fiber  dtbonds  completely  and  starts 
pulling-out  just  before  the  ultimate  strength  of  the  fiber  is  reached.  Data 
from  tests  on  deformed  fibers  [not  including  the  effect  of  fiber  elongation] 
showed  that  the  displacements  at  maximum  loads  were  15  to  20  times  those 
of  smooth  fibers.  If  fiber  elongation  is  accounted  for,  this  value  could 


PULL-  THROUGH  LOAD 


DEFORMED  FIBERS 


Fig.  6.10  -  Typical  pull-out  load  versus  end  slip  of  deformed  fibers 


PULL-DUT  LOAD  <lb.)  PULL-OUT  LOAD  (lb.) 
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Matrlx  Strengthi  Low 
PULL-DUT  Fiber'  Deforned 
Fiber  Volume  Fraction'  S’/. 
Additives'  None 


.04  ,06 

END  SLIP  (in.) 


Fig.  6.11  '  Pull-out  curve  of  deformed  fiber  up  to  the  peak  load 


Matrix  Strength'  Low 
PuU-Dut  Fiber'  Deforned 
Fiber  Volune  Fraction'  S'/. 
Additives'  None 


END  SLIP  (in.) 


Fig.  6.12  -  Pull-out  load  versus  end  slip  of  deformed  fibers 
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exceed  50  .  Table  6.7  shows  values  of  loads  and  displacements  of 
deformed  fibers  at  peak  load  and  the  load  at  which  the  second  LVDT 
(attached  to  the  free  end)  starts  to  move. 


Table  6.7  -  Loads  and  displacements  of  deformed  fibers. 


Series 

P2  (lbs.) 

mam 

Pmax  (lbs.) 

^max  (in.) 

AODN 

11.5 

0.00975 

35.4 

0.0511 

AIDN 

16 

0.00341 

33.8 

0.0822 

L2DN 

18 

0.00335 

43.3 

0.0660 

L3DN 

12 

0.00792 

40.0 

0.0423 

LOOM 

6.5 

0.00171 

21.6 

0.0674 

The  pull-out  behavior  of  deformed  fibers  for  large  slips  is  shown  in  % 

Figs.  6. 13  and  6. 14.  In  Fig.  6. 14,  which  shows  a  pull-through  curve,  it  is 
observed  that  the  curve  repeats  itself  every  0.2”  over  the  descending 
branch.  This  corresponds  to  the  size  of  the  deformed  portion  (indentation) 
of  each  segment  of  the  fiber.  The  load  cycles  reach  about  the  same  peak  • 

value  and  have  the  same  shape.  This  can  be  explained  as  follows:  as  slip 
increases,  the  already  disturbed  concrete  by  the  previous  indentation  is 
again  disturbed  by  the  following  indentation.  The  total  resisting  force  is 
thus  approximately  the  same.  • 

Fig.  6.13  shows  a  pull-out  test,  where  the  repeated  curves  are 
decreasing  in  amplitude.  This  decrease  is  due  to  the  fact  that  the  number  of 
deformations  decreases  with  a  decrease  in  the  embedment  length  during  ^ 

pull-out  . 


A  final  observation  on  deformed  fibers,  is  that  the  deformations  (or 
indentations)  on  the  fiber  surface  were  not  exactly  uniform  for  all  batches 
of  fibers  used.  For  some  fibers  surface  deformations  were  not  deep 


PULL-THROUGH  LOAD  (lbs.) 
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cnough  .  As  a  result  of  this,  these  fibers  were  pulled  out  rather  easily 
from  the  matrix  . 

6.7.3  Hooked  Fibers 

A  typical  pull-out  load  versus  end  slip  curve  for  hooked  fibers  is 
shown  in  Fig.  6.15.  As  was  the  case  for  deformed  fibers,  the  ascending 
branch  is  made  up  of  a  linear  and  a  non-linear  part.  Figs.  6. 16,  6. 17,  and 
6.18  show  actual  pull-out  load  versus  slip  behavior  of  hooked  fibers 
respectively  up  to  the  peak  load,  up  to  0.25  "  end  slip  .and  up  to  0.5  "  end 
slip  .  Fig.  6.19  shows  a  typical  pull-out  load  versus  end  slip  curve  up  to  1" 
end  slip. 

Although  hooked  fibers  used  in  the  tests  have  a  different  diameter 
than  either  smooth  or  deformed  fibers,  it  is  still  possible  to  compare  their 
behavior  with  the  smooth  or  the  deformed  fibers.  The  ascending  branch 
had  generally  a  steeper  slope,  due  to  the  presence  of  the  hook,  which  pro¬ 
vides  a  considerable  resisting  force  to  pull-out.  As  the  hook  is  being 
pulled  out  from  the  matrix,  it  pushes  on  the  surrounding  concrete  and  tends 
to  straighten  up,  hence  causing  the  slope  of  the  ascending  branch  to  start 
decreasing.  When  the  hook  has  been  partially  straightened  and  debonding 
has  occurred,  the  pull-out  load  starts  to  drop.  As  the  slip  increases,  the 
load  continues  to  drop  as  the  partially  straightened  fiber  tries  to  slip  further 
in  the  surrounding  concrete  tunnel.  Therefore,  the  shape  of  a  portion  of 
the  descending  branch  reflects  the  case  of  a  partially  straightened  fiber 
trying  to  pull  out,  hence  inducing  frictional  forces  along  its  surface.  When 
the  hook  becomes  straight,  the  descending  branch  drops  fast  and  the  last 
lower  portion  of  the  curve  reflects  mainly  the  frictional  forces  between  the 
fiber  surface  and  the  matrix. 

The  maximum  pull-out  load  for  hooked  fibers  varied  from  50  lbs.  to 
105  lbs.  depending  on  the  matrix  strength.  Displacements  at  peak  loads 
were  10  to  15  times  those  of  smooth  fibers,  not  including  the  effect  of  fiber 
elongation.  If  fiber  elongation  is  accounted  for, this  value  could  exceed  30. 
On  the  other  hand,  displacements  at  peak  loads  of  hooked  fibers  were  less 
than  those  of  deformed  fibers  .  Table  6.8  shows  the  peak  loads  and  the 
corresponding  displacements  of  hooked  fibers. 
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PULL-DUT  LOAD  (lbs.)  PULL-DUT  LOAD  (lbs,) 

20.  40.  60.  80.  100.  H  ^0'  QQ' 


Fig.  6. 16  -  Pull-out  curve  of  hooked  fibers  up  to  the  peak  load 


Fig.  6. 17  -  Pull-out  curve  of  hooked  fibers  up  to  0.25  "  end  slip  . 


PULL-DUT  LOAD  (lbs.)  PULL-DUT  LOAD  (Ibs.) 

10.  40.  60.  80.  100.  120.  20,  40.  60.  80,  100. 
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0  .1  .2  .3  .4  .5  .6 


END  SLIP  (in.) 

Fig.  6.18  -  Pull-out  curve  of  hooked  fibers  up  to  0.5"  end  slip  . 
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Table  6.8  -  Loads  and  displacements  of  hooked  fibers. 


Series 

Pmax  (lbs.) 

HOHN 

102.8 

0.0288 

AOHN 

80.3 

0.0354 

LOHN 

58.9 

0.030 

AlHN 

91.7 

0.0478 

LIHN 

48.8 

0.0437 

L3HN 

66.1 

0.0354 

HOHN  L=0.5" 

89.2 

0.0332 

AOHN  L=0.5" 

72.6 

0.0411 

It  was  observed  in  some  experiments  that  the  hook  had  broken  inside 
the  matrix  before  straightening  .When  the  hook  breaks  the  load  drops 
faster  and  the  pull-out  behavior  afterwards  resembles  that  of  a  smooth 
fiber.  Fig.  6.20  shows  end  hooks  before  and  after  pull-out.  Fig.  6.21(a) 
shows  three  pull-out  tests  of  hooked  fibers  up  to  1"  end  slip  while  figure 
6.21  (b)  shows  the  pull-out  load  versus  end  slip  of  hooked  fibers  where 
two  fibers  had  broken  inside  the  matrix. 


6.8  PULL-OUTWORK 

The  pull-out  work,  or  the  dissipated  bond  energy  is  defined  as  the 
area  under  the  pull-out  curve.  This  area  was  found  to  be  much  smaller  in 
the  case  of  smooth  fibers,  than  for  deformed  or  hooked  fibers. 

Typical  Pull-out  work  values  calculated  for  the  three  types  of  fibers 
from  test  data  are  shown  in  Table  6.9  .  The  table  shows  pull-out  work 
values  for  pull-out  and  pull-through  tests  at  peak  loads,  at  0.2"  end  slip, 
and  at  1"  end  slip. 


(  b  )  Straightened  hook. 


(c) 


Broken  hook. 


Fig.  6.20  -  End  hooks  before  and  after  pull-out. 


PULL-DUT  LOAD  (lbs.)  PULL-DUT  LOAD  (lbs.) 

pn  40.  60.  80.  100.  80.  40.  60.  80.  100. 
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0  .8  .4  .6  .8  1.0 


END  SLIP  (in.) 

(a) 


.00  .05  .10  .15 

EMD  SLIP  (in.) 
(b) 


Fig.  6.21  -  Pull-out  load  versus  end  slip  of  liooked  fibers 


.20 


.25 
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Tablc  6.9  -  Pull-out  work  for  different  types  of  fibers  . 


Series 

Test 

Wmax(lb-in) 

W  0.2"  (lb-in) 

W  l"(lb-in) 

H3SN 

0.0239 

1.25 

H2SN 

0.0481 

1.21 

- 

H2SM 

Pull-out 

0.0473 

0.925 

- 

LOSN 

0.0178 

0.761 

- 

HISN 

0.0243 

1.21 

4.29 

AOSN 

Pull-out 

0.0336 

1.29 

4.06 

AOSL 

Pull-out 

0.0694 

0.967 

4.40 

AOSN  D=o.oi" 

Pull-out 

0.0436 

0.761 

1.37 

Sifeon  D=o.oi 

0.0736 

1.32 

6. 12 

AOSN  D=o.03“ 

0.0342 

1.41 

- 

HIDN 

3.76 

5.13 

13.4 

AODN 

4.35 

6.03 

16.9 

L2DN 

2.19 

5.47 

14.8 

LOHN 

Pull-out 

1.21 

8.81 

18.7 

LIHN 

Pull-out 

1.69 

6.93 

L3HN 

Pull-out 

1.78 

8.56 

- 

AOHN 

Pull-out 

2.33 

10.9 

23.9 

AlHN 

Pull-out 

3.41 

- 

- 

AOHN  L=0.5" 

Pull-out 

2.29 

10.5 

16.8 

HOHN 

Pull-out 

3.06 

14.5 

29.2 

To  compare  pull-out  work  of  different  types  of  fibers,  pull-out  work 
values  of  some  scries  of  deformed  and  hooked  fibers  were  normalized  with 
respect  to  the  pull-out  work  values  of  smooth  fibers  at  peak  load,  at  0.2" 
end  slip,  and  at  1"  end  slip  .  Normalization  was  done  by  dividing  the  pull¬ 
out  work  of  deformed  and  hooked  fibers  by  the  pull-out  work  of  smooth 
fibers.  Table  6.10  shows  these  values. 


Table  6. 10  -  Pull-out  work  ratios 


Series 

Wmax/WmaxS 

W0.2"/W0.2"S 

Wi'vWfS 

L2DN 

123 

6.82 

4.03 

AODN 

129 

4.67 

4.  16 

LOHN 

67 

7.33 

4.  1 

AOHN 

69 

7.73 

4.2 

It  is  observed  in  Table  6. 10  that  the  pull-out  work  to  the  peak, 
for  smooth  fibers  is  about  0.7%  to  1.5%  of  the  work  to  the  peak  of 
deformed  or  hooked  fibers.  This  value  increases  to  about  15%  to  20%  at 
0.2  "  end  slip  and  to  about  25%  at  1"  end  slip. 

6.9  DEVELOPMENT  OF  BOND  SHEAR  STRESS  VERSUS  SLIP 
RELATIONSHIP  CURVES 

The  procedure  described  in  Chapter  V  was  used  to  predict  bond  shear 
stress  versus  slip  relationships  from  the  average  pull-out  curves.  Since  the 
development  of  the  theory  in  Chapter  V  was  based  on  the  assumption  that 
the  bond  shear  stress  versus  slip  curve  is  elastic  with  an  exponentially 
decaying  frictional  zone,  the  (t-S)  curves  predicted  had  the  same 

characteristic  shape.  Since  only  selected  pull-out  curves  were  used  to 
predict  the  (x-S)  curves,  the  bond-slip  curves  obtained  are  representative 

but  not  exact  for  all  cases.  A  more  accurate  way  of  getting  the  bond  slip 
curve  for  a  given  matrix-fiber  combination  would  be  to  perform  a  point  by 
point  analysis  of  the  pull-out  curve,  and  for  each  point  on  the  curve,  get  a 
corresponding  point  on  the  (t-S)  curve.  Such  a  method  would  be 

substantially  more  complicated  and  numerically  involved.  The  main 
parameters  found  for  individual  series  using  the  procedure  described  in 
Chapter  V,  namely  k,  Xmax>  ^  given  in  Table  6. 1 1. 
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Using  the  parameters  of  Table  6. 1 1 .  it  is  possible  to  predict  the  pull¬ 
out  load  versus  end  slip  as  described  in  Chapter  V.  The  predicted  curve  can 
then  be  compared  to  the  experimentally  observed  curve.  Such  comparisons 
are  shown  in  Figs.  6.22  to  6.28.  Figs.  6.28  (b),  (c),  and  (d)  show  a 
magnified  scale  of  the  ascending  branch  only.  It  can  be  observed  that  the 
predicted  curves  are  in  reasonably  good  agreement  with  the  experimental 
curves. 


Table  6.  1 1  -  Pull-out  Data 


Series 

KdO^lb/in^) 

^max  (psi) 

If  (psi) 

H2SF 

3.69 

255 

223 

0.053 

LOSN 

5.25 

211 

211 

0.910 

AISN 

5.00 

200 

180 

0.0337 

H2SL 

8.79 

1390 

259 

0.938 

H2SM 

1.51 

285 

275 

1.03 

HOSN 

0.841 

334 

257 

1.30 

L2SN 

1.65 

176 

172 

-6.50 

AOSN 

3.87 

210 

210 

-1.00 

HOSL 

3.13 

207 

202 

1.06 

H3SN 

2.21 

289 

289 

0.523 

A2SM 

1.24 

209 

209 

5.73 

H2SN 

3.00 

300 

280 

2.28 

6.  10  EFFECT  OF  FIBER  AND  MATRIX  PARAMETERS  ON  THE  LOAD- 
SLIP  RELATIONSHIP 

6. 10. 1  Effect  of  Matrix  Strength 

As  shown  in  Fig.  6.29,  for  the  case  of  hooked  fibers  with  different 
matrix  strength,  the  pull-out  load  increases  with  increasing  matrix  strength. 
For  smootn  fibers  this  increase  was  mostly  noticeable  in  the  ascending 
branch. 
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;.  6.22  -  Experimental  pull-out  curve  versus  analytical  prediction  of  pull-out 
curve  of  smooth  fibers  . 
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Fig.  6.23  -  Experimental  pull-out  curve  versus  analytu  '  prediction  of  pull-out 


curve  of  smootli  fibers  . 
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.6.24  -  Experimental  pull-out  curve  versus  analytical  prediction  of  pull-out 


curve  of  smooth  fibers . 
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Fig.  6.25  -  Experimental  pull-ou.  ,  rve  versus  analytical  prediction  of  pull-out 


curve  of  smooth  fibers  . 
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.6.26  -  Experimental  pull-out  curve  versus  analyUcal  prediction  of  pull-out 


curve  of  smooth  fibers . 
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Fig.  6.27  -  Experimental  pull-out  curve  versus  analytical  prediction  of  pull-out 


curve  of  smooth  fibers . 
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Fig.  6.28  -  Experimental  pull-out  curve  versus  analytical  prediction  of  pull-out 


curve  of  smooth  fibers  . 
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Fig.  6.28  -  Experimental  pull-out  curve  versus  analytical  prediction  of  pull-out 
curve  of  smooth  fibers  up  to  the  peak  load 
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Fig  6.29  -  Effect  of  matrix  strength  on  pull-out  load  end  slip  relationship. 
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6. 10.2  Effect  of  Fiber  Volume  Fraction 

The  effect  of  fiber  volume  fraction  is  shown  in  Fig.  6.30  .  This 
effect  seems  to  be  rather  insignificant.  It  can  be  observed  that  an  increase 
up  to  3%  in  fiber  volume  fraction  increases  the  pull-out  load  by  about  10%, 
while  it  leads  in  the  post-peak  range  to  a  slight  increase  in  the  fiber 
resistance  to  pull-out.  This  trend  was  not  clear  when  compared  to  the 
variability  of  results  observed.  It  seems  that  the  pull-out  fibers,  being 
much  smaller  in  diameter  than  ordinary  reinforcement,  probably  do  not 
cause  appreciable  damage  to  the  surrounding  concrete  upon  pull-out;  thus 
adding  some  fibers  to  the  matrix  may  not  add  much  to  the  interfacial  bond 
strength  between  the  fiber  and  the  matrix  . 

Wh;n  SIFCON  matrices  were  used  [  Vf  =  11%  ],  considerable 
improvement  in  the  interfacial  bond  strength  between  the  fiber  and  the 
matrix  was  observed,  both  before  debonding  and  after  debonding  .  An 
increase  of  20%  to  25%  in  the  pull-out  load  was  observed  before  debonding 
and  75%  to  80%  increase  was  observed  after  debonding  .  The  presence  of 
fibers  in  the  matrix  in  such  high  percentages  seems  to  provide  an  additional 
bond  component,  namely  fiber  interlock.  Such  an  arrangement  provides 
excellent  confinement  to  the  mortar  which  considerably  increases  the 
interfacial  bond  strength  and  the  pull-out  work.  Also  it  helps  reduce 
frictional  decay  in  the  descending  branch  of  the  curve.  Fig.  6.31  shows  a 
typical  pull-out  load  versus  end  slip  curve  using  plain  and  SIFCON 
matrices  . 

6. 10.3  Effect  of  Additives 

The  effect  of  additives  is  presented  in  Fig.  6.32.  In  terms  of 
maximum  pull-out  load  values.  Latex  was  found  to  yield  the  highest  values, 
followed  by  Fly  ash  .  Microsilica  did  not  seem  to  improve  the  interfacial 
bond  strength  and  the  test  series  was  repeated  twice  to  ascertain  this  result. 
With  Latex  a  50%  to  60%  increase  in  the  pull-out  load  was  observed  over 
plain  matrix;  this  remarkable  increase  in  pull-out  load  leads  to  a  rapid  drop 
in  the  load  at  debonding  as  was  observed  with  high  strength  matrices  . 
After  complete  debonding,  the  load  remains  almost  constant  in  the 


PULL-THROUGH  LOAD  (lbs.)  21  PULL-THROUGH  LOAD  (lbs.) 


-151- 
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6.30  -  Effect  of  fiber  volume  fraction  on  pull-out  load  slip  relationship  of  smooth 
fibers . 


Fig.  6.3 1  -  Effect  of  SIFCON  matrix  on  pull-out  load  slip  behaviour  of  smooth  fibers  . 
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Fig.  6.32  -  Effect  of  additives  on  pull-out  slip  relationship  of  smooth  fibers  . 
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frictional  phase.  Fly  ash  increased  the  maximum  pull-out  load  by  10%  to 
20%  and  was  more  effective  than  Latex  in  improving  friction  forces. 

6. 10.4  Effect  of  Fiber  Diameter 

Pull-out  load  versus  end  slip  curves  for  fibers  of  different  diameters 
are  shown  in  Fig.  6.33  .  A  change  in  fiber  diameter  did  not  affect  the  pul- 
out  load  -end  slip  behavior  of  smooth  fibers.  Elastic  strains  in  the  free 
portion  of  the  fiber  decreased  as  fiber  diameter  was  increased  although  the 
pull-out  load  increased  .  It  was  observed  that  the  load-slip  measurements 
became  more  sensitive  to  fiber  curvature  as  the  diameter  increases. 

6.  10.5  Effect  of  Fiber  Embedment  Length 

The  effect  of  fiber  embedment  length  on  bond  strength  is  shown  in 
Figs.  6.34  ,  6.35,  and  6.36.  It  was  observed  that  the  ascending  branch  of 
the  curve  is  less  sensitive  to  the  fiber  embedment  length  than  the 
descending  branch  for  smooth  fibers.  The  pull-out  load  end  slip 
relationship  exhibited  the  same  behaviour  in  both  cases  for  the  three 
different  types  of  fibers.  In  hooked  fibers,  since  the  hook  contributed  most 
to  the  pull-out  load,  changing  the  embedment  length  from  0.5"  to  1"  did 
little  to  the  ascending  portion  of  the  curve.  However,  in  deformed  fibers  , 
since  the  number  of  surface  indentations  increase  with  the  increase  in 
embedment  length,  changing  the  embedment  length  from  0.5"  to  1" 
approximately  doubled  the  value  of  loads  in  the  ascending  branch. 

6. 10.6  Contribution  of  Surface  Indentations  and  End  Hooks 

6.10.6.1  Surface  Indentations.  In  one  series  of  deformed  fibers,  the  pull¬ 
out  fibers  were  greased  on  the  surface  to  isolate  the  effect  of  surface 
indentations.  The  average  pull-out  curve  versus  end  slip  of  this  series  was 
then  superimposed  on  the  pull-out  curve  versus  end  slip  of  a  smooth  fiber 
of  same  diameter.  The  resulting  pull-out  curve  versus  end  slip  curve  was 
compared  to  the  experimentally  observed  curve  of  deformed  fibers.  Such  a 
comparison  is  shown  in  Fig.  6.37.  It  is  observed  that  the  two  curves  are 
quite  similar  .  The  difference  between  the  two  curves  can  be  attributed  to 
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g.  6.33  -  Effect  of  fiber  diameter  on  pull-out  slip  relationship  of  smooth  fibers  . 
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Fig.  6.34  -  Effect  of  fiber  embeddment  length  on  pull-out  load  slip  behaviour  of  smooth 
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Fig.  6.35  -  Effect  of  fiber  embeddment  length  on  pull-out  load  slip  behaviour  of 
deformed  fibers  . 


Fig.  6.36  -  Effect  of  fiber  embeddment  length  on  pull-out  load  slip  behaviour  of 


hooked  fibers  . 
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Fig.  6.37  -  Experimental  pull-out  curve  versus  superposition  pull-out 
curve  of  deformed  fibers. 
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Fig.  6.38  -  Experimental  pull-out  curve  versus  superpositior  /oll-out 
curve  of  hooked  fibers  . 


-157- 


the  variability  in  experimental  data  and  to  the  non-  uniformity  of  grease 
distiibution  on  the  fiber  surface. 

6.  10.6.2  End  Hooks.  In  one  series  of  hooked  fibers,  the  pull-out  fibers 
were  greased  on  the  surface  to  isolate  the  contribution  of  the  end  hook  to 
the  pull-out  load.  The  average  pull-out  load  versus  end  slip  of  this  series 
was  then  superimposed  on  the  pull-out  load  versus  end  slip  curve  of  a 
smooth  fiber  of  same  diameter.  Similar  to  deformed  fibers,  the  resulting 
curve  was  compared  to  the  experimentally  observed  curve.  Such 
comparison  is  shown  in  Fig.  6.38.  Here  also,  the  two  curves  are  in 
reasonably  good  agreement. 

6. 10.7  Stress  Levels  in  Pull-out  Fibers 

The  pull-out  load  of  the  average  curve  of  one  series  of  smooth,  deformed, 
and  hooked  fibers  was  divided  by  the  corresponding  fiber  area  [(n*d*d)/4] 
to  obtain  the  fiber  axial  stress.  Fig.  6.39  shows  the  values  of  fiber  stresses 
versus  slip  for  the  three  types  of  fibers.  It  is  observed  that  the  stress  levels 
in  smooth  fibers  are  much  smaller  than  those  of  deformed  or  hooked  fibers. 
Many  deformed  fibers  failed  during  pull-out  because  stress  levels  were 
close  to  the  failure  stress  of  the  fiber  (0~c>u  =  170  ksi).  For  hooked  fibers  , 
the  fiber  stress  varied  from  100  to  150  ksi  approximately. 


6.  1 1  GENERAL  OBSERVATIONS  AND  CONCLUSIONS 

From  the  results  of  this  limited  experimental  study,  the  following 
conclusions  can  be  drawn: 

1.  Hooked  fibers  and  deformed  fibers  have  higher  resistance  to  pull¬ 
out  than  smooth  fibers  because  of  the  mechanical  contribution  of 
end  hook  and  surface  deformations  respectively.  At  larger  slips, 
deformed  fibers  provide  oscillatory  resistance  to  pull-out  ,  while 
the  resistance  of  hooked  fibers  continues  to  decrease  as  the  end 
hook  straightens  out  .  Smooth  fiber  resistance  to  pull-out  also 
decreases  but  at  a  faster  rate  than  deformed  or  hooked  fibers  . 
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Fig.  6.39  -  Stress  levels  in  pull-out  fibers. 
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2.  As  matrix  strength  increases,  the  bond  between  the  fiber  and  the 
matrix  increases.  However,  debonding  following  the  peak  load  is 
much  faster  for  high  strength  matrices  . 

3.  Fiber  volume  fraction  in  the  matrix  has  little  effect  on  the  bond 
between  the  fiber  and  the  matrix  .  An  increase  of  10%  to  15%  in  the 
maximum  pull-out  load  was  observed  as  the  fiber  volume  fraction 
increases  from  0%  to  3%  .  A  similar  trend  was  observed  in  the  post 
peak  range. 

4.  The  presence  of  latex  emulsion  and  fly  ash  in  the  matrix,  was  found 
to  improve  the  bond  at  the  fiber  matrix  interface.  However,  latex 
was  much  more  effective  than  fly  ash.  The  addition  of  microsilica 
did  not  lead  to  any  noticeable  improvement  in  ihe  bond  strength  at 
the  interface  . 

5.  The  pull-out  work  of  smooth  fibers  is  significantly  smaller  than 
that  of  deformed  or  hooked  fibers.  It  was  also  observed  that,  in 
the  case  of  smooth  fibers,  the  amount  of  pull-out  work  to  the  peak 
load  was  a  small  fraction  of  the  total  pull-out  work.  This  means 
that,  in  the  case  of  smooth  fibers,  the  work  provided  by  frictional 
forces  is  larger  than  that  provided  by  the  forces  of  initial  adhesion. 
For  hooked  or  deformed  fibers  on  the  other  hand,  the  fraction  of 
the  pull-out  wor'"  to  the  peak  load  was  larger  than  for  smooth 
fibers,  and  the  work  contributed  by  the  mechanical  forces  was 
doir’nant.  Everything  else  being  equal,  the  pull-out  work  up  to  the 
complete  pull-out  of  smooth  fibers  is  only  about  15%  to  20%  of 
that  of  hooked  or  deformed  fibers. 

6.  Bond  shear  stresses  increase  with  an  increase  in  matrix  strength, 
and  are  affected  by  the  type  of  pull-out  fiber  used.  For  smooth 
fibers,  maximum  bond  stresses  were  found  to  range  from  176  to 
334  psi,  and  frictional  stresses  ranged  from  172  to  289  psi.  A 
maximum  bond  stress  of  up  to  1390  psi  was  observed  when  latex 
was  used. 
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The  fiber  diameter  and  fiber  embedment  length  did  not  affect  the 
characteristics  of  the  pull-out  load  versus  end  slip  relationship. 

The  bond  shear  stress  slip  relationship  of  hooked  and  deformed 
fibers  can  be  decomposed  into  two  parts:  one  due  to  the  interfacial 
bond  of  the  smooth  surface  of  the  fiber,  and  the  other  is  provided 
by  the  end  hooks  in  case  of  hooked  fibers  or  by  surface 
indentations  in  the  case  of  deformed  fibers. 


9 


-161- 


Appendix  VI  A 


"Additional  Experimental  Load-Slip  Curves" 


Notes 

1-  For  all  Smooth  Fibers,  the  Diameter  is  0.019”  and  the 
Embedded  Length  is  1"  except  where  indicated  otherwise. 

2-  For  all  Deformed  Fibers,  the  Diameter  is  0.018"  and  the 
Embedded  Length  is  1"  except  where  indicated  otherwise. 

3-  For  all  Hooked  Fibers,  the  Diameter  is  0.03"  and  the 
Embedded  Length  is  1"  except  where  indicated  otherw' 
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CHAPTER  VII 


BOND  STRESS  MODEL  FOR  FIBER-REINFORCED  CONCRETE  BASED 
ON  BOND  STRESS-SLIP  RELATIONSHIP 


7.1  INTRODUCTION 

There  is  a  basic  need  for  the  study  of  bond  in  fiber-reinforced  ce¬ 
mentitious  composites.  Materials  scientists  and  engineers  have  been  aware 
of  the  importance  of  bond;  the  problem  has  indeed  been  addressed  before. 
However,  few  thorough  studies  have  been  conducted  to  date,  and  to  the 
best  of  the  authors'  knowledge,  no  complete  analytical  model  of  bo"d  in 
cementitious  composites  has  been  developed. 

Existing  studies  are  incomplete,  in  a  sense  that  all  the  aspects  of 
bond  were  not  investigated.  Some  authors  assumed  non-linear  elastic  shear 
bond  without  debonding  [27,  42,  48,  54],  hence  ignoring  the  frictional 
bond  stresses  that  develop  after  debonding  occurs.  Aveston  et  al.  [3] 
assumed  frictional  shear  bond  with  no  debonding,  hence  ignoring  the 
elastic  bond  that  develops  for  shear  stresses  not  exceeding  the  bond 
capacity  of  the  interface. 

To  incorporate  both  the  adhesion  and  the  mechanical  components  of 
bond  introduced  crimped  or  deformed  fibers  requires  a  more  comprehensive 
study  of  bond.  Bond  studies  undertaken  on  reinforced  concrete  cannot  be 
directly  applied  to  fiber-reinforced  cementitious  composites.  The 
mechanical  and  adhesive  bond  have  the  effect  of  changing  both  the 
distribution  and  the  average  value  of  the  shear  bond  at  the  interface. 
Everything  else  being  equal,  a  crimped  fiber  for  example  can  develop  more 
bond  stresses  at  the  interface  than  a  smooth  fiber.  Likewise,  a  matrix  con- 
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taining  latex  could  develop  more  adhesion  to  the  fiber,  and  thus  induce 
more  bond  shear  stresses.  The  matrix  composition  is  expected  to  affect  the 
bond  properties  of  fiber-reinforced  cementitious  composites. 

Two  different  analytical  models  of  interfacial  bond  are  derived  in 
Chapters  7  and  8.  The  model  presented  in  this  chapter  is  based  on  an 
assumed  known  bond  shear  stress-slip  relationship,  which  is  used  to 
analytically  predict  the  interfacial  bond  shear  stress,  as  well  as  the  normal 
stress  distributions  in  both  the  fiber  and  the  matrix,  when  a  fiber-re¬ 
inforced  concrete  specimen  is  subjected  to  tensile  stresses.  The  model 
assumes  perfect  alignment  of  the  fibers,  as  well  as  a  square  packing.  The 
solution  applies  to  the  case  where  a  crack  doesnot  occur  along  the  fiber  or 
prism  considered. 

7.2  BASIC  ASSUMPTIONS 

a)  The  main  assumption  needed  for  the  development  of  the  model  is 

the  shear  stress-slip  relationship.  As  can  be  seen  in  Fig.  7.1,  the 
assumed  relationship  is  bi-linear  or  elastic-perfectly  frictional.  A 
region  of  linear  elastic  shear  bond  prevails  until  debonding  occurs 
when  the  shear  stress  reaches  the  maximum  shear  capacity  of  the 
interface  at  which  point  the  stresses  developed  are  of 

frictional  nature,  i.e.  are  independent  of  the  amount  of  slip. 

b)  It  is  also  assumed  that  the  fibers  are  aligned  within  the  specimen, 
and  that  they  are  squarely  packed.  A  square  packing  means  that  at 
any  given  section  perpendicular  to  the  line  of  stress,  the  existing 
amount  of  fibers  is  the  same.  Also,  within  a  given  section,  the 
spatial  distribution  of  the  fibers  is  uniform.  Fig.  7.2  shows  a 
square  packing  case. 

c)  Furthermore,  it  is  assumed  that  each  individual  fiber,  along  with 
its  share  of  matrix  act  and  behave  independently  of  other  fibers 
and  the  rest  of  the  matrix  body.  This  assumption  is  essential  to 
simplify  the  analysis  of  the  whole  composite  down  to  that  of  a 
representative  unit.  The  representative  unit,  which  is  a  square 
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7.1  -  Bilinear  bond-slip  relationship. 
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Square  packing  configuration. 
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prism  with  a  fiber  along  its  axis,  is  further  simplified  to  a  cylin¬ 
der  of  matrix  encasing  one  fiber  (see  Fig.  7.3). 

d)  Finally,  the  tensile  load  acting  on  the  representative  unit  is  as¬ 
sumed  shared  by  the  fiber  and  the  surrounding  matrix. 


7.3  MATHEMATICAL  DERIVATION  OF  THE  SHEAR  STRESS 
DISTRIBUTION 

Using  the  bond  stress-slip  relationship  assumed  previously  (see  Fig. 
7.  1),  the  following  equation  holds  for  the  elastic  bond  range; 

T(x)  =  tc  S(x)  (7.1) 


The  variation  of  slip  S  is  equal  to  the  difference  between  the  strain  in 
the  fiber  and  the  strain  in  the  matrix.  Hence: 


dx 


ef(x)  • 

ff  (x) 
Ef 

F(x) 
Af  Ef 


em(x) 


fm(x) 


T(x) 

Em 


(7.2) 

(7.3) 

(7.4) 


where: 

ff(x)  =  local  stress  in  fiber 

fm(x)  =  local  stress  in  the  matrix 

F(x)  =  local  force  in  fiber 

T(x)  =  local  force  in  the  matrix 

Ef  =  modulus  of  elasticity  of  the  fiber 

Em  =  modulus  of  elasticity  of  the  matrix  in  tension 

Af  =  area  of  one  fiber 
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a.  Longitudinal  section  for 
representative  unit 


2 

(d  /  D)  =  Vf 


c.  Assumed  cross  section 
for  representative  unit 


Fig.  7.3  -  Typical  representative  unit. 
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An,  =  area  of  the  matrix  within  the  representative  unit 

Assuming  that  the  representative  unit  is  subjected  to  a  tensile  force 
Fq,  then: 

T(x)  =  Fo-F(x)  (7.5) 

therefore: 


F(x) 
Af  Ef 


[Fq  -  F(x)] 

Am 


(7.6) 


Af  Ef 


Am  Er 


Am  ^m 


(7.7) 


Since  a  square  packing  is  assumed,  the  following  relationship  is  true: 

where  Vf  is  the  fiber  volume  fraction  in  the  composite. 

It  can  be  shown  that  Eq.  7.8  can  be  rewritten  as: 


F(x) 
Af  Ef 


1  +  n 


Vf 

1  -  Vf 


n  F„  Vf 
Af  Ef  1  -  Vf 


where  n  is  the  modular  ratio,  or: 


Eq.  "'.9  c;?..  be  written  as  a  function  of  the  stress  in  the  fiber  ff: 


(x) 


1  +  n 


Vf  >1  .  JLla  Vf 
1-Vf  J  AfEf  1-Vf 


(7.10) 


Differentiating  both  sides  of  Eq.  7. 10  leads  to: 


^  _  -L  r 

dx2  -  Ef  [ 


1  +  n 


1  -  Vf 


dff(x) 

dx 


(7.11) 
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but 

where 


dff(x) 

dx 


thus; 


V 

T(x)  = 

d2s 

dx2  - 


t(x)  V 

Af 


fiber  perimeter 


local  bond  shear  stress 


1 

Ef 


1  +  n 


1  -  V 


f- 


T(X)  W 


Recalling  Eq.  5. 1,  Eq.  5. 13  can  be  rewritten  as: 
d2S  4  ^  ^  Vf_ 


dx2 


d  E 


7( 


I  -t-  n 


1  -  Vf. 


K  srx) 


where  d  is  the  fiber  diameter. 


This  differential  equation  is  of  the  form: 
S"  -  K2  S  =  0 


where: 


K2 


4  r  ,  Vf _ 

Ef  (  ^  "  1  -  Vf. 


The  solrtion  to  this  differential  equation  is  of  the  form: 


S(x)  =  A  sinh  (Kx)  +  B  cosh  (Kx) 


(7.12) 


(7.13) 


(7.14) 


(7.15) 


(7.16) 


(7.17) 


The  constants  are  determined  from  the  boundary  conditions.  If  the 
origin  is  to  be  taken  at  the  center  of  the  fiber,  then: 


1.  @  X  =  0  ,  S  =  0  ,  hence  B  =  0 


2.  @x  =  |,F  =  xFo 
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where  X  is  the  fraction  of  the  total  force  acting  on  the  representative  unit 
that  is  taken  by  the  fiber  tips. 


The  factor  %  is  reflective  of  the  way  the  force  Fg  is  applied  onto  the 
representative  unit.  A  critical  value  for  x.  Xcrit>  would  be  the  one  for 
which  the  force  distributions  in  the  fiber  and  the  matrix  are  constant,  i.e. 
under  equal  strain  conditions.  If  the  force  in  the  fiber  is  to  be  maximum  at 
the  tips  and  minimum  at  mid-length,  x  is  to  be  larger  than  Xcrit>  i^^Jt  always 
less  than  1,  or  Xcrit  ^  X  -  i- 


Xcrii  ~ 

n  Vf 

1  +  (n  -  1)  Vf 

(7.18) 

Since  t(x)  = 

K  S(x) 

= 

A  X  sinh  (Kx) 

(7.19) 

,  ,  dff(x) 

but  - 

x(x)  v 

Af 

therefore: 

dff(x) 

dx 

hence: 

4 

j  A  K  sinh  (Kx) 

(7.20) 

ff(x)  = 

4 

A  K  cosh  (Kx) 

(7.21) 

F(x)  = 

A  K  cosh  (Kx) 

(7.22) 

II 

iiu 

A  K  cosh  (k|^) 

(7.23) 

but  F(|)  =  xFo 
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thus:  A 


n  d 


X  Fq  K 
K  cosh(Kj) 


Therefore; 

x(x)  = 


F(x)  = 


and  T(x)  = 


X  Fq  K 


7t  d  coshfKj) 


\ 

sinh  (K  x) 

/ 


'  X  Fq 
cosh(K|-) 


cosh  (K  x) 

j 


Fq  -  F(x) 


coshfK^) 


cosh  (Kx) 

) 


n.2A) 


(7.25) 


(7.26) 


(7.27) 


Eq.  7.25  is  valid  so  long  as  the  bond  stress  does  not  exceed  the  bond 
capacity  of  the  interface,  t^ax  (the  bond  strength),  also  assumed  to  be  the 
frictional  strength  between  the  matrix  and  the  fiber.  Fig.  7.4  shows  typical 
force  distributions  within  the  fiber  and  the  matrix,  while  Fig.  7.5 
represents  the  variation  of  the  bond  distribution  with  increasing  applied 
stresses  up  to,  and  then  past  the  stress  that  would  create  a  shear  stress 
equal  to  It  also  shows  how  debonding  occurs  in  similar  cases.  If  for 

a  given  tensile  force,  the  maximum  theoretical  shear  stress  in  a  fiber  is 
referred  to  as  xi,  then: 


X  Fq  K 
K  d 


tanh  (Kj) 


(7.28) 


If  Xj  is  larger  than  Xn,ax»  ^hen  debonding  occurs,  at  a  distance  from 
the  center  of  the  fiber  on  each  side  (Fig.  7.5).  Over  the  still  bonded 
region,  Eq.  7.25  applies  with  two  modifications,  namely  that  the  length  is 
2[i  instead  of  C,  and  the  force  is  (x  Fq  -  t^ax  V  ^d)  instead  (x  Fg),  hence: 


X 


0 


(1  -Z)F„ 


Fig.  7.4  -  Typical  force  distributions  within  the  fiber  F(x)  and  within 
the  matrix  T(x). 


centerline  of  fiber 


T  >  T 


Fig  7  5  -  Typical  shear  stress  distributions  at  the  interface  between  the 
fiber  and  the  matrix,  before  and  after  debonding. 
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7.4  NUMERICAL  EXAMPLE 

Let  us  consider  the  ideal  case  of  a  fiber-reinforced  concrete  tensile 
specimen  that  meets  the  conditions  mentioned  in  the  development  of  the 
model  (perfect  alignment  as  well  as  square  packing,  ...).  The  prism  has  an 
assumed  square  cross  section,  3”  x  3"  (76.2  mm  x  76.2  mm).  The  force 
applied  is  P  =  4500  lbs  (20000  N),  the  modular  ratio  n  =8  ,  the  modulus  of 
elasticity  of  the  steel  fibers  Ef  =  29,000  ksi  (2x105  MPa),  the  fiber  length  [ 
=  2in.  (50.8  mm),  the  factor  %  =  0.5  (arbitrarily  chosen).  The  bond  stress 
-  slip  relationship  is  assumed  to  be  as  in  Fig.  7.6,  hence  the  bond  modulus 
K  =  40,000  ksi/in  (10900  MPa/mm),  the  bond  strength  Tn,ax  equals  the 
frictional  bond  tf  =  400  psi  (2.76  MPa). 

Case  1: 


In  this  case,  it  is  assumed  that  the  diameter  of  the  fibers  d  =  0.01  in 
(0.254  mm),  and  that  the  fiber  volume  fraction  Vf  =  Xcrii  “  0.075. 

First,  the  applied  stress  is  found: 


SHEAR  STRESS  (ksi) 
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SLIP  (10 in.) 


Fig.  7.6  -  Assumed  bond  shear  stress  versus  slip  relationship  in 
numerical  example. 


\ - 1 

Fiber  mid-length.  Fiber  tip, 

X  =  0  X  =  1  in. 


Fig.  7.7 


Predicted  bond  stress  distribution  -  Case  1. 
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=  =  500  psi  =  0.5  ksi  (3.45  MPa) 

(3  X  j)  in^ 

The  force  Fg  applied  on  the  representative  unit  is  hence  equal  to: 

F  -  ^ 

^0  -  4  Vf- 

=  — O.o"'!  ~  =  0.003927  Kips  (17.5  N) 

From  Eq.  7. 16  ,  K  is  such  that: 


4  /<  .  Vf  ^  K 

Ef  (  '  j  d 


4 

29000 


1  +  8 


1  -  0.01 


40000 

0.01 


K  =  24.42  in'*  (0.961  mm'*) 

The  bond  distribution  is  described  by  Eq.  7.24,  or: 


t(x)  = 


X  Fg  K 
jc  d  cosh(K^) 


sinh  (Kx) 


0.5x0.003927x24.42  ^ 
- 2“  (24.42x) 

TC  (0.01)  cosh(24.42x2) 


=  7.572  X  10'**  sinh  (24.42  x)  (ksi)  where  x  is  in  inche 

=  5.221  x  10'***  sinh  (0.961  x)  (MPa)  where  x  is  in  mm. 

To  check  if  debonding  has  occurred,  the  maximum  theoretical  bond 
stress  Xj  is  obtained  @  x  =  |,  or: 


7.572  X  10'**  sinh  (24.42x2)  =  1.53  ksi  (10.5  MPa)  > 


0.4  ksi  (2.76  MPa) 
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therefore  debonding  occurs,  over  a  length  =  2  To  obtain  the  value  of 
Ti  has  to  be  solved  for  as  in  Eq.  7.30: 


X  Fq  ~  "^max  V  ^2 


jc  d 


K  tanh  (K  fj)  =  x^ax 


^0.5x0.003927  -  0. 4xjtx0. 01  xfl  ^ 


or 


hence  f  t 


7t  X  0.01 
0.885  in  (22.5  mm). 


I 


x24.42  tanh(24. 42xft)  =  0.4 


=  j  -  0.885  =  0. 115  in  (2.92  mm) 

The  bond  shear  stress  distribution  over  half  the  fiber  length  is  shown 
in  Fig.  7.7. 

C.as.s.li 

The  diameter  of  the  fibers  d  is  now  assumed  to  be  0.02  in  (0.508 
mm),  while  the  fiber  volume  fraction  Vf  is  taken  as  10%.  Xcrit  =  0.47. 

4500  1^  =  500  psi  =  0.5  ksi  (3.45  MPa) 

in^ 


(3  X  3)  in^ 


The  force  Fg  is  now  equal  to: 

Fn  = 


d2  ^ 


Vf 


^  =  0.001571  Kips  (6.988  N) 


K2 


-  J 


K 

d 


29000 


1  +  8 


0.1 


1  -  0 


.1  J  0. 


40000 
02 
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or 


22.83  in-1  (0.899  mm'l) 


The  bond  stress  at  any  point  can  be  found  from; 


T(x)  = 


7t  d  coshfKj) 


sinh  (Kx) 


(  0.5x0.001571x22.83  ^ 

7Cx(0. 02)xcosh(22. 83x^) 


sinh  (22. 83x) 


6.943  X  lO'll  sinh  (22.83  x)  (ksi)  where  x  is  in 
inches. 

4.787  X  lO'lll  sinh  (0.899  x)  (MPa)  where  x  is  in  mm. 


The  maximum  theoretical  bond  stress  Ti  is  obtained  @  x  =  y .  or: 


Xi  =  6.943  X  10-11  sinh  (22.83x|)  =  0.285  kd  (1.97  MPaX  Xmax  =  0-^ 

ksi  (2.76  MPa) 
therefore  there  is  no  debonding. 

The  bond  shear  stress  distribution  over  half  the  fiber  length  for  Case 
2  is  shown  in  Fig.  7.8. 


7.5  CONCLUDING  REMARKS 

1.  The  analytical  bond  model  developed  in  this  chapter  for  fiber 
reinforced  concrete  composites  is,  to  the  authors'  knowledge,  the 
first  attempt  to  predict  analytically  the  bond  stress  distribution  at  the 
interface  between  fibers  and  the  surrounding  matrix  using 
experimentally  obtainable  data,  namely  the  bond  shear  stress-slip 
relationship. 


2.  The  model  predicts  the  shear  stresses  at  the  interface  at  any  section, 
as  well  as  the  force  distributions  within  the  fibers  and  the  matrix. 


Fiber  mid-length.  Fiber 

X  =  0  X  =  1 


Fig.  7.8  - 


Predicted  bond  stress  distribution  -  Case  2. 
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3.  The  method  followed  in  deriving  the  model  analytically  can  be 
applied  to  different  assumed  bond  stress-slip  relationships,  whether 
linear  or  nonlinear.  Experimentally  derived  bond  stress-slip 
relationships  can  be  easily  accomodated. 

4.  The  model  presented  incorporates  all  the  important  geometric  and 
mechanical  properties  of  both  the  fibers  and  the  matrix. 

5.  In  the  course  of  the  development  of  the  model,  no  properties  specific 
to  the  matrix  or  the  fiber  were  made  (except  that  they  both  behave 
linear-elastically ).  Hence,  the  same  model  can  be  used  for  all 
fibrous  composites. 

6.  The  model  can  be  generalized  to  the  case  of  random  distribution  and 
orientation  by  introducing  adequate  correction  factors  to  the  value  of 
the  volume  fraction  used. 
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chapter  VIII 

MODELING  OF  BOND  IN  FIBER-REINFORCED  CONCRETE  BASED  ON 
FORCE  TRANSMISSION  MECHANISM 


8.1  INTRODUCTION 

The  model  to  be  developed  in  this  chapter,  just  like  the  one  derived 
in  Chap.  7  analytically  predicts  the  interfacial  bond  shear  stress,  as  well  as 
the  normal  stress  distributions  in  both  the  fiber  and  the  matrix,  when  a 
fiber-reinforced  concrete  specimen  is  subjected  to  tensile  stresses. 
However,  while  the  model  in  Chap.  7  is  based  on  an  assumed  bond  shear 
stress  versus  slip  relationship,  the  model  of  this  chapter  uses  an  assumed 
mechanism  of  force  transmission  between  the  fibers  and  the  surrounding 
matrix  to  predict  the  bond  and  normal  stress  distributions  at  the  interface 
and  in  the  fibers  and  matrix.  In  both  models,  a  perfect  alignment  of  the 
fibers,  as  well  as  a  square  packing  are  assumed. 

This  model  is  furthermore  used  to  analytically  predict  the  interfacial 
bond  modulus  (ratio  of  the  local  shear  stress  and  local  slip  under  elastic 
conditions),  which  is  shown  to  be  location  independent.  An  expression  for 
the  pre-cracking  strength  of  a  cementitious  fiber  composite  is  also  derived 
and  parametrically  analyzed.  Finally,  the  debonding  stress  of  the 
composite  is  derived  from  the  model. 

The  model  developed  here  is  essentially  based  on  Yankelevsky' s 
work  [65],  with  the  proper  modifications  and  additions.  Yankelevsky' s 
work  was  centered  on  the  problem  of  axially  loaded  specimens,  where  a 
single  deformed  bar  is  encased  concentrically  in  a  long  concrete  cylinder 
with  the  bar  ends  exposed. 
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The  present  model  is  based  on  an  assumed  interaction  between  un- 
cracked  concrete  and  a  deformed  bar  through  a  mechanical  system  that  uses 
the  same  approach  as  that  of  Yankelevsky  [65],  Hence,  it  assumes  the 
same  mode  of  load  transfer  between  the  reinforcement  and  the  matrix.  The 
model  applies  primarily  to  continuous  fibers  and  to  long  discontinuous 
fibers  where  it  can  be  observed  that  one  or  more  crack  will  cross  the  fiber 
so  that  the  load  is  applied  to  the  fiber. 


8,2  BASIC  ASSUMPTIONS 

a)  The  structural  model  used  by  Yankelevsky  for  a  reinforced  bar 
embedded  in  concrete  applies  to  fibers. 

b)  The  angle  a  that  is  used  in  the  derivation  is  constant. 

c)  The  fibers  are  asssumed  aligned  within  the  specimen,  and  are 
squarely  packed  (Fig.  7.2). 

d)  Each  individual  fiber,  along  with  its  share  of  matrix  acts  and 
behaves  independently  of  other  fibers  and  the  rest  of  the  matrix 
body.  The  representative  unit  is  further  simplified  to  a  cylinder 
of  matrix  encasing  one  fiber  (Fig.  7.3). 

e)  Finally,  the  tensile  load  acting  on  the  representative  unit  is  as¬ 
sumed  applied  directly  at  the  fiber  tip. 


8.3  MATHEMATICAL  DERIVATION  OF  THE  SHEAR  STRESS 
DISTRIBUTION 

A  typical  representative  unit  of  the  specimen  is  analyzed.  This  unit 
consists  of  one  fiber  and  a  cylinder  of  matrix  that  encases  it,  even  though 
when  a  square  packing  is  assumed,  a  typical  unit  would  be  expected  to  be  a 
fiber  encased  by  a  prism  of  matrix  with  a  square  cross  section.  However,  a 
circular  cylinder  with  a  cross  sectional  area  equal  to  that  of  the  prism  is 
chosen  to  ensure  that  all  the  matrix  area  is  accounted  for.  The  length  of  the 
cylinder  would  be  C,  sarne  as  the  length  of  the  fiber  (see  Fig.  7.3).  Since  a 
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uniform  packing  is  assumed,  the  exterior  diameter  of  this  cylinder  D  can  be 
directly  found  from  the  fiber  diameter  d  and  the  volume  fraction  of  the 
fibers  Vj.  The  area  of  the  fiber  Af  and  that  of  the  matrix  are  interre¬ 
lated  with  the  volume  fraction,  thus: 

Vf  =  _ ^ _ 

f  Af  -»■  Am 


thus: 

Vf 


D2 


(8.1) 


A  tensile  force  Fg  is  applied  to  the  representative  unit.  This  force  is 
a  fraction  of  the  total  tensile  force  P  applied  to  the  specimen.  A  fraction  of 
Fg  is  applied  to  the  fiber  end,  while  the  balance  of  it  is  applied  to  the  ma¬ 
trix  directly.  One  extreme  case  would  be  to  assume  that  the  whole  force  Fg 
is  applied  at  the  fibers'  ends,  while  nothing  is  given  to  the  matrix  (Fig. 
8. 1-a).  Under  this  assumption,  the  tensile  force  in  the  fiber  is  maximum  at 
the  fiber's  ends,  decaying  to  a  minimum  at  the  middle.  Thus,  the  tensile 
force  in  the  matrix  would  be  zero  at  the  ends,  and  maximum  at  the  middle. 
The  bond  shear  stress  and  the  slip  at  the  interface  between  the  fiber  and  the 
matrix  would  be  maximum  at  the  ends,  and  zero  at  the  middle.  The  other 
extreme  case  would  be  to  assume  that  the  matrix  takes  the  force  at  the  ends 
in  its  entirety  (Fig.  8.1-c).  In  that  case,  the  force  in  the  fiber  would  be 
zero  at  the  ends,  reaching  a  maximum  in  the  middle,  while  the  tension  in 
the  matrix  would  decrease  from  a  maximum  at  the  ends  to  a  minimum  at  the 
middle.  The  bond  stress  and  the  slip  would  behave  in  a  similar  way  as  in 
the  first  case.  If  one  is  to  assume  that  the  force  Fg  is  shared  between  the 
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Fig.  8.1  -  Typical  force  distribution  in  fiber  F(x)  and  matrix  T(x), 

assuming  that  the  unit  force  Fj,: 

(a)  is  applied  at  the  fibers’  ends 

(b)  is  shared  by  the  matrix  and  the  fiber  in  proportion  to  their 
relative  stiffnesses 

(c)  is  taken  entirely  by  the  matrix  ends. 
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matrix  and  the  fiber  in  proportion  to  their  relative  stiffnesses,  then  the  ten¬ 
sion  in  the  fiber  would  be  constant  throughout  the  fiber  length,  and  so 
would  the  force  in  the  matrix  (Fig.  8.1-b).  The  bond  shear  stress  as  well 
as  the  slip  between  the  fiber  and  the  matrix  would  both  be  nil.  The  model 
presented  herein  is  derived  based  on  the  first  assumption,  i.e.  that  the  force 
Fq  is  directly  applied  to  the  fibers'  ends  (Fig.  8.2)  as  in  Yankelevsky' s 
model  [65]. 


Defining  the  x-axis  as  being  the  fiber  axis  with  the  origin  at  one  of 
the  ends,  say  the  left  one,  a  differential  element  of  length  dx  of  the  fiber, 
located  at  a  distance  x  from  the  origin  is  considered  (Fig.  8.3-a).  The 
tensile  force  in  the  fiber  is  resisted  by  inclined  compressive  forces  dC,  thus 
reducing  the  axial  force  F  in  the  fiber  as  x  increases.  The  structural  model 
through  which  the  resisting  force  dC  is  transferred  to  the  matrix  is  shown 
in  Fig.  8.3-b.  The  compressive  struts  AB  and  A’B'  make  an  angle  a  with 
the  x-axis,  and  react  on  the  fiber  by  a  force  dC  per  unit  length  of  fiber 
perimeter,  hence  creating  a  compressive  stress  dC/dx  on  the  fiber  pe¬ 
riphery  that  makes  an  angle  a  with  the  fiber  axis.  As  was  mentioned 
before,  the  angle  a  is  assumed  to  be  location  independent,  a  is  indeed  as¬ 
sumed  to  be  a  material  property  directly  related  to  the  interfacial  properties 
of  the  composite.  The  circumference  BB'  is  the  locus  of  segment  centroids 
of  a  typical  disc  EE'  (Fig.  8.3-b).  This  disk  would  have  an  external 
diameter  D  and  an  internal  diameter  d.  The  forces  dN,  dC,  and  dT  as 
shown  in  Fig.  8.3-b  are  in  equilibrium.  By  such  a  mechanism,  the  tensile 
force  is  reduced  at  a  rate  of  dF/dx  in  the  fiber,  and  is  increased  at  an  equal 
rate  in  the  matrix.  The  centroid  location  r  (Fig.  8.3-b)  is  equal  to: 


r 


1  + 


(8.2) 


but  can  be  usually  taken  as  D/3,  especially  for  lower  reinforcement  ratios 
or  higher  D/d  such  as  in  the  case  of  FRC. 

Forces  are  then  interrelated  through  static  equilibrium  and  compati¬ 
bility.  A  differential  equation  in  F  (the  axial  force  in  the  steel)  is  obtained 


I 


a.  Actual  load  on  b.  Assumed  load  on 

representative  unit  representative  unit 


Fig.  8.2 


Load  on  representative  unit. 


to  be  analyzed 

;ed  in  the  development  of  the  model, 
evsky. 
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from  compatibility  equations  (for  more  details,  see  Ref.  6). 
form  of  the  differential  equation  was  given  as: 

d^F 

Ai  -  A2  F  =  -  A3 


The  general 


(8.3) 


The  coefficients  Aj,  A2,  and  A3  have  been  modified  from  their  val¬ 
ues  given  in  [65]  to  accommodate  the  problem  at  hand  (fiber-reinforced 
concrete),  as  well  as  to  include  the  Poisson’s  ratio  effect.  These  coeffi¬ 
cients  are  found  to  be: 


sin^  a 


-  ^  2  q(-n,2  .  1) 

4  7t  cos2  a 


^  fi  + 

M  1  * 


(8.4) 

(8.5) 

(8.6) 


((1  +  q2)  V  (I  *  q2)]  [(1  +  q2m2)  -h  v  ( 1  -  q2m2)]  (8.7) 
1  1 


II 

m 

> 

E 

Atn 

1 

(l  + 

11 

Ef  Af 

1 

"e^)  ( 

1  + 

Vf 

1  -  Vf 


_4 _ \  /I  +  (n  -  l)Vf\ 

d2  EfJ  I  1  -  Vf  J 


(8.8) 


where: 


modular  ratio 


(8.9) 
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V  =  Poisson  ratio  of  matrix 


in  which: 

Pc  =  Total  force  on  the  composite 

Ac  =  Cross-sectional  area  of  composite 

d  =  diameter  of  fiber 

The  solution  to  Eq.  8.3  is  of  the  form: 


F(x)  =  Cl  eP*  +  C2  c'P*  +  ^ 
where: 

P 

^  =  F  r _ _ "i  =  F* 

A2  0  (n-l)  Vf  -h  1  j 

P  can  also  be  expressed  in  the  following  form: 


(8.  10) 

(8.  11) 


(8.12) 

(8.13) 

(8.14) 

(8.15) 
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in  which: 


Zl  = 


sin 


M 

q  (m2  - 


Z2 


1  ^ 

n  ^  1  -  Vf 


(8.16) 

(8. 17) 


Factors  Zj  and  are  introduced  for  mathematical  convenience  and 
have  no  particular  physical  significance. 

*  •  •  • 
The  factor  is  the  product  of  the  overall  force  Fq  times  the  fiber's 

relative  stiffness  (Eq.  8. 14),  while  the  balance  of  Fq, 


portional  to  the  relative  stiffness  of  the  matrix. 


* 


-  F 


f 


1  -  V 


L 


(n-1)  Vf  +  1 


(8.18) 


To  coefficients  Cj  and  C2  are  obtained  from  the  following  boundary 
conditions: 

(1)  (2)  X  =  0,  F  =  Fq  (8. 19) 

Cl  C2  +  Fj.  =  Fq 

or  Cl  +  C2  =  (8.20) 

L  dF 

(2)  (2)  X  =  2-  =  0,  and  hence  ^  =  0.  (8.21) 

^  =  Cl  p  eP*  -  C2  |3  c-P*  (8.22) 

Cl  p  cP^/2  .  C2  P  c-P^/2  =  0 
Cl  =  C2e-P^ 


or 


(8.23) 
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(1)  +  (2)  =>  C.  e-PU  Co  =  F* 


therefore: 


F(x)  = 


and  T(x)  = 


X(x)  = 


m  1  +  e-3^ 


m  1  +  e-pr 

p*  f  — ^  ePx  +  - 1 — 

1  +  e-3^  1  +  e- 


1  +  e-P^ 


1  +  e-P^ 


-dF  *  _1_ 


(8.24) 


(8.25) 


J-Jj  e-Px  Uf;  (8.26, 


(8.27) 


-P  F 

^  /  e-P^ 

K  d  I  1  +  e-P^ 


eP*  - 


1  +  e'P^ 


P  Fo  I  -  Vf  e-P*  -  e-P^eP’^ 
j,  d  (n-1)  Vf  +  1  1  +  e-P^ 


(8.28) 


8.4  HBER  DEBONDING 


Eq.  8.28  gives  the  theoretical  bond  or  shear  stress  distribution  along 
the  fiber  interface.  However,  it  must  be  recognized  that  there  is  a  limit  that 
the  shear  stress  t  cannot  exceed,  namely  the  bond  shear  strength  of  the  in¬ 
terface,  tmax  (Fig.  8.4).  The  maximum  theoretical  bond  shear  stress, 
referred  to  as  if,  is  given  by  Eq.  8.28,  with  x  =  0. 


P  Fq  1  -  Vf  1  -  e-Pf 

7C  d  Vf  I  1  +  e-P^ 


(8.29) 


SHEAR  STR 


Fig.  8.4  -  Example  of  a  bond-slip  rclai 
exceed  the  bond  strength. 


-208- 


If  Ti  does  not  exceed  the  bond  distribution,  as  predicted  in  Eq. 

8.28  holds.  This  means  that  the  stress  applied  to  the  composite  does  not 
cause  any  fiber  debonding.  Debonding  will  occur  when  Tj  =  X[„ax-  ^^e 
force  Fq  that  would  cause  debonding  can  hence  be  found: 


^max 


(n-1)  Vf  +  1\ 
(1  -  Vf)  p  j 


(l  e'P^)  K  d 


(8.30) 


But  the  stress  in  the  composite  =  P/A^  can  be  related  to  Fq 

through  Eq.  8. 1 1.  The  applied  composite  stress  that  would  initiate  debond¬ 
ing  can  hence  be  found  as: 


.  ^  ^  X/  .  (n-1)  Vf  1  1  +  e-3f 

deb  ^  niax  ^  j  p  1  -  e'P^ 


(8.31) 


P  ^max  (8-32) 

In  Eq.  8.32,  p  is  a  dimensionless  factor  that  is  indicative  of  the 
strength  of  the  bond  within  the  composite.  Here,  p  will  be  termed  the 
Debonding  Index.  Fig.  8.5  shows  the  variation  of  this  factor  with  the 
volume  fraction,  for  different  values  of  a.  These  were  all  chosen  less  than 
45*  in  recognition  of  the  fact  that  p  is  almost  perfectly  symmetrical  about  a 
=  45*;  in  other  words,  p  will  have  almost  the  same  value  for  a  =  45’  -i-  aj 
and  a  =  45*  -  aj,  when  ai  ^  45*.  It  is  clear  from  the  figure  that  for  a 
given  angle  a,  the  Debonding  Index  p  increases  almost  linearly  with  the 
volume  fraction  Vf.  This  curve  would,  nevertheless,  only  apply  to  fiber- 
reinforced  composites  in  which  the  fiber  volume  fraction  is  relatively  low, 
say  Vf  ^  5%.  For  higher  volume  fractions,  the  interaction  between  the 
fibers  comes  into  the  picture,  hence  changing  the  mechanics  of  the  prob¬ 
lem.  Furthermore,  everything  else  being  equal,  the  closer  a  is  to  45*,  the 
lower  the  value  of  p. 


Fig.  8.6  shows  the  variation  of  the  Debonding  Index  with  the 
modular  ratio  n,  for  different  values  of  the  volume  fraction.  As  expected, 
higher  values  of  n,  meaning  stronger  relative  fiber  moduli,  lead  to  higher 
factors.  This  means  that  everything  else  being  equal,  the  higher  the  mod¬ 
ulus  of  the  fiber,  the  higher  the  debonding  stress. 
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Fig.  8.5 


Fig.  8.6 


Debonding  index  versus  volume  fraction  for  different  values 
of  a. 


MODULAR  RATIO 


Debonding  index  versus  modular  ratio  for  different  values  of  • 

the  volume  fraction. 
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In  the  majority  of  cases,  the  factor  e'P^  can  be  neglected  in  Eq.  8.31, 
as  the  Debonding  Index  turned  out  to  be  basically  independent  of  the  fiber 
aspect  ratio. 


8.5  DEBONDING  LENGTH 

If  the  applied  stress  on  the  composite  is  less  than  (ctc^deb’ 

bond  distribution  along  the  interface  between  a  given  fiber  and  the  matrix  is 
as  given  in  Eq.  8.28.  On  the  other  hand,  if  the  composite  stress  exceeds 
(a(.)deb>  fiber-matrix  interface  can  be  divided  into  two  parts:  a 

debonded  part,  2C^  long,  at  both  ends  of  the  fiber,  and  the  other,  [  -  2[^^ 
long,  that  is  bonded.  The  higher  the  composite  stress  Oq,  the  larger  the 
value  of  the  debonded  length  2/d.  There  is  no  closed  form  solution  for  the 
general  case  where  0  <  /j  ^  U2.  To  find  the  value  of  half  the  debonded 
length  /"(j,  the  same  differential  equation  described  in  Eq.  8.3  has  to  be 
solved,  with  the  following  differences: 

a)  The  length  f  is  to  be  replaced  with  /  -  2C^ 

b)  The  force  Fg  is  to  be  replaced  with  Fg  -  t^ax  V 

The  interfacial  shear  stress  distribution  in  the  bonded  zone  becomes: 


'Cb(x)  = 


^  Z'  I  -  Vf  e-P*  -  e-Pf'^-^^^eP* 

TT  1  1  +  e-P^^V 


(8.33) 


where 

^g  “  '  ^max  V  (d 

For  a  given  P,  hence  the  length  of  the  debonded  zone  can  be 
found  by  solving  the  equation: 

“  ^max 


or 
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^max 


P  (Fq  -  -tmax  ¥  ^d)  f  1  -  Vf  1  - 

7t  d  I  (n-1)  Vf  +  1  1  +  e-P^^Vj 


(8.34) 

Therefore,  to  find  the  debonded  length  for  a  given  Fg  such  that 
(Fo)dgb^Fo^('^inax  the  non-linear  Eq.  8.34  has  to  be  solved  by  trial 

and  error. 


The  bond  distribution  along  the  interface  between  the  reinforcement 
and  the  matrix  in  the  case  where  debonding  occurs  is  as  follows: 

''(’t)  =  ■'max  for  0  <  X  <  Id 


and 


T(x) 


n  d 


-  Vf 
Vf 


e-px  .  e-P^^  VeP* 

1  +  e-P^^^^^ 


for  \  <  f/2 


(8.35) 

Fig.  8.7  shows  two  cases  of  loading:  in  one  case,  tj  is  less  than 
"^max’  hence  no  debonding  occurs,  while  in  the  other,  the  applied  stress 
is  larger  than  (^rc)(jeb’  hence  the  bond  distribution  at  the  interface  is  de¬ 
scribed  by  Eq.  8.34.  Fig.  8.8  shows  the  variation  of  the  interfacial  bond 
shear  stress  distribution  as  the  applied  composite  stress  increases. 

8.6  APPLICATION  OF  THE  MODEL 
8.6.1  Bond  Modulus 

The  slip  S  is  the  relative  displacement  of  the  reinforcement  with  re¬ 
spect  to  the  matrix,  or: 
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Fig.  8.8 


Fig.  8.9 


T 

mu 


Variation  of  bond  distribution  with  increasing  applied 
composite  stress. 


Bond  shear  stress  versus  slip  relationship  based  on  predicted 
bond  modulus  and  on  bond  strength. 
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S(x)  =  5n,  -  5f 

=  J  dx 


(8.36) 


/l 


dx 


^  -  F  1  dx 


1 

Ef  Af 


F  dx 


(8.37) 


Replacing  F  and  T  by  the  expressions  in  Eqs.  8.26  and  8.27,  and  re- 

«  « 

placing  Ff  and  F^  by  their  respective  values  given  in  Eqs.  8. 14  and  8. 15, 


it  can  be  shown  that  the  slip  S  can  be  expressed  as  follows: 

_  4  Fq  e'Px  .  e'P^ePx 

”  p  7C  d2  Ef  1  + 


(8.38) 


It  can  be  noticed  from  Eqs.  8.59  and  8.69  that  the  ratio  of  the  shear 
stress  to  the  slip  in  the  elastic  range  (t-to-S  ratio)  is  independent  of  x,  and 
hence  of  the  location.  Furthermore,  this  ratio  is  constant  given  the  volume 
fraction,  the  type  of  fibers  and  the  matrix.  This  ratio,  the  bond  modulus, 
is  equal  to; 


p2  (1  .  Vf)  d  Ef 
4  [1  -h  (n  -  1)  Vf] 


(8.39) 


Recalling  the  expression  of  P  given  in  Eq.  8.17,  Eq.  8.39  can  be 


rewritten  as: 
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Z2  can  be  replaced  by  its  value  as  given  in  Eq.  8.  17.  Eq.  8.40  hence 
becomes: 


K 


sin2(2a)  E^ 
Zi  d 


(8.41) 


The  bond  modulus  K  is  the  slope  of  the  bond  stress  versus  slip  curve 
in  its  elastic  ascending  portion  (Fig.  8.9). 

Em 

Since  the  ratio  has  the  same  units  as  K,  Kj  is  hence  a  dimen¬ 
sionless  factor  referred  to  as  the  Bond  Modulus  Index  (BMI).  This  factor 
increases  almost  proportionally  with  increasing  volume  fraction.  Figs. 
8. 10  and  8. 1 1  show  the  variation  of  iCj  as  a  function  of  the  volume  fraction 
for  different  values  of  the  angle  a  and  as  a  function  of  the  angle  a  for 
different  volume  fractions  Vf,  respectively.  As  can  be  seen  in  Fig.  8. 10, 
for  a  given  volume  fraction,  the  highest  K\  is  obtained  for  an  angle  a  = 
45*.  For  angles  less  than  or  greater  than  45’,  Kj  is  less.  The  bond  modu¬ 
lus  K  is  hence  directly  proportional  to  the  elastic  modulus  of  the  matrix  En,, 
inversely  proportional  to  the  fiber  diameter  d,  almost  directly  proportional 
to  sin2(2a),  of  the  volume  fraction  Vf,  but  is  almost  insensitive  to  the 
value  of  the  Poisson's  ratio  of  the  matrix  v,  even  though  v  is  included  in 
the  expression  of  Kj. 


The  bond  modulus  k  is  usually  thought  of  as  being  an  interface 
property.  Therefore,  its  dependency  on  the  fiber  volume  fraction  Vf  is 
somehow  unexpected.  Nevertheless,  this  dependency  can  be  accounted  for 
by  relating  it  to  the  confining  effect  of  fibers.  Indeed,  the  higher  the 
volume  fraction,  the  closer  the  fibers  and  the  smaller  the  diameter  D  of  the 
representative  unit. 


BOND  MODULUS  INDEX 
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VOLUME  FRACTION 

Fig.  8,10  -  Bond  modulus  index  versus  volume  fraction  for  different 
values  of  a. 
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Fig.  8.11  - 


Bond  modulus  index  versus  (x  for  different  values  of  the 
volume  fraction. 
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8.6.2  Composite  Pre-cracking  Strength 


The  model  presented  here  can  be  used  to  predict  the  pre-cracking 
strength  in  tension  of  fiber-reinforced  cement  composites,  assuming  the 
fibers  to  be  perfectly  aligned  and  squarely  packed.  A  typical  stress-elon¬ 
gation  relationship  for  a  fiber-reinforced  cementitious  composite  would 
display  an  almost  linear  behavior  up  to  cracking.  Beyond  this  point,  a 
different  behavior  is  observed.  The  maximum  pre-cracking  stress  or  the 
pre-cracking  strength  is  hence  significant  and  is  addressed  in  what  follows. 


From  Eq.  8.27,  the  maximum  tensile  force  in  the  matrix  can  be  found 
by  setting  x  =  j  .  Or: 


T 


max 


2  c-m  X 

1  e-3fj 


(8.43) 

(8.44) 


where: 


f2  c-m 

(^1  +  e‘3^  J 


(8.45) 


The  maximum  tensile  stress  in  the  matrix  is  hence  the  ratio  of  this 
force  to  the  area  of  the  matrix,  A,„,  that  is: 


^m 


Inui 

^m 


Q  F 


m 


‘m 


(8.46) 


where: 


-218- 


Am  =  ^  (D^  -  d2) 

7t  d2  /D2 

4  [d2 

71  d2  -  Vf^ 


(8  '7) 


(8.48) 


The  composite  is  assumed  to  crack  when  the  tensile  stress  in  the  ma- 
trix  Fj  reaches  the  tensile  strength  of  the  matrix,  ft-  Setting  equal  to  ft, 

is  then  found: 
m 


(8.49) 


but  F^  is  related  to  Fq  as  in  Eq.  8. 18,  therefore  the  following  Fq  is  found: 


Therefore,  the  cracking  strength  of  the  composite  is  the  product  of  a 
factor  times  the  tensile  strength  of  the  matrix. 

It  is  worth  noting  that  in  the  development  of  this  equation,  it  was  as¬ 
sumed  that  the  matrix  is  weaker  than  the  reinforcement.  Eq.  8.52  would 
therefore  only  pertain  to  composites  that  meet  this  condition. 
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A  closer  look  at  the  factor  Cl  ought  to  be  taken, 
following  can  be  written: 


2  sin(2a) 


^rom  Fq.  8.15,  the 

(8.53) 

(8.54) 


where 

0)  =  2  sin(2a) 

Cl  can  now  be  expressed  as  follows: 


(8.55) 


Cl 


1  - 


2 

1  +  Expf-di^ 


(8.56) 


CJ  is  thus  a  function  of  the  fiber  aspect  ratio  rather  than  the  fiber 

length  [.  It  is  also  function  of  O),  which  in  turns  depends  on  the  angle  a, 
on  the  volume  fraction  Vf,  on  the  modular  ratio  n,  and  on  the  Poisson  ratio 
of  the  matrix  v. 


The  factor  Q,  which  will  be  referred  to  as  the  Bonding  Factor,  re¬ 
flects  the  quality  of  bond  that  exists  between  the  matrix  and  the  reinforce¬ 
ment.  In  the  theoretical  case  of  a  "perfect  bond",  Q  can  be  taken  equal  to 
1,  whereas  where  there  is  no  bond  at  all  between  the  fiber  and  the  matrix  at 
the  interface,  0,  in  which  case  Eq.  8.52  does  not  apply.  In  the  case  of 
a  perfect  bond.  Cl  would  be  equal  to  I  and  the  cracking  strength  of  the  com¬ 
posite  would  then  be: 

(Oc)^,=  [(n  -  1)  Vf -*■  11  ft  (8.57) 

which  is  the  same  as  the  equation  that  could  be  predicted  from  a  trans¬ 
formed  section  analysis,  assuming  perfect  strain  compatibility,  linear  elas¬ 
tic  materials,  continuous  reinforcement,  as  well  as  a  matrix  weaker  than  the 
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reinforcement.  Fig.  8.12  shows  the  variation  of  the  bonding  factor  Q 
versus  the  fiber  aspect  ratio. 


(Oc), 


In  Fig.  8.13,  the  cracking  strength  ratio,  — 7-^  is  plotted  against  the 

‘t 

fiber  volume  fraction  Vf,  for  different  values  of  the  aspect  ratio.  It  can  be 
seen  from  the  graph  that  for  ^  >  60,  the  cracking  strength  ratio  becomes 


insensitive  to  the  bonding  factor  and  approaches  the  values  predicted  in  Eq. 
8.57.  It  can  also  be  seen  from  Fig.  8.13  that  everything  else  being  equal, 
composites  reinforced  with  lower  aspect  ratio  fibers  have  higher  cracking 
strengths.  This  could  be  explained  by  observing  that  a  shorter  fiber  would 
display  a  more  uniform  bond  shear  stress  distribution  at  the  interface.  The 
fibers  are  hence  more  efficient  as  far  as  bond  goes.  However  this  result  is 
surprizing,  since  the  theory  on  which  the  model  is  based  assumes  that  the 
load  is  applied  to  the  fibers  thus  is  applicable  primarily  to  long  fibers. 
While  the  theory  predicts  a  peak,  followed  by  a  dip  in  the  case  of  fibers 
with  aspect  ratio  =  10  (Fig.  8.13),  no  reasonable  explanation  can  be  pro¬ 
vided  for  such  a  behavior. 


Fig.  8.14  shows  that  the  cracking  strength  is  not  very  sensitive  to 
the  angle  a,  so  long  as  a  does  not  approach  the  extremes  values  of  0  or 
90*.  Finally,  the  cracking  strength  ratio  as  per  Eq.  8.57  is  plotted  in  Fig. 
8.15  versus  the  modular  ratio  for  different  values  of  the  volume  fraction. 
It  can  be  observed  that  the  cracking  strength  increases  with  both  the  volume 
fraction  of  fibers  and  the  modular  ratio. 


8.7  CONCLUDIND  REMARKS 

The  following  remarks  can  be  made  about  this  model: 

1.  This  model  predicts  the  bond  shear  stresses  at  the  interface  between 
the  fibers  and  the  surrounding  matrix  in  a  pure  tension  specimen.  It 
also  predicts  the  tension  stresses  in  the  fibers  and  in  the  matrix. 

2.  The  model  was  used  to  determine  theoretically  the  bond  modulus, 
which  is  the  slope  of  the  bond  shear  stress-slip  curve.  It  was  shown 
that  the  bond  modulus  is  location-independent. 


-221- 


oa 

O 

H 

u 

< 

U> 

U 

Z 

O 

z 

o 

CQ 


Fig.  8. 12  -  Bonding  factor  versus  aspect  ratio  for  different  values  of  the 
modular  ratio. 


Fig.  8.13- 


Cracking  strength  ratio  versus  volume  fraction  for  different 
values  of  the  aspect  ratio. 
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Pig.  8.14-  Cracking  strength  ratio  versus  a  for  different  values  of  the 
volume  fraction. 
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Pig-  8.15.  Cracking  strength  ratio  versus  modular  ratio  for  different 


values  of  the  volume  fraction. 
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3.  The  model  was  also  used  to  predict  the  pre-cracking  strength  of  the 
composite  as  a  function  of  the  tension  strength  of  the  matrix.  The 
expression  derived  for  the  pre-cracking  strength  is  based  on  the  as¬ 
sumption  that  the  matrix  is  weaker  in  tension  than  the  fiber. 

4.  An  expression  for  the  debonding  stress,  which  is  the  stress  that 
needs  to  be  applied  to  the  composite  to  impend  debonding,  was  also 
derived  from  the  model.  It  was  a  direct  function  of  the  strength  of 
the  bond  between  the  fiber  and  the  matrix. 

6.  While  the  mathematical  expression  for  the  bond  modulus  was  derived 
based  on  the  assumption  that  the  force  Fq  is  applied  in  its  entirety  to 
the  fibers'  ends,  it  can  be  shown  that  this  expression  (of  bond 
modulus)  is  actually  independent  of  this  assumption.  This  result 
came  as  no  surprise,  since  the  bond  modulus  should  be  a  material 
property  that  is  independent  of  the  load,  or  its  mode  of  application. 
As  for  the  other  expressions  derived, i.e.  bond  stress,  tension 
stress, precracking  strength,  they  were  dependent  on  the  load  applica¬ 
tion  assumption. 
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